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1  3  AB5TSAC  T 


Processes  that  can  produce  optical  damage  in  transparent  insulators  are  studied 
in  the  present  work.  A  general  discussion  is  given  of  laser  induced  damage  which 
distinguishes  between  breakdown  from  intrinsic  processes  and  damage  from  absorbing 
contaminants  called  inclusions.  In  order  to  measure  properties  of  intrinsic  breakdown 
it  is  nec  .ssary  to  either  avoid  inclusion  damage  completely  or  to  distinguish  its 
presence  when  it  occurs.  New  and  previously  published  techniques  for  identifying 
damage  from  inclusions  are  discussed  and  applied.  Because  self-focusing  can 
develop  and  confuse  bulk  damage  measurements,  a  detailed  summary  of  self-focusing 
theory  is  given  which  establishes  the  conditions  for  minimizing  beam  distortion  from 
s elf-foc using  effects.  A  series  of  controlled  experiments  are  then  described  which 
represent  the  first  unambiguous  measurements  of  an  intrinsic  damage  mechanism  at 
near-infrared  and  visible  frequencies.  The  intrinsic  damage  process  in  the  alkali 
halides  at  1.  06  pm  is  identified  as  an  electron  avalanche  that  is  virtually  in¬ 
distinguishable  from  dc  dielectric  breakdown.  Measurements  are  discussed  which 
probe  the  frequency  and  time  dependence  of  avalanche  breakdown  and  from  which 
details  of  the  avalanche  mechanism  can  be  inferred.  The  effect  of  crystal  disorder  on 
damage  fields  is  also  considered.  A  statistical  character  is  observed  in  the  intrinsic 
breakdown  process  and  compared  to  previously  published  results  of  surface  damage. 

A  study  of  the  relationship  between  surface  and  bulk  damage  is  shown  to  confirm  a 
model  by  Bloembergen  of  field  enhancement  at  structural  defects.  Finally,  a  summary 
of  avalanche  breakdown  theories  is  presented,  and  the  implications  of  the  experimental 
data  for  avalanche  theory  are  discussed. 
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ABSTRACT 


Processes  that  can  produce  optical  damage  in  transparent  insulators 
are  studied  in  the  present  Work.  A  general  discussion  is  given  of  laser 
induced  damage  which  distinguishes  between  breakdown  from  intrinsic 
processes  and  damage  from  absorbing  contaminants  called  inclusions.  In 
order  to  measure  properties  of  intrinsic  breakdown,  it  is  necessary  to 
either  avoid  inclusion  damage  completely  or  to  distinguish  its  presence 
when  it  occurs.  New  and  previously  published  techniques  for  identifying 
damage  from  inclusions  are  discussed  and  applied.  Because  self-focusing 
can  develop  and  confuse  bulk  damage  measurements,  a  detailed  summary 
of  self -focusing  theory  is  given  which  establishes  the  conditions  for 
minimizing  beam  distortion  frum  self-focusing  effects.  A  series  of 
controlled  experiments  are  then  described  which  represent  the  first 
unambiguous  measurements  of  an  mlrimA  damage  mechanism  at  near- 
infrared  and  visible  frequencies.  The  intrinsic  damage  process  in  the 
alkali  halides  at  1,06  pm  is  identified  as  an  electron  avalanche  that  is 


virtually  indistinguishable  from  dc  dielectric  breakdown.  Measurements 
are  discussed  which  probe  the  frequency  and  time  dependence  of  avalanche 
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breakdown  and  from  which  details  of  the  avalanche  mechanism  can  be 
inferred.  The  effect  of  crystal  disorder  on  damage  fields  is  also 
considered.  A  statistical  character  is  observed  in  the  intrinsic  break¬ 
down  process  and  compared  to  previously  published  results  of  surface 
damage.  A  study  of  the  relationship  between  surface  and  bulk  damage 
is  shown  to  confirm  a  model  by  Bloembergen  of  field  enhancement  at 
structural  defects.  Finally,  a  summary  of  avalanche  breakdown  theories 
is  presented,  and  the  implications  of  the  experimental  data  for  avalanche 
theory  are  discussed. 
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INTRODUCTION 


A.  The  Problem  and  Its  History 

Intense  laser  light  is  known  to  damage  initially  trans¬ 
parent  solids.  Such  damage  may  be  a  cumulative  material 
deterioration1  or  it  may  appear  as  a  sudden  catastrophic 
material  failure  accompanied  by  a  spark  and  permanent  material 
damage.  This  second  type  of  optical  damage,  which  is  the 
subject  of  the  present  work,  represents  a  serious  limitation 
on  the  design  and  operation  of  high-power  lasers.  In  addition, 
damage  to  laser  components  such  as  windows,  rods,  and  Pockels 
cells  requires  replacement  and  realignment  which  can  be  both 
expensive  and  mime-consuming. 

A  plethora  of  published  experiments  and  theories  have 
appeared  since  1964  with  the  aim  of  elucidating  the  mechanisms 
of  optical  damage  and  determining  damage  fields  for  initially 
transparent  solids.*^  The  early  experimental  work  in  the  area 
v/as  confused  by  a  number  of  problems i  uncertainties  in  the 

4 

laser  mode  structure,  the  presence  of  small  absorbing 
particles,^  self-focusing  in  the  bulk,^  and,  perhaps  most 
important,  the  lack  of  a  theoretical  framework  for  designing 
and  interpreting  damage  studies. 

Despite  these  problems  a  number  of  useful  conclusions 
were  reached  concerning  optical  damage.  First,  it  was  found 
that  the  most  frequent  cause  of  laser-induced  damage  was 
extrinsic,  a  dirt  effect  caused  by  the  heating  of  small 
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absorbine  inclusions  with  diameters  between  about  0.01  pm  and 
1  pm.  Second,  confirming  the  expectations  of  a  number  of 

7 

scientists,  it  was  found  that  in  damage  experiments  the 
collapse  of  high-power  laser  beams  from  the  well-known  self- 

g 

focusing  nonlinearity  was  increasing  the  light  intensities 

g 

inside  solids  to  values  in  excess  of  the  incident  intensities. 
Damage  developed  after  this  intensification  so  that  apparent 
damage  thresholds  in  the  bulk  of  materials  were  actually 
thresholds  for  catastrophic  self-focusing.  Because  techniques 
had  not  yet  been  developed  for  avoiding  inclusion  damage  and 
catastroohic  self-focusing,  experimental  measurements  prior  to 
1971  gave  virtually  no  information  on  the  nature  of  intrinsic 
damage  processes.  The  seemingly  unavoidable  self- focusing 
problem  compelled  many  damage  experimenters  to  focus  their 
attention  and  their  laser  beams  on  the  surfaces  of  transparent 
solids  where  self- focusing  could  not  develop.10  Work  commenced 
on  thin  film  damage11  where  numerous  materials  problems  added 
new  complexities. 

In  1971  two  important  experimental  studies  of  intrinsic 

12 

optical  damage  were  reported.  The  first,  by  Bass  and  Barrett, 

observed  statistics  at  1.06  pm  and,  in  a  later  work,  at 

0,69  pm1^  in  the  optical  strengths  of  surfaces  for  several 

materials.  These  statistics  could  not  be  explained  by  laser 

fluctuations,  but  they  were  explained  bv  a  simple  model  for 

electron  avalanche  breakdown. 

The  second  set  of  measurements,  performed  by  Yablonovitch 
1 4 

at  Harvard,  was  the  first  unambiguous  study  of  optical 
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intrlnsic  bulk  damage.  Yablonovitch  measured  the  damage 
proDerties  of  ten  alkali  halides  at  10.6  (in  and  found  that  the 
damaging  process  appeared  to  be  identical  to  dc  dielectric 
breakdown.  Self- focusing  was  avoided  by  making  use  of  an 
observation  already  oresont  in  the  literature,^1  namely  that 
catastrophic  self- focusing  cannot  occur  below  a  material- 
dependent  critical  power.  By  confining  the  laser  input  powers 
to  values  which  were  orders  of  magnitude  less  than  the  critical 
power  for  self-focusing  while  strongly  focusing  the  laser  beam 
with  external  optics,  Yablonovitch  was  able  to  achieve  the  high 
intensities  necessary  to  damage  his  crystals  without  the  con¬ 
fusing  presence  of  self- focusing.  In  designing  his  experiment, 
he  was  aided  by  the  fact  that  the  self- focusing  threshold 
Increases  with  wavelength  as  and  that  typical  Gaussian  beam 
diameters  for  the  COj  laser  are  about  an  order  of  magnitude 
larger  than  the  Gaussian  beam  diameters  for  YAGtNd  and  ruby 
lasers.  Inclusion  damage  was  distinguished  from  intrinsic 
breakdown  by  examination  of  the  damage  morphology. 

The  year  following  this  work,  Bloembergen  completed  a 
simple  but  intriguing  model1 ^  which  predicted  that  the  apparent 
optical  strength  of  surfaces  would  normally  be  less  than  the 
optical  strength  of  the  bulk.  He  recognized  that  the  di¬ 
electric  discontinuity  occurring  at  structural  defects  can 
enhance  the  electric  field  over  its  average  value.  Such 
defects  are  likely  to  be  created  at  surfaces  by  abrasive 

1  £  ^  n 

polishing.  Crisp,  et  al. ,  had  shown  earlier  that  surface 
field  enhancement,  in  this  case  from  Fresnel  reflection,  could 
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explatn  the  paradoxical  difference  between  the  damage  resis¬ 
tance  of  exit  and  entrance  surfaces  on  the  same  sample, 

A  number  of  important  experimental  questions  remainedi 
While  the  process  for  intrinsic  damage  at  10,6  ji  m  is  apparently 
the  same  as  at  dc,  namely  avalanche  ionization,  what  intrinsic 
damage  process  is  operative  in  the  near  infrared  and  visible? 
What  are  the  thresholds  for  this  intrinsic  process  in  various 
materials?  Are,  in  fact,  the  thresholds  well-defined?  What 
are  the  frequency  and  time  dependences  to  the  intrinsic  process? 
Is  the  damage  resistance  of  a  single  crystal  any  different  from 
that  of  a  polycrystal  or  an  amorphous  form  of  the  same  material? 
Is  Bloembergen* s  model  correct,  or  does  surface  absorption 
play  a  role  in  surface  damage?  In  short,  when  damaging  con¬ 
taminants  are  eliminated,  what  magnitudes  of  optical  fields 
can  initially  transparent  solids  withstand? 

3.  Organization 

The  present  work  addresses  itself  to  these  questions. 
Chapts.  i-3  contain  a  general  description  of  the  optical 
damage  process,  a  detailed  study  of  self-focusing,  and 
descriptions  of  the  experimental  arrangements  and  techniques. 
The  remaining  chapters  summarize  the  experimental  results  and 
theory  relevant  to  optical  electron  avalanche  breakdown. 

In  Chapt.  1,  the  results  of  a  simplified  model  of  in¬ 
clusion  damage  are  discussed.  Of  particular  importance  is  a 
discussion  of  the  pulse-width  dependence  to  inclusion  damage. 
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It  is  9hown  that  although  the  intensity  threshold  for  inclusion 
damage  will  always  increase  with  decreasing  laser  pulse  width, 
the  magnitude  of  this  increase  depends  on  the  inclusion  size. 
The  qualitative  frequency  dependences  of  the  two  important 
intrinsic  damage  processes,  electron  avalanche  breakdown  and 
multiphoton  ionization,  are  also  described  in  Chapt.  1.  Al¬ 
though  multiphoton  ionization  has  not  been  observed  in  the 
present  work,  a  brief  discussion  of  the  multiphoton  process 
is  given. 

Self- focusing  is  described  in  Chapt.  2  within  the  context 
of  the  paraxial  ray  approximation.  In  addition  to  clarifying 
the  conditions  under  which  catastrophic  self-focusing  can  be 
avoided,  the  discussion  also  includes  a  novel  technique  for 
measuring  self- focusing  parameters  using  intrinsic  optical 
damage. 

Chapt.  3  summarizes  the  experimental  arrangements  for  the 
measurements  of  intrinsic  optical  damage.  A  new  technique  for 
longitudinal  mode  selection  of  Q-switched  lasers  is  described, 
and  the  various  measurements  characterizing  the  two  types  of 
lasers  used  in  the  present  work  are  reported.  A  summary  is 
given  of  new  and  previously  published  techniques  for  avoiding 
and  identifying  inclusion  damage. 

A  series  of  controlled  experiments  on  intrinsic  optical 
damage  at  1,06  gm  and  0.69  gm  are  described  in  Chapts.  4-7. 

The  first  experiment,  conducted  at  1.06  gm  with  a  Q-switched 
laser,  extends  Yablonovitch ’ s  infrared  measurements  on  the 
alkali  halides  to  shorter  wavelengths.  It  is  found,  contrary 
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to  published  predictions,  that  electron  ava'Janche  breakdown 
at  3  x  lO1*1  herz  is  virtually  indistinguishable  from  dc  break¬ 
down  and  intrinsic  laser  breakdown  at  10.6  nm.  This  observation 
has  important  implications  relating  to  the  electron-phonon 
Interaction  in  high  fields.  A  series  of  time  related  obser¬ 
vations  give  additional  evidence  that  an  electron  avalanche 
Is  resoonsible  for  intrinsic  damage. 

Damage  measurements  at  visible  frequency  (0.69  gm)  are 
described  in  Chapt.  5  where  evidence  of  frequency  dispersion 
in  the  threshold  fields  for  avalanche  breakdown  is  Presented. 

The  dependence  of  the  damage  field  on  laser  pulse  duration, 
which  gives  information  on  the  temporal  development  of  the 
electron  avalanche,  has  been  studied  in  NaCl  at  1.06nm  using 
a  mode-locked  laser.  It  was  found  that  the  rms  optical  damage 
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field  increased  by  almost  an  order  of  magnitude  to  over  10 

volts/cm  as  the  laser  pulse  width  was  decreased  from  10  ns  to 

15  ps.  The  laser  data  is  used  to  calculate  a  field-dependent 

ionization  rate  that  is  found  to  agree  at  least  qualitatively 

with  estimates  of  the  dc  ionization  rate  as  determined  by 

1 8 

Yablonovitch  and  Bloembergeri.  Chapt.  5  concludes  with  an 
exDeriment  and  discussion  of  the  effects  of  crystal  disorder 
on  the  damage  field. 

Laser  damage  statistics  is  the  subject  of  Chapt.  6.  It 

is  found  that  the  bulk  damage  process  in  a  number  of  materials, 

including  NaF,  has  a  statistical  character.  Measurements  of 

the  distribution  of  avalanche  starting  times  is  compared  to 

1 2 

the  experiments  of  Bass  and  Barrett  on  surface  damage  and 
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found  to  be  consistent  with  those  experiments,  A  direct 
comparison  between  surface  and  bulk  damage  on  the  same  sample 
of  fused  quartz  shows  virtually  identical  statistics.  A  few 
speculative  ideas  are  presented  concerning:  the  source  of  the 
statistics . 

Chapt.  ?  summarizes  a  short  experiment  which  confirms 
Bloembergen' s  model  of  surface  field  enhancement.^  Three 
surface  finishing  techniques  are  evaluated  in  terms  of  their 
ability  to  produce  damage  resistant  surfaces. 

Chapt.  8  summarizes  the  uncertain  state  of  avalanche 
breakdown  theory.  The  implications  of  our  experimental  results 
for  avalanche  theory  are  discussed.  Some  suggestions  for 
future  work  in  oDtical  damage  are  made  in  Chapt.  9. 


CHAPTER  1 


LASER  INDUCED  DAMAGE  IN  SOLIDS 

A.  Generalized  Damage  Process 

Optical  damage  in  solids  occurs  when  an  intense  laser  beam  causes 
an  irreversible  alteration  of  the  medium.  This  alteration  can  be  a 
thermally-induced  fracture  without  a  phase  change,  or  more  commonly,  it 
is  a  phase  change  which  may  or  may  not  be  accompanied  by  thermally 
induced  fractures. 

The  damage  process  can  be  conceptually  divided  into  two  steps.  In 
the  first  step  heat  is  deposited  by  some  absorption  process  which  may  be 
extrinsic  in  origin  or  which  may  be  a  property  of  the  pare  crystal.  In  the 
second  step,  this  heat  deposition  leads  to  a  temperature  elevation  and  a 
material  irreversibility.  Damage  is  normally  only  di  rected  by  the  material 
irreversibility  so  that  the  observer  is  always  a  step  removed  from  the 
physics  of  the  heat  deposition.  A  complete  description  of  optical  damage, 
therefore,  must  in  principle  include  more  than  an  analysis  of  how  the 
material  gains  energy.  It  must  also  include  a  modeling  of  the  thermal 
problem  of  how  that  heat  produced  a  material  irreversibility.  Conditions 
under  which  this  second  step  can  be  ignored  in  analyzing  damage  are 
discussed  below. 

In  order  to  establish  approximate  minimum  requirements  for  the 
heat  deposition  process,  the  details  of  the  second  step  are  replaced  with  an 
explicit  model.  It  is  assumed  that  heat  is  deposited  adiabatically  and  that 
damage  occurs  when  the  solid  is  melted.  The  amount  of  heat  necessary  to 
incude  damage  can  then  be  approximated  in  a  simple  manner. 

According  to  the  law  of  Duiong  and  Pettit,  the  specific 
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heat  of  a  solid  at  room  temperature  is  of  the  order  of 
1  joule/em^~°K.  In  order  to  raise  the  temperature  of  a 
typical  transparent  insulator  to  its  melting  point  of  about 
1000°  K,  about  10^  joules/cm^  must  therefore  be  supplied. 
Melting  of  the  lattice  is  equivalent,  typically,  to  heating 
by  another  1000°  K  so  that  a  phase  change  in  an  insulator 
requires  on  the  order  of  10^  to  10^  joules/cm^.  If  this  heat 
deposition  is  accomplished  by  the  absorption  of  1  pm  light, 
approximately  1022  photons  must  be  absorbed  per  cm^  in  order 
for  damage  to  be  produced. 

Thl3  energy  density  requirement  can  also  be  expressed  in 
terms  of  an  absorption  coefficient  if  uniform  bulk  absorption 
is  assumed.  Because  the  laser  intensity  varies  with  propaga¬ 
tion  distance  z  as  IQ  exp  (-az),  it  is  easy  to  show  that  for 
a  typical  damage  intensity  of  1010  watts/cm2  (see  Chapt,  4) 
and  for  a  laser  pulse  width  of  10  n.s ,  damage  can  occur  only 
when  a  =  10  cm*1.  Such  a  large  absorption,  if  present  for 
weak  light  signals,  makes  the  solid  opaque. 

On  the  basis  of  this  simple  calculation,  we  can  reach  a 
useful  conclusion  concerning  the  relationship  of  absorption 
by  atomic  impurities  to  optical  damage.  The  number  of  im¬ 
purity  states  in  a  high-purity  insulator  is  normally  less  than 
1017  /cm\  If  each  impurity  absorbs  one  photon,  the  amount 
of  energy  per  unit  volume  removed  from  the  laser  beam  is  less 
than  O.Ol'fo  of  the  energy  density  necessary  for  damage.  If 
absorption  is  occurring  because  of  impurity  excitation  into 
a  bound  state,  therefore,  each  impurity  must  be  excited  and 
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must  subsequently  relax  nonradiatlvely  a  total  of  at  least 
lO^  times  during  the  transit  of  the  laser  pulse  If  damage  is 
to  occur.  Although  systems  can  perhaps  be  found  where  this 
requirement  can  be  met,  it  seems  safe  to  conclude  that  im¬ 
purity  excitation  into  a  bound  state  cannot  normally  produce 
optical  damage  in  initially  transparent  solids. 

As  an  application  of  this  observation,  consider  the  rare 
earth  dysprosium  in  the  crystalline  host  yttrium  aluminum 
garnet  (YAG).  Dysprosium  is  known  to  strongly  absorb  1.06  pm 
radiation  through  an  electronic  excitation  into  a  bound  state 
having  a  lifetime  longer  than  1  ns.1^  An  interesting  experi¬ 
ment  might  appear  to  be  a  measurement  of  the  damage  threshold 
of  YAG  at  1,06  pm  as  a  function  of  dysprosium  concentration. 

On  the  basis  of  the  present  discussion,  however,  we  can  predict 
that  the  dysprosium  concentration  will  have  very  little 
effect  on  the  breakdown  strength  of  YAG.  Under  the  intense 
irradiation  necessary  for  damage,  the  dysprosium  absorption 
will  simply  bleach. 

An  absorption  process  must  be  found  which  can  produce  an 
a  =  10  cm  1  with  reasonable  electron  densities.  Joule 
heating  of  quasi-free  electrons  is  such  a  process.  Relaxation¬ 
time  models  of  the  electron- phonon  interaction  predict  that 
the  rate  of  energy  deposition  per  unit  volume  is  given  by 

dW  ..  J/,  2  2. 

-gf  -  N  e  p  E  /(1+  u)  t  )  (1-1) 

where  N  is  the  ccnduction-band  population  density,  e  is  the 


-12- 

electronic  chares,  g  is  the  electron  mobility  in  the  presence 

of  a  field  E  with  frequency  co,  and  r  is  the  eloctron-phonon 

2  2 

collision  time.  Letting  <u  r  <<  1  and  inserting  the  mobility 

o 

value  for  NaCl  (10  cm  /volt-sec)  and  damage  field  value 

(10^  volts/cm),  we  find  that  a  10  ns  (nanosecond)  laser  pulse 

can  deposit  a  damaging  amount  of  heat  when  N  >  lO^-7  /cm'. 

The  excitation  of  tnis  density  of  electrons  will  involve  a 

negligible  loss  of  energy  to  the  laser  beam. 

Let  us  apply  this  description  to  intrinsic  damage  cro- 

cesses.  Intrinsic  breakdown  in  transparent  solids  can  be 

caused  in  principle  by  one  of  two  processes,  either  electron 

2  0 

avalanche  ionization  or  multiphoton  absorption.  The  funda¬ 
mental  difference  between  these  processes  is  not  in  the  manner 
in  which  heat  is  deposited.  In  both  processes  energy  is 
transferred  to  the  electron  population  through  joule  heating 
of  conduction-band  electrons,  and  the  lattice  is  heated 
through  the  mechanism  of  electron-phonon  collisions.  The 
processes  are  fundamentally  different  only  in  the  manner  in 
which  they  produce  a  high  density  of  conduction  band  electrons. 
Electron  avalanche  breakdown  is  a  chain  reaction  generation 
of  quasi-free  electrons  through  impact  ionization  of  hot 
electrons  into  the  conduction  band,  f^ultiphoton  ionization, 
on  the  other  hand,  is  the  direct  excitation  of  electrons 
across  the  band-gap  because  of  the  simultaneous  absorption 
of  more  than  one  photon.  It  is  the  limiting  case  of  frequency 
dependent  tunneling. 

Having  considered  some  general  aspects  of  heat  deposition, 
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we  now  return  to  a  discussion  of  the  second  stage  of  optical 
damage  that  establishes  the  link  between  absorption  and 
actual  material  damage. 

As  the  conduction-band  population  of  the  lattice  or  of 
an  absorbing  inclusion  is  being  heated  by  joule  heating, 
energy  is  very  rapidly  transferred  to  the  lattice.  The  re¬ 
sulting  local  temperature  rise  induces  heat  flow  which  is 
described  by  the  geometry-d eDendent  thermal  diffusion  equation. 
When  the  heated  region  is  very  small,  thermal  diffusion 
effects  can  be  extremely  important  in  describing  the  tempera¬ 
ture  rise  because  the  thermal  diffusion  time  scales  with 
radius  R  of  the  heated  region  as  l/R  .  Thermal  loss  in 
principle  raises  the  damage  threshold  ano  alters  the  depen¬ 
dence  of  damage  fields  on  laser  pulse  width. 

The  thermal  asDects  of  intrinsic  damage  are  difficult 
to  consider  because  R  becomes  a  complicated  function  of  time 
for  intrinsic  Drocesses.  For  a  highlv  nonlinear  intrinsic 
process  such  as  avalanche  breakdown,  however,  the  thermal 
aspects  of  damage  need  not  be  calculated  because  they  have 
little  or  no  measurable  effect  on  damage  threshold.  If  the 
focusing  condition  in  laser  damage)  or  elect: ode  design  (in 
dc  damage)  of  one  experiment  make  breakdown  less  likely  oy 
altering  the  processes  leading  to  material  damage  thereby 
requiring  a  higher  rate  of  energy  deposition,  this  higher  rate 
can  be  achieved  by  an  immeasurably  small  change  in  the  elec¬ 
tric  field  strength.  The  damage  fields  for  avalanche  break¬ 
down,  therefore,  will  depend  only  on  the  parameters  which 
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describe  the  time  evolution  of  the  conduction-band  population 
N(t)  (and  therefore  the  heat  deposition)  and  will  have  no 
measureable  connection  with  the  geometry  of  the  field  con¬ 
figuration  or  the  morphology  of  the  damage. 

Inclusion  damage,  on  the  other  hand,  is  not  highly  non¬ 
linear  so  that  the  details  of  the  thermal  process  can  be 
extremely  important  in  determining  the  dependence  of  damage 
thresholds  on  laser  pulse  width.  This  dependence  will  be 
discussed  in  some  detail  in  the  next  section. 

Two-  and  three-photon  ionization  is  another  possible 
source  of  intrinsic  damage.  It  is  expected  that  the  sensi¬ 
tivity  of  the  damage  threshold  to  laser  pulse  width  will  be 
intermediate  between  that  of  avalanche  breakdown  and  inclusion 
damage,  A  precise  calculation  of  thermal  effects  is  difficult 
because  th .*  size  of  the  heated  region  is  time  dependent  unless 

the  multiphoton  ionization  is  occurring  in  a  uniform  manner 

2 1 

inside  a  microscopic  inclusion.  Bliss  has  presented  the 
results  of  a  calculation  for  the  pulse-width  dependence  of 
N^-order  multiphoton  breakdown  which  ignores  these  thermal 
effects.  His  results  follow  from  Eq.  (1-1)  with  a  simple 
modeling  of  N(t)  and  should  be  regarded  as  approximate. 

In  summary,  then,  laser  damage  most  likely  occurs  from 
the  joule  heating  of  conduction  band  electrons  in  a  bulk 
solid  or  in  an  absorbing  inclusion.  This  absorption  produces 
a  temperature  rise  and  eventually  a  material  irreversibility 
that  actually  defines  the  onset  of  damage.  A  proper  des¬ 
cription  of  the  damage  process,  especially  of  its  pulse-width 
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dependence,  must  Ln  principle  include  a  discussion  of  the 
manner  in  which  the  heat  deposition  produces  the  material 
irrevers ibility. 

B.  Inclusion  Damage 

Optical  damage  can  be  caused  by  absorbing  microscopic 
contaminants  called  inclusions  if  the  temperature  of  an  in¬ 
clusion  can  be  elevated  sufficiently  to  induce  a  material 
irreversibility  before  an  Intrinsic  damage  mechanism  develops. 
As  already  noted,  inclusion  damage  is  geometry  dependent  be¬ 
cause  the  total  process  is  sensitive  to  the  manner  in  which 
the  neat  deposition  produces  a  temperature  rise.  We  will 
eventually  use  a  simplified  model  to  study  this  damage  process 
and  to  determine,  in  particular,  its  pulse-width  dependence. 

The  temperature  rise  experienced  by  metallic  spherical 
inclusions  under  Q-switched  irradiation  has  been  treated  in 
detail  by  Hopper  and  Uhlman. ^  Damage  was  assumed  to  occur 
when  a  fracturing  stress  was  thermally  induced.  Although  this 
type  of  material  irreversibility  is  not  physically  equivalent 
to  a  phase  change,  the  energy  deposition  requirements  are 
about  the  same.  It  was  found  that  since  the  surface-to- 
volume  ratio  of  spherical  inclusions  varies  as  l/R ,  the  tem¬ 
perature  rise  induced  by  a  fixed  laser  pulse  incident  onto  an 
opaque  inclusion  will  initially  increase  as  the  size  of  the 
particles  is  decreased.  For  very  small  particles,  however, 
heat  will  begin  to  diffuse  away  from  the  inclusion  during 
the  transit  of  the  laser  pulse.  Since  the  thermal  diffusion 


-16- 


2 

time  is  of  the  order  of  D/R  »  where  D  is  the  thermal  dif¬ 
fusivity  of  the  host  (assumed  to  be  much  greater  than  the 
thermal  diffusivity  of  the  inclusion),  the  temperature  rise 
as  a  function  of  R  will  be  peaked  when  D/r  is  on  the  order 
of  the  laser  pulse  width  t  .  This  behavior  is  illustrated  in 
Pig.  la  where  it  is  seen  that  inclusions  having  diameters 
<  0.1  nm  will  probably  not  induce  damage  when  irradiated  by 
a  Q-switched  laser  pulse  having  an  energy  density  of  20  joules/ 
cm  and  a  pulse  duration  of  30  ns. 

For  long  pulses  and  fixed  inclusion  size  the  damage 
threshold  should  be  dependent  on  the  laser  intensity  or, 
equivalently,  since  the  focusing  conditions  are  fixed,  on  the 
laser  power.  When  tD  is  less  than  the  thermal  diffusion  time, 
the  threshold  for  damage  will  depend  on  the  laser  pulse  energy 
so  that  the  damage  power  will  vary  as  tp”1.  These  simple 
predictions,  which  are  al  so  valid  for  semiconducting  inclusions, 
explain  the  experimental  data  of  Danlleiko  et  al.*°  (see 
Fig.  lb)  who  irradiated  the  surfaces  of  sapphire  samples  with 
laser  pulses  of  varying  time  duration. 

We  now  consider  a  simplified  model  of  inclusion  damage 
that  displays  the  qualitative  features  of  Hopper  and  Uhlman's 
more  accurate  calculation  for  metallic  inclusions.  This 
model  will  be  used  to  treat  semiconducting  inclusions  and  to 
account  in  an  approximate  manner  for  the  finite  wavelength  of 
the  incident  laser  light.  Spherical  inclusions  are  considered 
with  the  thermal  diffusivity  of  the  inclusions  assumed  to  be 
much  larger  than  that  of  the  host.  Damage  is  assumed  to 
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occur  when  the  temperature  of  the  surrounding  host  material 
reaches  a  critical  value  ATC,  if  the  material  damage  is  a 
phase  change,  then  this  critical  temperature  rise  will  be 
given  by  the  difference  between  room  temperature  and  the 
melting  point  with  an  approoriate  addition  to  account  for 
the  heat  of  fusion. 

In  the  absence  of  thermal  dissipation,  the  temperature 
gain  AT  of  a  metallic  inclusion  with  radius  R  subject  to  a 
square  laser  pulse  of  intensity  I  is  given  by 


C  at 

V 


(area)I 
X  (volume) 


(1-2) 


where  Cv  is  the  specific  heat,  is  the  spectral  emissivity, 
and  K  is  a  correction  factor  to  account  for  the  effects  of  the 
finite  wavelength  of  the  irradiating  light  on  the  absorption 
cross-section.  For  metallic  particles,  M  is  approximately 
given  by31,  21 


M  = 


1  when  R  *  x/20 

20R/X  when  R  <  x/20 


(1-3) 


Exact  calculations  of  M  are  very  difficult,  especially  for 
absorbing  semiconductors,  and  only  a  few  special  cases  have 
been  considered  in  detail.  We  will  use  the  simplified  form 
of  Eq.  (1-3)  to  describe  the  diffraction  effects  for  both 
semiconducting  and  metallic  inclusions. 

After  the  laser  pulse  has  passed,  the  inclusion  will 
cool  and  AT  will  decay  approximately  as 
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AT0  e 


-Dt/Rl 


(1-4) 


where  D  is  the  thermal  diffusivity  of  the  host.  Eqs.  (1-2) 
and  (1-4)  together  imply  that  the  time-dependent  temperature 
rise  of  a  metallic  inclusion  during  the  transit  of  a  square 
laser  pulse  is 


'X  IRM 


e.Dt/R 


(1-5) 


The  damage  intensity  Id  is  defined  as  the  intensity  of  a 


laser  pulse  of 

duration  tv 

r 

which 

induces  a  critical 

ture  rise  ATC 

.  From  Eq. 

(1-5)  , 

we  find 

T(m) 

Xd 

/ 4C  DaT 
/  v  c 

\  .1 

1 

'  3ex  M 

/  R 

1  -  exp(-  Dtp/RZ) 

(1-6) 


If  the  inclusions  are  made  of  a  contaminated  semicon¬ 
ductor,  they  can  still  absorb  largp  amounts  of  light  energy 
and  cause  optical  damage.  Residual  inclusive  in  high  puri”y 
optical  crystals  are  normally  less  than  a  micro  in  diameter 
so  that  the  absorption  depth  for  a  semiconducting  inclusion 
will  normally  be  greater  than  the  particle  diameter.  For  such 
inclusions  in  Eq,  (1-2)  is  replaced  by  a  term  proportional 
to  aF  where  a  is  the  frequency-dependent  absorption  constant. 
The  temperature  rise  and  damage  intensity  then  become 


aT(s,«  ir2m  (1  -  e'Dt/R  j 


(1-7) 
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and 


AT 


_ 1 _ 

1  -  exp  (-5T7?) 

P 


(1-8) 


Eq.  (1-5)  has  the  same  quallt  clve  form  as  the  more 
accurate  solution  for  metallic  inclusions'’  as  displayed  in 
Fig.  la.  An  optimum  radius  Rm  exists  for  metallic  particles 
in  the  sense  that  any  other  value  of  R  will  require  a  higher 
light  intensity  in  order  to  produce  the  same  temperature  rise 

AT.  R  can  be  shown  to  be 
m 

Rm=  0.9(Dyl/2  (1.9) 


when  M  =  1  for  the  simple  model  we  have  considered.  Heating 

of  semiconducting  inclusions  produces  a  different  behavior  as 

seen  by  Eq .  (1-7).  When  aR  <<  1  as  we  have  assumed,  the 

( s ) 

temperature  rise  AT'  increases  monotonically  with  increasing 
radius  so  that  an  optimum  radius  for  semiconducting  inclusions 
does  not  exist. 

Eqs.  (1-6)  and  (1-8)  have  the  same  limiting  behavior  as 
the  experimental  data  of  Danileiko  et  al .  as  summarized  in 
Fig.  lb.  When  Dt^R2  <<  1,  the  damage  threshold  is  defined 
by  the  energy  in  the  pulse  reaching  a  critical  value.  I .  e.  . 


R/M,  metallic  inclusions 
l/M,  semiconducting  inclusions. 


For  long  pulse  durations  (or  small  particles)  such  that 
nt  /r2  >>  1,  the  damage  threshold  is  indeneno »nt  of  the  laser 
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pulse  duration.  In  particular, 


Xd  “ 


1/MR,  metallic  Inclusions 
1/MR2;  semi  conducting'  inclusions. 


(1-11) 


Sincb  the  laser  input  power  is  proportional  to  the  laser 
intensity,  the  pulse-width  dependence  to  the  damage  power 
displayed  in  Fig.  lb  is  to  be  expected  if  damage  is  occurring 
from  absorbing  inclusions  which  all  have  the  same  radius  R. 

In  \\n  actual  damage  experiment  inclusions  present  in  the 
focal  volume  may  be  different  sizes  so  that  the  experimental 
results  of\  Danileiko  et  al.  represent  a  special  case.  Let  us 
consider  th\e  more  general  problem  An  experiment  is  per¬ 
formed  in  whVch  the  threshold  for  damage  is  determined  in  a 
single  sample,,  for  two  different  pulse  durations.  Using 
methods  to  be  \liscussed  in  Chapt,  3,  optical  damage  is  shown 
to  result  from  inclusion  absorption.  How  does  the  damage 
intensity  compare  for  the  two  pulse  durations?  It  is  assumed 
either  that  large\  focal  volumes  are  irradiated  so  that  each 

focal  volume  is  effectively  the  same  or  that  a  large  number 

\ 

of  sites  are  probed  and  meaningful  averages  obtained. 

The  solution  to  this  problem  depends  on  the  average  dis¬ 
tribution  of  inclusion  sizes  in  the  focal  volume.  We  will 
study  a  few  special  cases. 

Consider  first  a  solid  in  which  the  probability  of 
finding  a  metallic  inclusion  with  R  >  X/20  -  0.05  pm  is  unity 
for  all  R.  If  D  =  3  x  10"2  cm2/sec  (appropriate  for  NaCl), 
then  an  optimum  radius  R  can  he  found  from  Eq.  (1-9) 
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whenever  tp  >  1  ns.  By  Inserting  Eq.  (1-9)  into  (1-6),  we 
can  show  that  for  such  pulse  durations  the  damage  intensity 
1^  will  scale  with  laser  pulse  duration  as 


V 


1/2 


(1-12) 


Eq.  (1-12)  predicts  a  behavior  which  is  different  from  either 
Eq.  (1-10)  or  (l-ll). 

When  semiconducting  inclusions  are  present  in  the  focal 
volume  with  the  same  distribution  of  particle  sizes  and  the 
same  pulse-width  range  assumed  for  Eq.  (1-12),  the  damage 
intensity  is  proportional  to  t^-1  as  in  Eq.  (1-10).  This 
result  assumes  that  Dt^R  <<  1,  a  condition  that  is  valid 
when  the  largest  particle  radius  present  is  1  jim  and  tp<< 

330  ns.  The  reason  that  the  pulse-width  dependence  is 
different  for  the  two  types  of  Inclusions  is  that  the  damage 
intensity  has  a  minimum  as  a  function  of  R  only  for  metallic 
inclusions.  When  semiconducting  inclusions  are  present,  the 
damage  intensity  will  decrease  with  increasing  R  until 
aR  -  1. 

Consider  now  an  experiment  similar  to  that  of  Sect.  C 
in  Chapt.  5  in  which  damage  is  induced  with  laser  pulses 
having  durations  of  15  ns  and  15  ps  (picoseconds).  We  wish 
to  determine  how  the  thresholds  for  inclusion  damage  compare 
for  the  two  pulse  durations.  Assume  first  that  all  inclusion 
sizes  with  R  <  1  pm  are  present.  According  to  Sq.  (1-8)  the 
damage  intensity  for  semiconducting  particles  will  increase 
by  a  factor  of  10^  as  the  pulse  width  is  decreased  from  1 5  ns 
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to  1 5  ps.  The  situation  is  somewhat  more  complicated  when 
metallic  inclusions  are  present.  If  M  *  1,  an  optimum  radius 
Rm  can  be  found  as  before  and  the  relative  damage  fields 
calculated  from  Eqs.  (1-6)  and  (1-9).  When  t  <  1  ns,  however, 
Eq.  (1-9)  no  longer  applies.  For  such  short  pulse  durations, 
diffraction  effects  as  summarized  by  the  parameter  M  in  Eq. 
(1-3)  introduce  an  additional  factor  of  R  in  the  expression 
for  the  damage  field  Id^m'.  With  this  change,  I^m^  has  the 
same  radius  dependence  as  the  damage  field  for  large  semi¬ 
conducting  particles.  The  inclusion  size  which  is  most 
potentially  damaging  when  t  <  1  ns,  therefore,  is  R  =  A/20. 
Using  this  fact,  the  damage  field  is  found  from  Eqs.  (1-6) 
and  (1-9)  to  increase  by  a  factor  of  about  150  when  the  laser 
pulse  duration  is  decreased  from  15  ns  to  15  p3.  This  increase 
is  almost  an  order  of  magnitude  less  than  the  corresponding 
increase  for  semiconducting  inclusions. 

When  only  particles  with  fixed  radius  are  found  in  the 
focal  diameter,  the  damage  intensity  for  both  types  of  in¬ 
clusions  scales  with  pulse  duration  as 

Id(tj)  1  -  expl-Dty'R^)  (1-13) 

If  we  let  tj  =  15  ns  and  tg  =  15  ps ,  the  calculated  ratio  of 
damage  intensities  for  particles  with  R  =  1  gm  is  10-*,  and 
for  particles  with  R  =  0.01  pm,  the  ratio  is  2.8, 

It  is  thus  seen  that  the  change  in  damage  intensity  with 
decreasing  pulse  duration  depends  strongly  on  the  distribution 
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of  particle  sizes  in  the  focal  volume.  Under  certain  con¬ 
ditions,  damage  from  semiconducting  inclusions  can  have  a 
different  apparent  pulse-width  dependence  than  damage  from 
metallic  inclusions. 

Particles  as  small  as  0.01  gm  may  not  cause  optical 
damage  in  an  actual  experiment  with  Q-switched  laser  pulses, 
because  the  damage  intensity,  especially  for  semiconducting 
inclusions,  may  exceed  the  intrinsic  breakdown  strength  of 
the  host,1 ^  The  intrinsic  breakdown  field  increases  for 
subnanosecond  laser  pulses,  however,  as  shown  inSect.  C  of 
Chapt.  5.  As  the  laser  pulse  width  is  decreased  from  about 
10  ns  to  15  ps,  the  intrinsic  damage  intensity  increases  by 
a  factor  of  about  35  in  NaCl,  and  it  may  be  possible  to 
observe  inclusion  damage  from  small  particles  with  a  15  ps 
mode-locked  laser  pulse  when  damage  produced  in  the  same 
sample  by  a  a-switched  laser  is  observed  to  proceed  by  an 
intrinsic  mechanism.  Bliss  and  Milam,11  in  fact,  have  ob¬ 
tained  evidence  for  this  presumably  rare  situation  in  their 
studies  of  ootical  damage  to  thin  film  optics,  and  we  have 
found  during  the  subnanosecond  damage  study  of  Chapt.  5  that 
inclusion  damage  can  occur  more  frequently  compared  to  in¬ 
trinsic  damage  when  15  ps  light  pulses  are  used  to  induce 
damage  than  when  5  ns  light  pulses  are  used.  This  situation 
is  unusual,  however.  It  is  normally  found  that  inclusion 
damage  can  be  avoided  by  shortening  the  laser  pulse  duration. 

The  model  of  inclusion  damage  that  we  have  applied  in 
this  section  is  phenomenological  and  should  be  considered 
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to  be  a  qualitative  description.  It  was  presented  largely 
to  demonstrate  the  manner  in  which  the  second  step  of  the 
damage  process--the  step  that  describes  the  sample’s  tempera¬ 
ture  elevation — affects  the  pulse-width  dependence  of  inclusion 
damage.  This  dependence,  which  for  spherical  inclusions  can 
be  translated  into  an  inclusion-size  dependent  damage  inten¬ 
sity,  can  be  treated  accurately  by  using  a  thermal  diffusion 
equation.  Any  prediction  of  the  pulse-width  dependence  must 
Include  an  assumption  about  the  distribution  of  inclusion 
sizes. 

The  frequency  dependence  to  inclusion  absorption,  which 
is  determined  by  the  details  of  the  heat  deposition  (step  1), 
has  been  ignored  to  this  point  in  the  discussion.  This  de¬ 
pendence  occurs  because  the  absorption  cross-section  changes 
with  frequency  due  to  the  frequency  dependences  of  joule 
heating,  impurity  ionization,  band-to-band  absorption,  optical 
scattering,  etc.  Even  without  attempting  to  quantify  this 
frequency  dependence,  it  is  clear  that  normally  inclusions 
will  become  more  absorbing  with  increasing  laser  frequency 
in  the  visible  and  UV.  This  simple  observation  suggests  that 
inclusion  damage  may  become  more  of  a  practical  problem  as 
high-power  lasers  with  shorter  wavelengths  become  available. 

C .  Intrinsic  Damage  Processes 

Avalanche  breakdown  and  multiphoton  ionization  were 
introduced  in  Sect,  A,  In  this  section  we  will  briefly  re¬ 
consider  these  intrinsic  breakdown  mechanisms  from  the 
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viewpoint  of  their  qualitative  frequency  dependences.  This 

discussion  will  be  extended  somewhat  in  Chapt.  5  for  avalanche 

breakdown.  The  laser  pulse-width  dependence  of  avalanche 

breakdown  will  also  be  considered  in  Chapt.  5* 

Crosc-sections  for  multiphoton  ionization  have  been 

22 

calculated  from  perturbation  theory.  It  is  found  that  is 
the  laser  frequency  increases,  the  cross-section  for  the 
ionization  of  valence-band  electrons  will  increase  in  a  com¬ 
plicated  fashion.  Successively  lower-order  multiphoton 
processes  become  possible  as  the  frequency  is  raised  until 
single- photon  absorption  or  direct  band-to-band  excitation 
becomes  possible.  At  low  frequencies,  where  only  high-crder 
multiphoton  processes  can  occur,  the  excitation  process 

changes  in  a  continuous  manner  into  Zener  tunneling  as  ‘-ho wr. 
20 

by  Keldysh.  It  is  expected,  then,  that  multiphotor,  absorp¬ 
tion  will  have  the  behavior  shown  schematically  in  Fig.  ?, 
Quantitative  calculations  of  multiphoton  damage  inten¬ 
sities  cannot  be  accurately  obtained  from  perturbation  them*.  »« 
Not  only  do  the  various  matrix  elements  become  large  in 
typical  damage  fields  of  10^  volts/cm,  but  the  electron1,  c 
levels  within  a  solid  become  distorted.  When  bound  srareo 
are  involved  in  the  ionization  process,  say  from  impuri  t 
the  distortion  is  particularly  serious.  This  distortion 
develops  because  the  state  is  life-time  broadened  ou«'  to 

large  probability  of  transitions  into  the  conduction  V-ani  or 

2 1 

into  other  states  J  and  because  the  Stark  interaction  ms*.;:, 
3tates  and  thereby  increases  the  energy  uncertain1- v . 


A  SCHEMATIC  COMPARISON  OF  THE  DAMAGE  THRESHOLDS  FOR 
AVALANCHE  BREAKDOWN  AND  MULTIPHOTON  IONIZATION  AS  A 
FUNCTION  OF  LASER  WAVELENGTH 


Q10HS3dHl  NMOOXV3d8 


At  low  frequencies,  avalanche  breakdown  normally  has  a  lower  threshold  than  multiphoton 
ionization  and  will  thus  be  the  cause  of  intrinsic  electric  damage.  As  the  frequency  is  raised 
however,  breakdown  from  multiphoton  ionization  rapidly  becomes  more  probable  and  eventua 
dominates  avalanche  breakdown. 
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Electron  avalanche  ionization  involves  an  energy  ex¬ 
change  between  the  electric  field  and  the  electrons  described 
qualitatively  by  the  joule  heating  formula  of  Eq.  (1-1).  In 
the  absence  of  losses  from  electron  diffusion  and  trapping, 
joule  heating  completely  describes  the  energy  input  into  the 
conduction  band  population.  Eq,  (1-1)  shows  that  this  energy 
input  scales  with  frequency  as  E2/( l+a'2 r2 ) ,  and  because  the 
details  of  the  energy  input  should  determine  the  electron 
distribution  function  and  hence  N(t),  the  threshold  for  damage 
will  scale  in  the  same  manner.  More  sophisticated  calculations 
of  breakdown  include  an  enersry  dependence  in  the  collision 
time  t  and  average  Eq.  (l-l)  over  the  high-field  electron 
distribution. 

It  is  shown  in  Chapt3.  4  and  5  that  avalanche  breakdown 
in  the  alkali  halides  appears  to  have  an  effective  collision 
time  on  the  order  of  10  1  sec.  w  r  will  be  much  less  than 
1,  therefore,  until  frequencies  near  the  visible  spectrum. 

When  the  laser  frequency  becomes  comparable  to  the  collision 
frequency,  the  damage  field  should  begin  to  increase  as  shown 
schematically  in  Fig.  2, 

On  the  basis  of  this  qualitative  discussion,  we  can  pre¬ 
dict  the  relative  importance  of  the  two  intrinsic  damage 
mechanisms  as  a  function  of  laser  frequency.  While  it  will  be 
shown  in  Chants.  4  and  5  that  avalanche  breakdown  is  the  most 
likely  cause  of  optical  damage  in  transparent  insulators  near 
1  gm,  breakdown  from  multiphoton  ionization  will  become  in¬ 
creasingly  more  important  at  high  frequencies  and  should 
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eventually  dominate  avalanche  breakdown  as  the  frequency  is 
raised*  Fig*  2  summarizes  the  relative  importance  of  these 
two  processes  as  a  function  of  laser  wavelength.  Unambiguous 
experimental  confirmation  has  not  yet  been  made  of  the  com¬ 
petition  of  avalanche  breakdown  and  multiphoton  ionization  so 
that  future  damage  studies  in  the  visible  and  UV  are  highly 
desireable. 


CHAPTER  2 

SELF-FOCUSING  AND  ITS  RELATIONSHIP 
TO  OPTICAL  DAMAGE 

A.  Introduction 

A  laser  beam  propagating  in  a  transparent  medium  induces  an 

increase  in  the  index  of  refraction  by  an  amount  proportional  to  the  laser 
6  7 

intensity.  ’  At  powers  in  excess  of  a  critical  power,  this  nonlinearity 
causes  the  intensity  distribution  to  become  unstable  and  a  catastropic 
beam  collapse  occurs  in  a  finite  distance  z^..  In  solids  the  nonlinearity 
normally  results  from  electrostriction  or  from  a  nonlinear  electronic 
poiarizability. 

A  physical  description  of  this  process  is  summarized  in  Fig.  3 
where  the  propagation  of  a  high-power  laser  beam  having  a  Gaussian 
intensity  distribution  is  followed.  Diffraction  effects  are  ignored.  Along 
the  beam,  axis  the  light  intensity  is  highest  so  that  the  index  of  refraction 
is  greater  on  axis  than  it  is  in  the  wings.  Because  the  propagation 
velocity  is  given  by  c/n  ,  where  n  is  the  total  index  of  refraction  in¬ 
cluding  the  nonlinear  term  (see  Fig.  3),  light  in  the  wings  of  the  beam  will 
move  faster  than  light  along  the  axis,  thereby  causing  initially  plane  phase 

fronts  to  become  concave  and  the  beam  to  focus.  Beam  collapse  continues 

18  2  5 

until  some  nonlinear  process  limits  the  intensification.  ’  The  manner 
in  which  diffraction  effects  modify  this  description  will  be  discussed  in 

detail  below. 

Self  -  focusing  is  important  to  studies  of  optical  damage, 
because  unless  great  care  is  taken  in  the  design  of  damage 
experiments,  catastrophic  self  -  focusing  will  occur  and  cause 
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the  input  laser  team  to  intensify.  The  high  intensity  that 
results  from  self- focusing  will  then  Induce  damage,  but  this 
high  intensity  may  have  no  known  relationship  to  the  incident 
intensity.  When  catastrophic  self-focusing  is  present,  there¬ 
fore,  apparent  damage  thresholds  will  actually  be  thresholds 
for  catastrophic  self- focusing. 

In  this  chapter  the  effects  of  the  self-focusing  non¬ 
linearity  on  the  propagation  of  high  power  laser  beams  will  be 
discussed  within  the  framework  of  a  paraxial  ray  approximation. 
This  discussion  provides  the  theoretical  basis  for  the  design 
of  the  bulk  damage  measurements.  Some  important  results  of 
this  discussion  arei 

1.  The  tendency  of  light  to  diffract  and  the  tendency  to 
self-focus  act  in  opposite  directions  in  the  sense  that 
diffraction  effects  try  to  bend  light  rays  away  from  the 
axis  of  propagation  and  cause  initially  plane  phase  fronts 
to  become  convex  while  self- focusing  effects  try  to  bend 
light  rays  towards  the  axis  and  cause  initially  plane 
Phase  fronts  to  become  concave. 

2.  The  relative  balance  between  self- focusing  and  diffraction 
does  not  change  as  the  beam  propagates  through  the  non¬ 
linear  medium.  As  a  result,  if  diffraction  dominates  at 
the  entrance  plane,  it  will  dominate  everywhere. 

3.  There  exists  a  critical  power  Pc  such  that  when  the  laser 
input  power  P  is  equal  to  Pc,  the  focusing  effects  from  the 
nonlinear  index  approximately  balance  diffraction  effects 
over  the  entire  laser  beam  and  an  initially  colummated 
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laser  beam  propagates  with  essentially  no  change  in  size. 

When  P  >  Pc,  the  laser  beam  will  collapse  after  propagating 
some  finite  distance  inside  the  nonlinear  medium.  When 
P  <  Pc,  catastrophic  self-focusing  is  absent  but  beam 
distortion  from  the  index  nonlinearity  is  still  present. 

4.  A  second  power,  Pcr,  i-S  more  useful  in  describing  this 
beam  distortion  than  is  Pc.  Pcr  is  defined  as  that  value 
of  input  power  for  which  diffraction  and  self-focusing 
effects  balance  on  axis.  For  Gaussian  beams  its  value 

is  0.273  Pc.26 

5.  Transient  effects  make  self-focusing  more  difficult. 

6.  When  a  laser  beam  with  P  <  Pcr  is  focused  inside  a  non¬ 
linear  medium,  the  focal  diameter  will  be  smaller  than 
its  diffraction  limit  even  though  catastrophic  self- 
focusing  is  absent.  Because  the  change  in  beam  diameter 
can  be  calculated  to  high  precision  from  the  paraxial  ray 
approximation  when  the  index  nonlinearity  is  in  its  steady 
state,  self- focusing  parameters  can  be  calculated  as 
demonstrated  in  Sect.  F. 

B.  General  Formalism  for  a  Steady-State  Nonlinearity 

The  effect  of  the  nonlinear  index  of  refraction  on  the 
propagation  of  light  beams  can  in  principle  be  determined 
exactly  from  Maxwell’s  equations.  In  practice,  however,  an 
analytical  solution  is  impossible  unless  a  number  of  simplifying 
assumptions  are  made.  One  such  set  of  assumptions  leads  to  a 
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paraxial  ray  approximation  that  describes  the  qualitative 
behavior  of  self- focusing  at  all  power  levels  and,  in  addition, 
results  in  a  good  quantitative  fit  to  numerical  calculations 

27 

at  low  laser  powers  where  catastrophic  self-focusing  is  absent. 
Since  we  are  primarily  concerned  with  such  low  powers,  we  will 
use  the  paraxial  ray  approximation  to  describe  self-focusing. 

A  paraxial  ray  equation  for  a  nonlinear  medium  was  first 
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derived  by  Talanov  and  was  used  in  detail  by  Wane  and 

Wagner  et  al.^°  to  describe  catastrophic  self-focusing. 

7 

Zverev  and  Pashkov  applied  this  formalism  to  a  focused,  low- 
power  laser  beam.  The  applicability  of  the  paraxial  ray 
formalism  is  restricted  to  laser  beams  which  have  ray  slopes 
much  less  than  one  and  which  have  diameters  2a  much  larger 
than  a  wavelength--conditions  fulfilled  in  our  damage  experi¬ 
ments. 

In  order  to  describe  the  effects  of  the  self-focusing 
nonlinearity  on  light  propagation,  we  begin  with  Maxwell’s 
equations.  Wang2<^  has  shown  that  when  the  electric  field 
vector  E(r)  is  written  as  A(r)  exp  i  lk<£(r)  -  co  t]  ,  where 

2 

kQ  =  w/c  and  the  index  of  refraction  is  written  as  n0  +  n2  A  /2, 
Maxwell's  equations  reduce  to  the  simplified  form 

[<nQ  +  n^A  /Z)  -  (grad  6)  ]  A  +  — =j-  v  A  =  0. 

k 

o 

2  2 

Terms  of  order  (n2  A  )  have  been  ignored.  When  the  scalar 

product  of  this  equation  is  taken  with  A,  a  term  containing 

/  4  2 

the  factor  i/k0  dropped,  and  a  cross  term  with  n2  A 

i 

neglected,  this  result  leads  directly  to  an  effective  eikonal 
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equation 


2  2 
nl  -  (grad  t)  - 


where 


nl  =  no  +  2  n2  a2  4- 


1 


v2a 

„  2  - 

2k  n  A 
o  o 


(2-1) 


(2-2) 


In  the  limits  of  zero  nonlinearity  and  infinitesimal  wave¬ 
length,  n1  =  nQ  and  Eq.  (2-1)  is  just  the  basic  equation  of 
geometrical  optics.  Results  derived  from  the  usual  eikonal 
equation  can  be  used  for  the  nonlinear  problem  if  the  index  of 
refraction  is  replaced  by  Eq.  (2-2).  In  particular,  it  can  be 
shown^1  that  the  curvature  d2r/d p2  for  a  pencil  of  rays  with 
position  vector  r  and  with  p  the  coordinate  along  the  ray  path 
is  given  by 


d 


2-» 

r 


nl 


grad  n1  - 


dr 

dp 


(2-3) 


Eq.  (2-3)  can  be  simplified  by  applying  several  approximations. 

We  restrict  the  treatment  to  cylindrically  symmetric  beams  and 

assume  that  the  maximum  ray  slope  is  small  compared  to  unity. 

(In  our  experiments  the  maximum  slope  inside  the  sample  and 

before  the  focus  is  less  than  0.05.)  Both  the  second  term  on 

the  right  in  Eq.  (2-3)  and  the  longitudinal  component  of  the 

Laplacian  in  Eq ,  (2-2)  are  negligible.  The  curvature  is  now 

2  8 

expressed  in  a  form  first  derived  by  Talanov. 


d2r 

dz2 


grad 


n  i 


n , 


(2-4) 
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This  result  is  important  to  the  study  of  self- focusing 
effects  because  the  sign  of  d  r/dz*  Indicates  whether  or  not 
the  beam  Is  converging  and  its  magnitude  is  a  quantitative 
measure  of  that  convergence  or  divergence.  A  positive 
curvature  results  in  an  increase  in  the  slope  of  the  ray  path 
with  resoect  to  the  propagation  direction  and  thus  represents 
a  divergence  from  the  axis.  Diffraction  alone  will  produce  a 
positive  curvature  in  an  isotropic  medium.  A  negative 
curvature,  on  the  other  hand,  will  cause  convergence  of  the 
beam  towards  the  axis  and  indicates  the  dominance  of  the  self- 
focusing  nonlinearity. 

Eq.  (2-4)  can  be  written  in  a  useful  form  by  using  a 

2.6)  P 

"constant  shape  approximation"  '  which  we  now  describe.  A  (r,z) 
's  prooortional  to  the  light  intensity,  and  where  the  beam 
propagates  with  little  or  no  change  in  shape,  the  intensity  at 
r  and  z  is  equal  to  a  function  that  describes  the  beam  shape 
divided  by  the  beam  area  rra  (z).  The  shape  function  also 
depends  on  r  and  z,  but  in  a  very  simple  manner.  We  denote  it 
by  p.  This  function  may  be  Gaus3ian,  for  example,  with 
p  =  exp  (-2r  /a  ) .  If  only  the  beam  radius  changes  with  i 
and  the  beam  shape  does  not  change,  then  p  is  actually  a 
function  of  a  reduced  radial  coordinate  and  is  independent 
of  beam  size.  In  particular,  we  can  define  the  reduced  co- 
ordinate  x  as  x  =  r/a(z)  and  write  A  (r,z)  in  the  simplified 
form 


3l 

a 


A 


2  2 
A  (r,  z)  -  p(x)/n  a  (z) 


(?-5) 


a 


1 
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The  approximation  of  Eq.  (2-5)  is  normally  used  in  numerical 
calculations.  Along  with  Eq.  (2-4)  it  provides  the  basic 
relationship  needed  to  calculate  the  critical  power,  the  self- 
focusing  length,  and  the  quantitative  influence  of  the  index 
nonlinearity  when  diffraction  dominates. 

By  expanding  in  cylindrical  coordinates  r  and  z  it  is 

o  2 

easily  seen  that  A/A  *s  a^so  proportional  to  a  .  We  can 

then  write  Eq.  (2-2)  as 


n 


1 


n  + 
o 


a2(z) 


fix) 


(2-6) 


where  f(x)  is  the  sum  of  contributions  from  both  diffraction 
and  the  index  nonlinearity.  Gradx  is  just  a(z)-1  cWx  .  and 
n^  in  the  denominator  of  Eq .  (2-4)  can  be  replaced  by  nQ, 

Eq.  (2-4)  thus  gives  us  the  basic  equation  for  light  propaga¬ 
tion  within  the  paraxial  ray  approximation,  namely 


1  1 

3,  , 

a  (z) 
o 


—  lf(x)] 
dx 


(2-7) 


The  importance  of  neglected  terms  has  been  shown  to  be 
negligible29  except  near  a  catastrophic  self-focus  or  under 
experimental  condition  where  extreme  external  focusing  is  used. 

We  thus  have  our  first  important  conclusion  about  self- 
focusingj  the  derivative  of  f(x),  which  contains  no  dependence 
on  z,  determines  the  relative  importance  of  diffraction  and 
self- focusing.  Once  this  relative  importance  is  determined 
for  one  value  of  z,  such  as  z  =  0  at  the  entrance  plane,  it 
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is  determined  for  all  z,  The  critical  power  Pc  is  in  principle 
calculated  from  the  requirement  that  the  derivative  f(x) 
vanish  for  all  x,  leading  to  a  detailed  balance  of  self- 
focusing  and  diffraction.  Since  f(x)  contains  no  dependence 
on  beam  diameter,  the  critical  power  will  not  be  dependent  on 
beam  diameter  for  a  steady-state  nonlinearity. 


C,  Effects  of  Nonlinearity  on  Focused.  Low-Power  Laser  Beams 


The  index  nonlinearity  leads  to  intensity  distortions  even 
below  the  critical  power  for  catastrophic  self- focusing.  Using 
Eq.  (2-7),  quantitative  corrections  from  the  nonlinearity  can 
be  derived. 

The  coordinate  r  can  be  replaced  by  x  a(z)  on  the  left 
side  of  Eq.  (2-7)  to  obtain  a  differential  equation  for  the 
beam  radius  a(z),  namely 


Wagner  et  al 
near  x  =  0, 


d2a 


dz 

30 


a3(z) 


1  d  ,,  , 

- —  f(x) 


n  x 
o 


dx 


have  evaluated  Eq.  (2-8) 
They  obtain,  to  first  order 


(2-8) 

for  a  Gaussian  beam 


^4  =  [l/kVtz)j  U  -  P/P  ) 
dz 


(2-9) 


where  P  is  the  total  power  in  the  beam,  k  =  2n/\,  and  Pcr  is 
given  by 

P  =  cV2  (2-10) 

cr  32nn2 
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When  Eq.  (2-9)  is  multiplied  by  z(da/dz),  it  can  be  integrated 
subject  to  the  boundary  conditions  that  (da/dz)2  _  0  »  RQ 
and  a(0)  =  RQ.  When  the  laser  beam  is  focused  sufficiently 
far  inside  the  material,  R0  is  the  tangent  of  the  far-field 
diffraction  angle  whose  exact  value  will  be  given  later.  We 
find 

/da\ 2  “-P/Pc,'  .  ,2  ,  »-p/Pcr> 

\di  =  -  ~tt—  +  Ro  +  — -r~r  (2-u) 

K  a  (z)  k  R 

o 

The  minimum  beam  radius  ( 1 , e.  .  the  radius  of  the  focal 
point)  is  found  by  setting  da/dz  =  0  and  solving  the  resulting 
equation  for  a^.  This  orocedure  eivee  the  focal  radius  a-^  as 
a  function  of  input  power  and  self-focusing  parameters.  That 
is , 

a  --  ll-P/F  )1/2[k2[R2  +  (1-P/P  )/k2R2]  f1/2 

1  cr  1  1  o  cr  '  o  1  (?-'\  ?) 


For  the  power  range  of  interest  to  our  damage  experiments 
(P  <  pcr^»  Eq.  (2-12)  can  be  considerably  simplified  because, 
as  we  now  demonstrate,  the  power  dependent  term  in  the  brackets 
can  be  neglected.  We  want  to  show  that 


■  ?  ,  2  7 

R  +  l/k  R  >  > 

o  o 


P/P 


cr 


2  2 


(2-13) 


To  simplify  the  demonstration  that  this  equality  is  valid, 
the  positive  term  l/k2  RQ2  on  the  left  of  Eq.  (2-13)  is  ignored 
and  P/Pcr,  which  is  less  than  unity  in  the  region  of  interest, 
is  replaced  by  1.  Eq .  (2-13)  thus  reduces  to 
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lRJ>>1/kRo  •  u-,6) 

The  rate  of  change  of  the  Gaussian  beam  radius  at  the  entrance 
plane  Is  given  by 


I R  I  =  tan  p 

1  o 


I  =  o 0  for  email  angles  (weak  focusing} 

)  >  p  0  for  large  angles  (  s  t  r  c  *  c  u.  ti  *  j  * 


Assuming  the  worse  and  letting  |Rq  =  0o  where 


eo  =  -  <'/?-V 


for  a  Gaussian  beam,  we  find  that  the  inequality 
becomes 

(2/rr)  (\/2  a  )  >  >  \/2n  R 
o  '  o 


or 


2  R  >  ■>  a 
o  o 


Since  in  our  experiments  the  diffraction-limited  Gaussian 
beam  radius  a0  is  at  least  an  order  of  magnitude  smaller  then 
the  incident  Gaussian  radius  R0,  this  inequality  }  u  eas.liy 
satisfied. 

Because  we  can  neglect  the  power  dependent  term  in  the 
bracket  of  Eq.  (2-12),  that  equation  takes  the  Pintle  form 


*i  =  "  -p/pcr|i/2n 
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Eq .  (2-15)  summarizes  the  "behavior  of  a  laser  beam  with 
P  <  ?cr  that  is  focused  into  a  nonlinear  medium.  The  radius 
of  the  focus  is  reduced  below  its  diffraction  limit  aQ  because 
of  the  self-focusing  nonlinearity.  This  behavior  is  illus¬ 
trated  in  Fis:.  k.  As  the  laser  power  P  approaches  Pcr,  the 
focal  radius  decreases  until  it  becomes  so  small  that  the 
aporoximations  used  to  derive  Eq.  (2-15)  are  no  longer  valid. 
Eq.  (2-15)  appears  to  be  a  useful  approximation  for 
P  5  0.9  Pcr.  In  any  nonlinear  medium,  therefore,  the  laser 
intensity  on-axis  is  increased  above  its  diffraction  limit 
by  the  factor  (aj/a^j)  where  a0,  the  diffraction-limited  focal 
radius,  is  the  same  for  ll^ht  focused  in  air  or  in  a  solid. 

Pcr  can  be  shown  to  be  the  laser  input  power  for  which 

diffraction  and  self-focusine  effects  balance  at  the  center  of 

a  Gaussian  lisht  beam.  For  InDut  powers  greater  than  Pcr 

but  less  than  the  critical  Dower  Pc ,  diffraction  dominates 

•  27 

everywhere  except  near  the  beam  axis.  Numerical  calculations 
show  that  a  collimated  beam  will  initially  intensify  at  3uch 
powers  until  the  diffraction  of  the  winas  forces  the  on-axis 
intensity  to  drop.  Pcr  differs  from  Pc  because  the  latter  is 
a  quantity  averaged  over  the  entire  beam  while  Pcr  is  deter¬ 
mined  by  the  behavior  near  the  center.  In  fact,  Pc  is  not  a 
precisely  defined  quantity  because  it  is  not  possible  to 
exactly  balance  diffraction  and  self-focusinc  over  the  entire 
beam  cross-section.  At  an  input  power  of  Pc,  therefore,  a 
oropaaatin?  beam  will  not  change  its  size  measurably  and  so 
not  experience  a  catastrophic  self-focus,  but  its  intensity 
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Fig.  4.  THE  EFFECT  OF  THE  INDEX  NONLINEARITY  ON  THE 
MINIMUM  DIAMETER  OF  A  FOCUSED  LASER  BEAM 
WHEN  P  <  P 

cr 


PBN-72-128? 


Light  ray  paths  with  n2  s 0 

Light  ray  paths  with  * 0 

and  P  .<  P, 

input  cr 


The  tendency  to  self-focus  acts  in  opposition  to  diffraction.  Whereas  diffraction 
produces  a  positive  ray  curvature,  self-focusing  produces  a  negative  ray  curvature. 
The  focal  diameter  d --  2a  inside  a  nonlinear  medium  (n2=0)  will  be  effectively 
reduced,  therefore,  from  dQ--its  value  v/ith  02  ~  0--to  dj,  and  the  on-axis  light 
intensity  at  the  focus  will  be  correspondingly  increased.  As  shown  by  Wagner  et  al. 
(see  Ref.  30),  the  presence  of  the  index  nonlinearity  also  causes  the  point  of 
minimum  beam  radius  to  be  displaced  to  the  right  toward  the  geometrical  focus. 

This  latter  observation  is  net  important  to  the  discussion  in  the  text. 
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distribution  will  be  distorted.  Pc  has  the  same  functional 
form  as  Pcr  and  differs  by  just  a  numerical  factor  as  Pcr  = 
0.273  Pc  for  Gaussian  beams. 

D.  Transient  Effects 

The  above  analysis  is  correct  only  in  the  steady  state. 

In  solids  the  dominant  nonlinearity  is  normally  electro- 
strictlon,^  and  if  the  process  is  transient,  it  is  no  longer 
true  that  the  relative  balance  between  diffraction  and  self- 
focusing  is  independent  of  propagation  distance  and  that  the 
critical  powers  are  independent  of  beam  diameter.  The  changes 
occur  because  electrostriction  becomes  non-local  in  both  a 
temporal  and  a  spatial  sense.  Although  a  susceptibility 
approach  such  as  used  here  is  no  longer  strictly  correct,  it 
is  nonetheless  useful  for  establishing  functional  dependences 
for  self- focusing  parameters  and  approximate  quantitative 
values . 

Two  results  from  a  transient  analysis  are  important.^2 
The  first  is  that  transient  effects  decrease  the  effective 
nonlinear  index  n£  and  thus  make  self-focusing  more  difficult. 
The  second  important  result  involves  the  dependence  of  the 
critical  power  on  laser  pulse  width  and  on  beam  diameter.  In 
the  steady  state,  the  nonlinear  index  nj  is  given  by^ 
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where  nG  i9  the  index  of  refraction  in  the  absence  of  the 
nonlinearity,  p  is  the  material  density,  and  v  is  the  acousti¬ 
cal  sound  velocity.  The  quantity  (p  <Pr\0/<9p)  for  cubic 
materials  such  as  the  alkali-halides  may  be  found  approxl- 

-a*j 

mately  by  differentiating  the  Clausius-Mosottl  equation.  J 
The  remaining  constants  are  tabulated  in  handbooks. 


When  the  laser  pulse  width  tp  is  shorter  than  the 
electrostrictive  response  time  tp  =  a/v,  ^  is  decreased  in 
value,  thereby  increasing  the  critical  power.  For  a  triangular 
pulse,  Kerr^2  has  shown  that 


2*transient  ^n2^steady-state  ^  vp  ^  a  ^  (2-l6) 


where  D(-^)  is  Dawson's  integral  with 

cx 


D(^)  -  exp(-§)  j  eXp  rj 2  d  p 

0 


When  tp  1  tp/2 


2  transient 


k(n,)  .  ,  t  , 

2  steady-state 


This  result  is  valid  for  more  general  and  realistic  pulse- 
shapes  with  the  numerical  constant  k  being  of  order  unity  and 
having  a  value  dependent  on  the  precise  time-structure  of  the 
pulse.  When  this  result  is  inserted  into  Eq .  (2-10),  the 
critical  power  becomes 
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For  short  laser  pulses,  therefore,  when  the  process  is 
transient,  there  exists  more  properly  a  critical  intensity 
rather  than  a  critical  cower,  and  fairly  small  changes  in 
pulse-width  will  have  a  significant  effect  on  the  critical 
power. 

E.  Applications  of  Self-Focusing  Theory  to  Optical  Damage 

When  catastrophic  self- focusing:  occurs  in  solids,  the 
results  are  dramatic — an  extended  spark  and  local  disruption 
of  the  crystal.  Because  such  optical  damage  apparently 
always  accompanies  catastrophic  self- focusing  in  crystals,  it 
has  become  customary  in  the  damage  literature  to  associate 
damage  with  total  collapse  of  the  laser  beam.”^  Results  of 
the  preceding  sections  show,  however,  that  it  is  possible  to 
induce  optical  damage  without  measurable  beam  distortion  from 
self- focusing.  This  is  possible  because  catastrophic  self- 
focusing  and  beam  distortion  depend  on  the  laser  incident 
power  whereas  intrinsic  optical  damage  depends  on  intensity 
for  a  fixed  laser  pulse  width. 

These  observations  form  the  basis  for  the  design  of  our 
optical  damage  experiments.  A  laser  bean  with  input  power 
much  less  than  ■  £<?r  —  focused  tightly  with  external  optics  in 
order  to  achieve  the  high  intensities  necessary  to  induce 
intrinsic  damage.  Great  care  is  taken  to  avoid  measurable 
aberrations  in  the  focusing  optics. 

In  order  to  apply  this  technique,  estimates  of  Pcr  are 
required.  The  nonlinear  index  from  electrostriction  can  be 


-47- 


calculated  in  the  manner  discussed  in  the  previous  section, 
and  the  electronic  self- focusing  parameters  can  be  calculated 
from  the  third-order  electronic  susceptibility."^  In  parti¬ 
cular, 


where  -the  real  part  of  the  electronic  nonlinear 

( ' 

susceptibility  x  can  be  obtained  from  measurements 

of  third-harmonic  generation. Typically,  Pc  from  steady- 
state  electrostriction  at  1.06  pm  is  on  the  order  of  1  MW 
(106  watts)  while  the  electronic  critical  power  is  about  an 
order  of  magnitude  larger.  Because  of  transient  effects, 
electrostriction  will  weaken  as  the  laser  pulse  width  is  low¬ 
ered  below  about  5  ns  for  a  50  ^m  beam  diameter  and  will 
generally  be  unimportant  for  subnanosecond  laser  pulses. 

Experimental  checks  are  necessary  to  confirm  the  absence 

of  significant  beam  distortion.  Previous  experimenters  had 

q  37 

connected  the  damage  morphology  with  self-focusing,7’  J  but 
morphology  is  not  a  reliable  indicator  of  the  presence  of  self- 
focusing.  Intrinsic  damaee  in  the  absence  of  self-focusing 
effects  can  have  a  morphology  similar  to  that  seen  with 
optical  damage  preceded  by  catastrophic  self- focusing.  (See 
Chapt.  7.) 

The  simple  test  used  in  the  uresent  work  to  confirm  the 
absence  of  self-focusing  consists  of  measuring  the  damage 
field  with  lenses  having  different  focal  lengths.  If  the  input 
power  required  to  induce  damage  scales  with  the  square  of  the 
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focal  length,  as  expected  for  an  intensity  dependent  process 
such  as  intrinsic  breakdown,  then  catastrophic  self-focusing 
from  a  steady-state  nonlinearity  is  absent.  In  addition, 
using  Eq.  (2-15)  we  can  show  that  significant  beam  distortion 
is  also  absent. 

If,  on  the  other  hand,  the  laser  pulse  width  is  less  than 
the  electrostrictive  response  time  (a/v) ,  Eqs.  (2-15)  and 
(2-17)  predict  just  the  results  obtained.  It  may  be  desired 
to  demonstrate  the  absence  of  transient  electrostrictive  self- 
focusing  by  changing  the  laser  pulse  width  and  comparing  the 
input  damage  powers  with  those  expected  from  Eq.  (2-17). 

Because  intrinsic  breakdown  is  weakly  pulse-width  dependent 
as  shown  in  Chant.  5,  care  must  be  used  in  interpreting 
results. 

F.  The  Measurement  of  Self-Focusing  Parameters  Using  Intrinsic 
Optical  Damage 

Experimental  values  of  ng  for  various  materials  ha^e  been 

q  a 

determined  from  measurements  of  induced  birefringence,  by 

Q  TO 

studies  of  the  movement  or  positions  of  optical  damage 
tracks  in  crystals,  or  by  observing  the  changes  in  the  intensity 
distribution  of  the  transmitted  laser  beam  as  the  input  power 
is  increased. 2+0 

A  new  technique  for  measuring  self-focusing  parameters 
was  developed  from  the  results  of  Sects.  B-D.  It  consists  of 
determining  the  powers  necessary  to  induce  intrinsic  bulk 
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damage  as  lenses  with  different  focal  lengths  are  used  to 
focus  the  radiation  and  then  comparing;  these  results  to  pre¬ 
dictions  based  on  Eq.  (2-12)  •  High  input  powers  are  not 
needed  so  that  relatively  small,  stable  lasers  can  be  used. 

This  technique  is  described  in  detail  below,  and  as  an  illus¬ 
tration  of  its  application,  the  self- focusing  parameters  for 
sapphire  are  determined.  Some  experimental  results  and  ideas 
developed  in  later  chapters  will  be  briefly  introduced. 

Fig.  4  illustrates  schematically  the  effect  of  the  index 
nonlinearity  on  a  TEM00  laser  beam  when  P  <  Pcr  <  Pc.  As 
shown  by  Eq.  (2-15),  the  focused  beam  radius  is  decreased  from 
its  diffraction  limit  of  aQ  to  a  smaller  value,  a.,  because  of 
the  nonlinear  index.  As  a  result,  the  on-axis  intensity  is 
increased  by  the  factor  (a^/aj)2.  If  I0  is  the  intensity 
when  n2  =  0,  then  from  Eq .  (2-15)  we  find  that  the  intensity 
Id  in  the  presence  of  the  nonlinear  index  is  given  by 

*d  .  lo/(l  -  P/Pcr)  (2-19) 

where  P  is  the  input  Dower  of  the  laser  beam  and  P  in  esu 

cr 

units  is  given  by  Eq.  (2-10).  As  the  input  power  is  increased 

to  a  value  within  a  few  percent  of  Pcr,  corrections  to  Eq. 

(2-19)  due  to  diffraction  of  the  wings  of  the  laser  beam  may 
26 

become  important  and  numerical  calculations  may  be  required 
to  accurately  establish  the  relationship  between  Id  and  IQ. 

Although  we  have  considered  only  Gaussian  beams,  Eq.  (2-19) 
is  valid  for  any  smoothly  varying  intensity  distribution  which 
has  its  most  intense  region  along  the  axis  and  which  is 
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quadratic  in  r  near  the  axis.  The  precise  value  of  Pcr, 
however,  depends  on  the  coefficient  of  the  quadratic  term.  It 
was  shown  in  Sect.  C  that  corrections  to  Eq ,  (2-19)  are  unim¬ 
portant  when  P  <  Pcr  and  the  beam  is  focused  sufficiently  to 
induce  intrinsic  damage. 

The  steady-state  analysis  summarized  by  Eq,  ( >—19)  is 
appropriate  for  tightly  focused  laser  beams,  because  the 
electronic  contribution  to  n2  is  in  its  steady  state  for  all 
attainable  laser  pulses-^®  and  electrostriction  is  in  its 
steady  state  whenever  the  laser  pulse  width  tp  is  greater  than 
a  response  time  t^  eiven  by  the  ratio  of  the  radius  of  the 
focused  laser  beam  and  the  radial  compression  wave  velocity 
v.6  (In  sapphire,  where  v  =  1.1  x  10^  cm/sec  for  light  pro¬ 
pagation  along  the  crystalline  c  axis,  tp/lR «  1  when  tD  = 

Q 

10"  sec  and  the  focused  beam  diameter  is  220  inn.)  If  electro¬ 
striction  is  sliahtly  transient,  Eq.  ( 2 - i 9 )  may  still  be  used 
with  Pcr  corrected  in  the  manner  discussed  in  Sect,  D. 

Having  thus  determined  the  intensification  produced  by 
the  index  nonlinearity,  we  look  for  a  process  which  is  sensitive 
to  a  constant  value  of  on-axis  laser  intensity.  The  powers 
necessary  to  reach  this  intensity  will  be  recorded  as  lenses 
with  different  focal  lengths  are  used.  The  basic  requirements 
for  this  Drocess  are,  first,  that  it  be  easily  observable  with 
a  reasonably  well-defined  intensity  threshold  and,  second, 
that  its  threshold  can  be  reached  with  strong  external  focusing 
when  F  is  constrained  to  values  less  than  Pcr.  Intrinsic  opti¬ 
cal  damage  from  avalanche  breakdown  will  be  shown  in  later 
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chapters  to  meet  these  requirements. ^  Breakdown  leads  to  a 
spark  and  local  melting  of  the  material  and  is  therefore 
easily  observed. 

Pcr  is  determined  from  measurements  of  the  damage  power  P 
made  with  various  lenses  having  different  focal  lengths.  Be¬ 
cause  P  must  be  less  than  Pcr,  short  focal  lengths  are  required. 
Complications  'rom  aberrations  in  the  focusing  optics  and  in¬ 
accuracies  in  determining  the  input  powers  can  be  greatly 
reduced  by  the  use  of  a  reference  crystal  such  as  NaCl  that  can 
be  damaged  over  a  large  range  of  input  powers  before  self- 
focusing  becomes  important.  (See  Chapts.  4  and  5*)  Approxi¬ 
mate  values  of  Pcr,  which  are  useful  for  the  initial  selection 
of  focusing  optics,  can  he  calculated  as  discussed  in  the 
previous  section. 

The  focal  lengths  of  the  various  lenses  can  be  related  to 
the  focal  soot  3ize  in  air  aQ.  (The  focal  radius  in  air  anc 
the  focal  radius  in  the  crystal  are,  neglecting  aberrations 
and  self- focusing,  identical.)  A  measurement  is  thus  made  of 
the  incident  Dower  necessary  to  induce  damage  as  a  function  of 
the  focal  radius  in  air  a0.  In  Eq.  (2-19),  I0  is  replaced  by 
kP/a0  where  k  is  a  constant  that  depends  on  the  index  of 
refraction  nQ  and  the  intensity  distribution  of  the  laser 
pulse.  If  I(r)  =  I0  exp  (-2r2/ac2),  k  has  a  value  of  (2 /n)  nQ. 
Eq.  (2-19)  then  becomes 
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The  damage  intensity  1^  does  not  depend  on  aQ  for  avalanche 
breakdown.  A  plot  of  l/P  as  a  function  of  l/a0  yields  a 
straight  line  with  l/Pcr  as  the  intercept  and  2/mI^  as  the 
slope. 

Two  precautions  should  be  taken  when  using  this  measure¬ 
ment  technique.  The  first  is  to  avoid  damage  from  absorbing 
Inclusions,  and  the  second  is  to  account  for  any  dependence 
of  avalanche  breakdown  on  the  laser  pulse  width.  Methods  for 
distinguishing  between  inclusion  and  intrinsic  damage  are 
discussed  in  detail  in  Chapt.  3.  Inclusion  damage  can  often 
be  completely  avoided,  because  the  use  of  strong  external 
focusing  which  is  needed  for  the  present  technique  results  in 
the  sampling  of  volumes  that  are  so  small  that  damaging  in¬ 
clusions  are  normally  avoided.  The  laser  pulse  width  dependence 
to  the  damage  field  will  be  shown  in  Chapt.  4  to  be  small  for 
Q-switched  laser  pulses.  Measurements  conducted  with  mode- 
locked  lasers,  however,  cannot  be  directly  compared  to  those 
made  with  Q-switched  lasers. 

As  an  illustration  of  the  use  of  this  technique,  we  con¬ 
sider  data  accumulated  in  two  different  experiments  on  optical 
damage  in  sapphire.  In  the  first  experiment-^  a  Q-switched 
TEMoo  ruby  laser  was  used  by  Giuliano  et  al.  to  damage  sapphire 
with  three  different  lenses.  In  the  second,  reported  in 
Chaot.  7,  a  measurement  was  made  with  two  different  lenses  and 

a  Q-switched  TEMqq  YAGiNd  laser.  The  pulse  widths  were  20  ns 

4? 

and  5  ns,  respectively.  Table  I  summarizes  the  relevant 
results  of  these  studies. 


TABLE  I 


Damaffe  Power  in  Sapphire 
for  Different  Beam  Sizes 


Wavelength 

a0(nm)* 

P  (MW) 

0.6943  11  m 

20 

0.51  ±  .04 

(ref.  34) 

53 

0.?4  ±  .0? 

140 

1.23  ±  .10 

1 . 064  [i  m 

9.0 

0.24  ±  .05 

(Chant.  4) 

4.1 

0.05  +  .02 

The  beam  radius  is  defined  to  be  the  l/e  radius  for  the 
intensity.  In  both  experiments  Gaussian  laser  beams  were 
used. 


In  order  to  compare  the  data,  the  wavelength  dependence 
of  must  be  considered.  It  can  be  seen  from  Eqs.  (2-10) 
and  (2-20)  that  the  data  can  be  plotted  together  if  the 
coordinate  axes  for  the  YAG  data  are  properly  scaled  relative 
to  those  for  ruby.  In  the  plot  of  l/P  versus  l/a0  in  Fig.  5i 
the  axes  are  scaled  by  £x0  (YAG)/  A0  (ruby)]2  =  2.35.  The 
intercept  on  the  vertical  ruby  axis  is  Pcr  at  0,6943  gm,  and 
the  Blooe  measured  relative  to  either  pair  of  axes,  is 
(2Ald)  n0. 

Pcr  is  found  to  be  O.83  x  10^  watts  at  ruby  wavelength, 

Pc  =  3.0  x  I0b  watts,  n2  =  O.56  x  lO"1^  esu,  and  Id  is 

1.1  x  1011  watts/cm2.  The  value  for  Pcr  lies  within  the 
range  of  values  reported  in  Ref.  >,  and  Id  agrees  with  the 
results  of  Chapt.  7. 

The  plot  iM  Fig.  5  curves  downward  when  l/w  <  0.02  x 
-2  -2 

10“  (gm)  because  the  electrostrictive  contribution  to  the 
critical  power  is  becoming  transient.  Even  though  Eq.  (2-20) 
is  not  valid  for  transient  effects,  we  can  modify  it  to 
approximate  slightly  transient  behavior  by  using  Eq .  (2-l6), 

If  it  is  assumed  that  the  electrostrictive  nonlinearity  is 
much  larger  than  the  electronic  effect,  a  damage  power  of  about 

1.1  x  106  watts  is  predicted  for  w  =  140  jt  m.  Alternately, 
if  the  electronic  and  electrostrictive  contributions  are 
equal,  the  damage  power  is  estimated  to  be  1.0  x  10^  watts. 

Both  values  compare  favorably  to  the  measured  value  of 

1.2  x  10^  watts  in  Table  I. 

q 

Giuliano  and  Marburger7  have  inferred  a  range  of  values 
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A  plot  of  1,/p  as  a  function  of  l/w^  produces  a  straight  line  according  to  F.q.  (2-20). 

In  this  figure  we  have  used  w  to  denote  the  Gaussian  beam  radius  a,,,  for  large 
wll/w^  o,  02  x  10  lp.m)^),  departures  from  a  linear  relationshi  p  a  re  expected  because 

e le ctros t r icti on  becomes  transient.  This  is  experimentally  observed  as  seen  in  the 
insert  where  the  scales  for  t  c  ruby  data  have  been  expanded. 
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for  the  nonlinear  inoex  in  sapphire  of  0.31  to  9.5  x  10"*^  esu. 
It  is  clear  from  the  more  recent  measurements  of  Ref.  3**» 
however,  that  n 2  cannot  be  as  large  as  9.5  x  10"*  ^  esu. 

Giuliano  and  Marburger's  analysis  was  based  on  the  central 
assumption  that  there  exists  a  one-to-one  correspondence  be¬ 
tween  a  damage  spark  and  a  moving  self- focus  where  the  latter 
develops  from  catastrophic  self- focusing.  While  this  is  almost 
certainly  true  at  high  input  powers,  the  present  analysis  and 
the  experiments  in  later  chapters  show  that  a  spark  and  optical 
damage  can  occur  at  low  powers  where  a  moving  self-focus  is 
absent.  It  is  very  likely  that  :heir  plot  of  versus  l/z^-, 

where  Zf  is  the  self-focusing  di  ‘  -ce,  deviates  from  a  linear 
relationshi n  at  low  powers  becau.'  damage  is  occurring  without 
catastrophic  self- focusing  in  this  region  and  not  because  of 
transient  effects  as  they  suggest.  Only  data  obtained  with 
large  values  of  P,  therefore,  can  be  used  within  the  framework 
of  their  analysis  to  predict  n2  •  The  lower  value  of  their 
range  of  values  for  n2  (O.31  x  10"* 3  esu)  is  probably  the  most 
meaningful.  It  is  somewhat  smaller  than  the  present  estimate 
of  n2 ,  because  under  their  focusing  conditions,  the  electro- 
strictive  contribution  to  the  nonlinearity  nay  have  been 
transient  and  therefore  smaller  than  under  steady-state 
conditions.  It  thus  appears  that  the  present  measurement  of 
self-focusing  parameters  in  sapphire  is  consistent  with  the 
results  of  Giuliano  and  Warburger. 

The  electronic  contribution  to  n2  can  be  determined  by 
the  presen  :  technique  if  conditions  of  measurement  are  chosen 
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such  that  tj/r^  <<i.  Because  of  the  requirement  that  the 
input  power  be  less  than  Pcr,  tp/tR  <<  1  cannot  be  achieved  by 
increasing  through  the  use  of  weakly  focusing  external 
optics.  This  extreme  transient  limit  can  be  achieved,  however, 
by  using  single,  mode-locked  laser  pulses  with  time  durations 
of  the  order  of  10  ps.  Such  pulses  must  be  temporally  smooth 
and  hence  transform-limited  so  that  a  YAGiNd  system  is  a 
natural  choice  for  the  laser. 

In  summary,  the  analysis  of  Sect,  C  leads  to  a  useful 
technique  for  measuring  self- focusing  parameters.  The  ideas 
which  we  have  developed  here  not  only  provide  the  basis  for 
this  technique,  but  they  also  help  clarify  the  relationship 
between  self-focusing  and  optical  damage,  a  relationship 
which  has  often  been  confused  in  the  literature. 


CHAPTER  3 


THE  LASER  SOURCES  AND  EXPERIMENTAL 
TECHNIQUES 


A.  Introduction 

In  this  chapter  we  will  describe  the  equipment  and  procedures  used 
in  the  dama^  ixperiments  in  Chaps.  4-7.  Sect.  B  describes  the  two 
Q-switched  laser  systems,  developed  by  D.  Bua  and  M.  Bass  of  Raytheon 
Corp.  ,  which  were  used  for  the  nanosecond-pulse  studies  of  damage. 
Measurements  made  by  the  author  to  characterize  these  systems  will  be 
presented  along  with  an  alignment  technique  that  is  useful  for  longitudinal 
mode  selection.  The  experimental  procedures  used  in  the  damage 
studies  will  also  be  summarized  in  this  section.  In  Sect.  C  a  mode -locked 
laser  developed  by  J.  P.  Leteilier  and  J.  McMahon  of  Naval  Research 
Laboratories  will  be  described  and  measurements  to  characterize  its 
performance  summarized. 

In  the  studies  of  avalanche  breakdown,  inclusion  damage  was 
treated  as  an  experimental  problem  to  be  avoided.  For  this  reason 
techniques  for  identifying  inclusion  damage  will  be  considered  in  this 
chapter  in  Sect.  D.  Without  these  techniques,  several  of  our  damage 
experiments  could  not  have  been  completed.  The  convincing  identifica¬ 
tion  of  intrinsic  processes,  in  fact,  requires  that  either  inclusion 
damage  be  known  to  be  absent  or  that  inclusion  damage  be  distinguishable 
from  intrinsic  damage  at  each  damage  site. 

B .  The  Q-Switched  Lasers  and  Beam -Handling  Optics 

Fig.  6  shows  schematically  the  principal  features  of  the 
Q-switched  laser  sources.  These  systems  were  pulse  pumped  and 

Preceding  page  blank 
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electrc-optically  Q-switched,  Some  important  properties  of 
these  devices  are  summarized  in  Table  IS. 

In  order  to  make  meaningful  measurements  of  absolute 
damage  fields,  it  was  essential  to  restrict  the  laser  trans¬ 
verse  structure  to  a  single  radial  mode.  This  was  accomplished 
by  inserting  small  circular  apertures  into  the  plane-mirror 
laser  cavities.  It  is  well  known  that  diffraction  losses  are 
large  for  a  plane  mirror  cavity  so  that  the  use  of  a  circular 
aperture,  which  reduces  the  cavity  Fresnel  number,  can  be 
effective  in  preventing  oscillation  in  all  but  the  lowest  order 
(TEWqo)  mode.^  In  a  confocal  cavity,  on  the  other  hand, 
diffraction  losses  will  be  relatively  unimportant  compared  to 
mirror  and  scattering  losses,  and  an  aperture  of  reasonable 
size  cannot  normally  produce  TEMqo  output. 

Two  techniques  were  used  to  verify  TEMqo  operation.  In 
the  first,  a  1-mil  pinhole  was  moved  across  the  YAGiNd  beam 
in  1-mil  translation  steps  and  the  intensities  of  successive 
shots  recorded.  Fig.  7  summarizes  the  result  of  one  measure¬ 
ment.  Each  data  point  is  a  visual  average  of  about  5  laser 
shots.  The  small  dip  visible  near  the  center  of  Fig.  ?a 
could  not  be  reproduced  and  therefore  indicates  the  experimental 
errors  in  the  measurement.  By  plotting  the  logarithm  uf  the 
laser  intensity  as  a  function  of  the  radisl  distance  squared, 
it  was  found  that  the  output  of  the  beam  had  a  Gaussian  profile 
of  the  form  T0  exp  (-2r  /c  )  where  a  is  measured  from  Fig.  7b 
to  be  about  0.49  mm. 

Because  the  ruby  laser  was  fired  less  frequently  than  the 
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TABLE 1 1 

LASER  PARAMETERS 


Wavelength 

Nd:YAG 

Ruby 

1.  Ofyzm 

0. 694/im 

Energy  TEMqq  Mode 

1.5  mj 

2. 0  mj 

Beam  Diameter  at  Output 

Mirror 

0. 8  mm 

0.7  mm 

TEMqq  Mode 

Polarization 

Linear 

Linear 

Pulse  Repetition  Rate 

1  pps 

i  pulse/5  sec 

Pulse  Duration  in  TEMqq  Mode 

4. 7  nsec  (FWHP) 

14  nsec  (FWHP) 

Pulse  to  Pulse  Energy 

±  7% 

+  10% 

Reproducibility 


I  1  -  -L  -  _  .  ! - 1 i  J  .  .  .  .  J 

0  200  400  600 


Square  of  radial  distance  from  beam 
center  (mils)2 

Fig.  7.  INTENSITY  DISTRIBUTION  OF  THE  YAG  :Nd  LASER  AS  A  FUNCTION 

OF  RADIAL  DISTANCE  FROM  THE  BEAM  CENTER  AT  THE  POSITION 
OF  THE  FOCUSING  LENS 
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YAG  and  was  generally  less  stable,  a  second  technique  for 
measuring  the  transverse  mode  structure  was  used  with  the  ruby 
system.  A  lensless  camera  was  placed  in  front  of  the  ruby 
laser  with  a  coated  glass  plate  positioned  just  before  a 
photographic  transparency.  The  plate  was  coated  in  such  a 
manner  that  the  red  ruby  light  was  partially  reflected  at  each 
surface  and  the  light  transmitted  to  the  undeveloped  trans¬ 
parency  consisted  of  descrete  images  each  attenuated  by  a 
factor  of  2  from  the  image  above  it.  Fig.  8  shows  the  inten¬ 
sity  distribution  of  a  single  ruby  pulse  obtained  in  this 
manner.  It  is  jeen  that  the  intensity  distribution  of  the 
ruby  laser  more  nearly  resembled  the  far-field  pattern  of  a 
plane  wave  than  a  Gaussian.  By  imaging  the  focal  plane  of 
the  lens  used  in  the  damage  studies,  it  was  possible  to 
directly  measure  the  intensity  distribution  of  the  laser  beam 
when  it  was  tightly  focused.  This  technique  was  not  used  with 
the  YAG  laser  because  of  the  difficulties  involved  in  obtaining 
and  preparing  photographic  emulsions  sensitive  to  1.06gm 
radiation. 

Measurements  of  the  transverse  distribution  summarized  in 

Fies.  ?  and  8  integrate  the  transverse  intensity  distribution 

over  the  full  time  duration  of  the  laser  pulse.  Depletion 

effects  in  the  rod,  however,  can  cause  this  distribution  to 
,  .  44 

change  with  time.  To  verify  that  such  an  effect  was  not 
occurring  on-axis  where  the  laser  intensity  was  greatest  and 
damage  was  known  to  commence,  the  center  of  the  beam  was 
sampled  with  a  25-micron  pinhole  and  found  to  have  the  same 
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time  structure  as  the  entire  beam. 

The  lasing  band-width  of  the  lasers  was  restricted  to  at 
most  a  few  adjacent  longitudinal  modes  by  using  a  particularly 
simple  technique  that  does  not  require  any  special  optical 
cavity  componants  such  as  an  output  resonant  reflector.  The 
technique  consists  of  aligning  the  front  surface  of  a  plane- 
cut  laser  rod  parallel  to  the  100  percent  reflecting  mirror. 

Fig.  9  indicates  the  resulting  intracavity  resonator  Mj-A. 

The  anti-reflection  coating  on  each  end  of  the  laser  rod  has 
a  residual  reflectivity  of  about  0.25  percent  and  so  the 
passive  finesse  of  resonator  Mj-A  is  very  low.  When  lasing 
occurs,  however,  F^-A  contains  a  medium  with  gain  and  its 
effective  finesse  can  be  very  high.  In  our  ruby  laser  with  a 
10  cm  rod,  for  example,  the  combination  of  a  0.25  percent 
reflection  surface.  A,  a  laser  gain  coefficient  of  =  0.25  cm  ”  , 
and  the  99.9  percent  cavity  reflector,  Mj  ,  results  in  a  steady- 
state  reflectivity  finesse  of  about  6.  All  other  intracavity 
elements  are  placed  between  surfaces  A  and  M2  to  minimize 
reflection  and  scattering  losses  In  the  subcavity  M^-A, 

In  figures  10  and  11  the  effects  of  rod  alignment  on  the 
temporal  outputs  of  the  NdiYAO  and  ruby  laser  are  illustrated 
as  measured  by  a  photodiode-oscilloscope  combination  having 
a  measured  risetime  of  0,5  nsec.  For  the  ruby  laser,  additional 
mode  selection  was  achieved  by  using  a  single  sapphire  etalon 
as  the  output  mirror,  M2. 

The  detection  of  a  smooth  pulse  by  a  limited  bandwidth 
detector  maj  be  the  result  of  either  single-frequency  output 
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or  the  simultaneous  oscillation  of  many  randomly  phased  modes. 
In  order  to  resolve  this  ambiguity  the  lasers*  output  spectra 
were  studied  with  Fabry-Perot  interferometers.  The  lasing 
bandwidth  narrowed  significantly  as  the  rod  was  aligned  to 
achieve  smoother  pulses.  Unfortunately,  our  interferometers 
lacked  sufficient  finesse  to  resolve  individual  laser  cavity 
modes  (A  v  ~  300  MHz).  The  longitudinal  mode  content  of  a 
perfectly  smooth  laser  pulse  can  be  inferred,  however,  by 
combining  the  interferometer  and  pulse  waveform  data.  The 
observed  interferometer  limited  bandwidth  of  a  smooth  ruby 
laser  pulse  was  =1.2  GHz,  implying  that  a  maximum  of  4 
adjacent  modes  could  have  been  oscillating.  Since  the  photo- 
diode-oscilloscope  combination  could  detect  frequencies  as 
high  as  -  3  Ghz,  the  presence  of  four  adjacent  oscillating 
modes  would  have  been  detected  through  the  observation  of  mode 
beating  as  in  Fig.  lib.  Therefore,  our  smooth  ruby  laser 
pulses,  as  shown  in  Fig.  11c,  were  caused  by  single  longi¬ 
tudinal  mode  lasing.  In  support  of  this  conclusion  we  note 

i4,C 

that  Giuliano  et  al.  resorted  that  perfectly  smooth  pulses 
correspond  to  single  mode  oscillation  while  non-smooth  wave¬ 
forms  were  always  associated  with  two  or  more  oscillating 
modes.  Replacing  the  sapphire  etalon  with  a  wedged  output 
mirror  did  not  qualitatively  change  the  results. 

We  have  found  that  slight  misalignment  of  a  multi¬ 
component  resonant  reflector  or  of  one  which  employs  high 
reflectivity  mirrors  often  causes  transverse  mode  distortion. 
The  simple  technique  for  obtaining  smooth  pulses  described 
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above  does  not  have  this  problem  and  is  quite  compatible  with 
TEMqq  node  operation. 

When  the  product  of  the  double-pass  gain  and  the  output 
reflector  is  less  than  about  40  percent,  it  is  no  longer 
possible  to  operate  the  laser  in  a  single  longitudinal  mode. 

A  simple  modification  of  the  present  technique  is  to  misalign 
the  rod  and  to  insert  a  coated  optical  flat  parallel  to  the 
cavity  mirrors.  One  surface  of  the  flat  functions  as  the  rod 
face  Mj  above,  and  its  reflectivity  can  be  chosen  to  be  high 
enough  to  give  high  mode  discrimination  while  still  low 
enough  so  that  there  is  no  lasing  from  that  surface.  This 
i:onf!  juration  is  identical  in  basic  form  to  conventional 
resonant  reflectors  except  that  the  first  surface  has  a  low 
reflectivity  to  minimize  degradation  of  the  transverse 
structure.  As  in  the  previous  technique,  mode  discrimination 
is  achieved  by  providing  a  subcavity  containing  an  active 
medium. 

We  thus  see  that  the  Q-switched  lasers  used  in  the  present 
study  produce  TEP'oo  output  with  at  most  a  few  simultaneously 
lasing  longitudinal  modes.  Because  these  systems  are  well- 
characterized  and  stable,  it  is  possible  to  use  them  to  make 
accurate  determinations  of  damage  fields. 

The  actual  damage  measurements  were  made  by  focusing  the 
laser  beams  through  a  1 4  mm  focal  length  lens  to  approximately 
2  mm  inside  the  samples.  C'  e  was  taken  to  insure  that 
spherical  aberrations  from  both  the  lens  end  the  plane  en¬ 
trance  surface  of  the  sample  being  tested  we~e  unimportant , 
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A  brief  discussion  of  aberrations  is  given  in  the  Appendix. 

A  fast  photodiode  was  used  to  monitor  the  transmitted  light, 
and  an  energy  monitor  recorded  the  energy  in  each  laser  pulse. 
The  combination  of  one  rotatable  and  one  fixed  polarizer 
resulted  in  a  variable  light  attenuator  which  was  highly 
sensitive,  quite  renroducible ,  and  which  did  not  affect  the 

1 2 

laser  pulse’s  polarization,  spatial  distribution,  or  duration. 

If  the  fixed  polarizer  is  oriented  to  transmit  the  laser 
polarization  and  if  d  =  0°  is  the  angle  of  the  rotating 
polarizer  which  gives  maximum  transmission  through  this 
attenuator,  then  the  transmitted  intensity  at  any  other  angle 
of  rotation  about  the  beam  axis  is 

1(8)  -  b  I  cof.”1  8 
o 

where  i'0  is  the  incident  light  intensity  and  b  is  the  fraction 
transmitted  when  6  =  0°.  Calibrated  neutral  density  filters 
were  often  used  in  conjunction  with  the  variable  Gian 
attenuator. 

C.  The  Mode-Locked  Laser 

When  a  large  number  of  randomly  phased  longitudinal 
modes  oscillate  simultaneously,  the  laser  output  has  an 
irregular  spirted  envelope.  A  phase  correlation  CL_n  be  forced 
on  adjacent  modes,  however,  and  a  train  of  short  pulses  pro¬ 
duced.  The  technique  most  often  used  for  this  phase-  or  mode- 
locking  is  to  place  a  saturable  organic  dye  solution  inside 
the  cavity. ^  Because  the  absorption  of  the  dye  decreases 


r 


I 

fc 


for  hl?h  intensities,  the  field  in  the  cavity  is  forced  into 
a  high-intensity,  short-duration  wave-packet,  or  equivalently, 
a  large  number  of  phase-correlated  side-bands  are  generated 
which  coincide  with  the  longitudinal  mode  resonances  of  the 
cavl  ty  and  which  together  produce  the  desired  output. 

The  mode-locked  laser  used  for  the  subnanosecond  damage 
study  of  Chapt.  5  was  a  YAfJiNd,  passively  mode-locked  laser 
with  kodak  dye  97^0  used  as  the  bleachable  dye.  This  system, 
as  opposed  to  a  glasstNd  laser,  has  the  important  feature  that 
the  individual  spikes  in  the  mode-locked  train  are  free  of 
time  and  carrier-frequency  substructure.®®  It  is  possible, 
in  other  words,  to  establish  a  Dhase  correlation  over  the 
entire  lasing  band-width  with  YAC:Nd  as  the  active  medium. 

In  order  to  allow  a  large  number  of  longitudinal  modes 
to  reach  oscillation  simultaneously,  the  ends  of  the  rod  were 
cut  at  Brewster  angle.  A  circular  aperture  restricted  thp 
radial  distribution  to  a  T^Vqq  outnut,  and  the  tempera  lur  e 
stabilized  dve  was  carefully  circr'  i’  .  •  •;  ough  a  1  mm  thick 

cell  containing  the  rear  99.9%  refxt  '  r  Die  system  was 
operated  at  a  firing  rate  of  1  pps  with  an  amplitude  stability 
of  about  2  03  in  intensity. 

A  laser- triggered ,  electro-optic  shutter,  positioned  just 
after  the  laser  outnut  mirror,  was  used  to  select  a  single 
pulse  from  the  mode-locked  train.  The  operation  of  this 
device  is  summarized  in  'j  2,  Alter  the  nulse  train 

travels  through  an  initially  unbiased  KDF  Pcckols  cell,  it  is 


reflected  from  a  fixed  Plan  polar*  zer  arid  focused  onto  the 
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high-voltage  electrode  of  a  pressurized  laser-triggered  spark 
gap.  When  the  gas  between  the  electrodes  breaks  down,  a  5  KV 
voltage  spike  is  transmitted  down  a  coaxial  line  to  the 
Pockels  cell.  This  voltage  spike  has  a  time  duration  of  about 
7  ns  which  coincides  with  the  round-trip  cavity  time  and 
therefore  the  temporal  separation  of  the  pulses  in  the  train. 
During  this  7  ns  time  interval,  one  pulse  will  enter  the 
Pcckels  cell.  Because  of  the  high  voltage  applied  longitu¬ 
dinally  across  the  KDP  crystal,  the  plane  of  polarization  of 
that  pulse  is  rotated  approximately  90°  and  allowed  to  pass 
through  the  Olan  polarizer. 

With  this  desien  we  were  able  to  obtain  single  pulses  of 
approximately  15  ps  duration  for  the  damage  studies.  By 
replacing  the  cavity  output  mirror  with  a  sapphire  etalon,  the 
lasing  bandwidth  was  restricted,  and  pulses  of  just  under 
300  ds  duration  were  obtained.  Two-photon  absorption  measure¬ 
ments  established  the  pulse  durations. 

The  energy  calibration  posed  a  problem  because  a  single 
15  ns  pulse  caused  optical  breakdown  Inside  a  i  bra ted 
ballistlcal  thermal-pile  of  the  type  that  had  beer,  used  to 
measure  the  energy  of  the  Q-switched  lasers.  A  planar  r.-noio- 
voltaic  cell  was  finally  chosen  as  an  energy  nx.'rior.  :;v 
operating  it  in  an  integrating  mode  and  cal  r.  *  v;  \  r.r  it  wi  r.h 
a  commercial  Q-switched  YAO  laser,  we  were  able  to  esrar]  ,  s  n 
the  energy  calibration  to  within  about  2  5*?. 

This  system  had  been  known  to  operate  in  a  mode 

with  an  intensity  distribution  beet  describe-  r.s  r  :  •.;? 
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Oausslan.  (The  distortion  of  the  mode  probably  resulted  from 
the  optical  characteristics  of  the  dye.)  We  verified  that 
the  mode  pattern  had  not  changed  by  recording  successive  burn 
patterns  on  a  commercial  blackened  paper  as  the  light  attenu¬ 
ation  was  varied  in  fixed  steps. 

The  actual  beam  handling  optics  beyond  the  electro-optic 
shutter  were  the  same  as  with  the  Q-switehed  lasers  except 
that  laser  attenuation  was  accomplished  by  calibrated  filters 
alone. 


D.  The  Distinction  Between  Intrinsic  and  Inclusion  Induced 
Damage 


Transparent  materials  often  contain  absorbing  inclusions 
which  are  too  small  to  be  easily  detected  by  optical  in¬ 
spection.  Calculations  in  Chapt.  1  have  shown  that  when  tneir 
diameters  are  greater  than  about  0.1  gn,  such  inclusions  can 
produce  local  melt!  ne-  o*'  the  surrounding  material  when  ir¬ 
radiated  with  Q-swi  tch«*d  laser  ruls*‘r. 
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utilizes  the  damage  statistics  described  in  Chapt.  6  and  can 
be  applied  to  both  surface  and  bulk  studies.  Its  use,  however, 
requires  a  large  number  of  l:.r'nt  probes  by  a  highly  stabilized 
laser.  A  third  technic  up ,  which  was  developed  during  the 
present  study,  allows  an  immediate  determination  of  the 
general  source  of  opt:,  cal  damage  at  individual  sites  both  in 
th*»  bulk  and  on  tr.e  s;u-  fuco.  This  technique  does  not  require 
an  amolitud-5  stabi  1 :  ted  laser  and  thus  simplifies  the  problem 
of  ’  r .  v  i  ‘v  iv:  i  r trlnsi  c  surface  and  bulk  damage.  Finally,  a 
nuv.’oer  of  qualitative  observations  to  be  described  below  are 
part.’.  :-j  1  ar'i y  useful  as  indicators  of  inclusion  damage  when 
sutnano>?-»ro::i  i.»u  loss  are  u:.n-i  and  the  first  and  third  tests 
cannot  be  a  nulled. 
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in  the  bulk  by  examining  the  damage  morphology.  Figs.  13  and 
14  show  examples  of  this  morphology  as  seen  in  the  damage 
study  described  in  Chapt.  k. 

The  second  technique  for  determining  the  absence  of 
inclusion  damage  is  based  on  the  statistical  nature  of  intrinsic 
damage  as  first  observed  by  Bass  and  Barrett  '  and  extended  in 
Chapt,  6.  A  precisely  defined  threshold  for  intrinsic  damage 
does  not  exist.  As  a  result,  intensity  levels  can  be  found 
such  that  damage  occurs  on  each  pulse  with  some  finite 
probability  pj  defined  by 


p  _  total  number  of  damages 
H1  total  number  of  pulses 


(3-D 


Bv  fixing  the  intensitv  of  the  highly-stabilized  laser  source, 
z  distribution  function  fN  can  be  obtained  which  describes  the 
fractional  number  of  times  N  puises  are  required  to  produce 
damage.  If  laser  fluctuations  can  be  neglected  arid  p-^  is  the 
same  for  each  laser  oulse  ( p^  =  0  or  1  for  inclusion  damage), 
f jvj  is  given  by  the  probability  that  no  damage  has  occurred  on 

the  first  N-l  pulses  times  the  probability  it  has  occurred  on 
jth 


the  N' 


i  t  0  • 


r  -  /■  ,N  -  i 

1  “  P  |  )  Pi 


(3-2) 


Since  the  only  way  this  distribution  of  N  values  can  be 
obtained  is  if  is  constant,  the  agreement  between  measured 
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Geometric  Focal 
Plane 


ig.  14.  INCLUSION  DAMAGE  IN  RbCl 


Damage  which  we  regarded  as  due  to  absorbing  inclusions  consisted 
of  small  spheres  randomly  distributed  about  the  geometric  focal 
plane.  The  photograph  was  taken  during  the  experiment  oi  Chapt.  4. 
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values  of  f]\j  and  those  computed  by  Fq.  (3-2)  indicates  that 
material  inhomogeneities  were  not  measurably  affecting  the 
damage  data  and  that  the  laser  was  sufficiently  stable  to 
observe  damage  statistics.  An  example  of  the  use  of  this 

s'?- 

technique  from  the  measurements  of  Chapt.  6  is  given  in  Fig.  15 
where  the  measured  distribution  for  bulk  damage  in  fused 
quartz  and  Fq .  (3-2)  are  seen  to  agree. 

A  third  method  for  distinguishing  between  inclusions  and 
intrinsic  damage  was  develoned  and  used  in  the  oresent  work. 

It  consists  of  observing  the  light  transmitted  through  the 
sample  by  using  a  fast  photodiode-oscilloscope  combination. 

When  damage  occurs,  the  transmitted  laser  pulse  is  attenuated 
in  a  manner  which  is  found  to  be  characteristic  of  the  source 
of  damage. 

Pulses  which  are  attenuated  very  raoidly  form  damage 
regions  in  the  bulk  which  are  characteristic  of  the  intrinsic 
mechanism.  In  addition,  for  materials  which  are  virtually 
threshold-like  such  as  those  studied  in  the  present  work,  the 
intensity  at  the  instant  of  attenuation  for  such  pulses  varies 
by  no  more  than  about  25  nercent  abjut  some  average  value. 

An  example  of  this  type  of  attenuated  pulse-shape  is  given 
in  Fig,  l6b. 

It  was  found  in  the  experiments  of  Chapts.  4-7,  however, 
that  a  number  of  damaging  pulses  were  attenuated  in  a  very 
different  manner  and  produced  damage  regions  that  could  be 
identified  as  resulting  from  inclusion  absorption.  For  this 
class  of  damage  events,  the  attenuation  was  not  as  rapid  and 
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'ig.  15.  COMPARISON  OF  MEASURED  f  WITH  p^l-p^ 
INTRINSIC  BULK  DAMAGE  IN  FUSED  QUARTZ 


BARS  ARE  MEASURED 
DISTRIBUTION  OF  N  VALUES 
.5-  FOR  48  MEASUREMENTS  OF 
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is  the  fractional  number  ol  Limes  N  pulsus  were  required 
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as  coinplete  as  with  intrinsic  damage  nor  were  the  transmitted 
pulse  shape  and  intensity  at  the  first  instant  of  attenuation 
as  repeatible.  Examples  of  this  type  of  attenuation  are  given 
in  Figs.  16c  and  d. 

There  is,  in  addition,  a  small  percentage  of  damage 
events  which  could  not  be  unambiguously  identified  by  this 
scheme  or  by  microscopic  investigation  of  the  damage  site. 

In  the  laser  damage  studies  data  obtained  from  such  pulses 
were  not  used. 

Another  distinction  between  inclusion  and  intrinsic 
breakdown  involves  the  intensities  of  the  sparks  created  by 
damaging  laser  pulses.  In  general,  intrinsic  events  produced 
bright  sparks  while  damage  from  inclusions  occasionally 
resulted  in  sparks  which  were  barely  visible. 

Except  for  extremely  high  laser  intensities,  intrinsic 
breakdown  occurs  near  the  peak  of  the  laser  pulse  so  that  the 
energy  which  is  available  to  take  part  in  the  damage  process 
is  consistently  about  half  the  total  pulse  energy.  This  was 
observed  experimentally  and  can  be  explained  by  theoretical 
considerations  of  avalanche  breakdown  (see  Chapt.  8).  It  is 
thus  to  be  expected  that  the  energy  deposited  in  the  focal 
region  is  reasonably  constant  for  intrinsic  breakdown  and  that 
the  intensity  of  the  sparks  which  apparently  always  accompany 
such  damage  is,  therefore,  also  reasonably  constant. 

Figs.  l6c  and  d  suggest  that  these  same  conclusions  do 
not  hold  for  inclusion  damage,  and,  indeed,  there  is  no  reason 
why  they  should.  Damage  from  inclusion  absorption  is  energy- 
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dependent  for  many  different  inclusion  sizes,  and  actual 
material  damage  may  not  develop  until  very  late  in  the  laser 
pulse  ( i ,e.  ,  until  most  of  the  optical  energy  has  passed 
through  the  sample).  The  focal  region  may  thus  be  melted  but 
not  heated  sufficiently  to  produce  a  bright  spark.  A  number 
of  experimenters  have,  in  fact,  reported  material  damage 
without  sparks,^*  ^9  <p^e  ar^uments  advanced  here  indicate 
that  some  type  of  extrinsic  absorptive  mechanism  such  as 
inclusion  absorption  very  likely  caused  the  damage  in  their 
work. 

Two  of  the  techniques  discussed  here — examination  of  the 
damage  morphology  and  examination  of  the  transmitted  laser 
light--could  not  be  used  in  the  subnanosecond  study  of  in¬ 
trinsic  damage  in  NaCl  (Chant.  *5)  because  of  the  small 
volumes  of  the  damaged  sites  and  the  short  durations  of  the 
laser  pulses.  But  oecause  breakdown  is  virtually  threshold¬ 
like  in  NaCl,  ?.  number  of  qualitative  observations  were  used 
to  indicate  the  absence  of  inclusion  damage.  For  the  sample 
?d  to  obtain  the  actual  damage  data,  it  was  found  that  the 
damage  field  was  reasonably  well-defined  and  did  not  change  as 
different  regions  of  the  sample  were  probed  and  lenses  with 
different  focal  lengths  were  used  to  focus  the  radiation.  Also, 
only  one  faint  spark  occurred  with  each  damaging  laser  pulse, 
and  the  spark  always  appeared  to  form  at  the  geometrical  focal 
plane.  These  observations  contrasted  with  those  obtained 
under  conditions  where  inclusion  damage  had  been  seen  in 
another  NaCl  crystal  irradiated  by  the  same  laser  and  in  other 
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samples  irradiated  by  Q-switched  lasers.  It  was  therefore 
concluded  that  except  possibly  for  occasional  damage  sites, 
inclusion  damage  was  absent  in  highly  pure  NaCl  under  the 
conditions  of  measurment  described  in  Chapt.  5. 

In  summary,  then,  several  techniques  have  been  used  to 
either  confirm  the  absence  of  inclusion  damage  or  to  dis¬ 
tinguish  between  inclusion  and  Intrinsic  breakdown  when  in¬ 
clusion  damage  occurs.  Intrinsic  breakdown  cannot  be 
effectively  studied  without  such  techniques. 


CHAPTER  4 


INTRINSIC  OPTIMAL  DAMAGE  IN  THE  ALKALI 
HALIDES  INDUCED  BY  1.06  pm  RADIATION 


A.  Introduction 


Electron  avalanche  breakdown  had  been  mentioned  by  a  number  o£ 

authors  ^  as  a  possible  mechanism  for  damage.  Until  the  work  of 

14 

Yablonovitch,  *  however,  no  simple  technique  had  been  suggested  for 

identifying  the  occurrence  of  avalanche  breakdown  in  solids.  The  situation 

was  much  different  in  gas  studies.  ^  Not  only  had  the  details  of  dc 

5  3 

avalanches  been  better  understood  in  gases  than  in  solids,  but  the  ability 
to  vary  pressure  allowed  investigators  of  gas  breakdown  to  convincingly 
identify  the  development  of  avalanche  breakdown  ar.d  to  ascertain  the 
conditions  under  which  multiphoton  ionization  could  develop  before  an 
electron  avalanche. 

Yablonovitch 's  basic  approach  was  to  select  transparent  materials 

that  had  been  well-studied  for  their  dc  breakdown  characteristics  and 

for  which  damage  from  multiphoton  absorption  could  not  occur.  The 

alkali  halide  family  was  a  natural  choice  for  laser  breakdown  studies. 

Because  they  have  large  band-gaps,  these  compounds  are  transparent 

from  about  15  pm  to  about  0,  2  pm  and  optical  damage  in  them  from 

multiphoton  absorption  is  highly  improbable  until  frequencies  have  been 

reached  that  arc  well  into  the  visible  spectrum.  The  alkali  halides  have, 

54 

in  addition,  been  studied  extensively  in  dc  damage  experiments. 
Yablonovitch  found  that  the  relative  root-me an- square  damage  fields  for 
nine  alkali  halides  at  10.  6  pm  was  the  same  as  had  been  measured  at  dc. 
On  the  basis  of  this  experimental  measurement  of  breakdown  strength 
plus  the  prediction  that  multiphoton  absorption  could  not  cause 
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damage  at  10.6  pm,  it  was  concluded  that  the  intrinsic  damage 
process  was  an  electron  avalanche,  similar  in  basic  character 
to  dc  dielectric  breakdown. 

The  discovery  that  the  avalanche  breakdown  caused  by  a 
laser  pulse  at  10.6  pm  in  the  alkali  halides  appeared  to 
proceed  by  the  same  mechanism  and  to  have  the  same  threshold 
as  dc  electric  breakdown  can  be  explained  if  simple  impact 
ionization  theories^0*  55  with  electron  collision  times 
shorter  than  the  light  period  are  correct.  As  discussed  in 
Chapt.  1,  the  electron  avalanche  develops  because  electrons 
are  heated  by  an  electric  field  in  the  presence  of  phonon 
collisions  that  primarily  change  the  electron  momentum.  Al¬ 
tering  the  frequency  of  the  applied  field  should  have  little 
effect  on  this  heating  until  the  fiela  reverses  direction  in 
a  time  comparable  to  the  electron- phonon  collision  time.  If 
only  the  rate  of  conduction-band  electron  heating  determines 
the  avalanche  characteristics,  the  avalanche  should  therefore 
appear  to  be  in  its  dc  limit  as  Ions:  as  cv<r-^  where  o»  is  the 
laser  angular  frequency  and  r  is  the  average  electron-phonon 
collision  time.  (This  behavior  is  displayed  by  Eq .  (1-1)  of 
Chapt.  1.)  An  upper  limit  r  x  to  the  collision  time  in  NaCl 
can  be  estimated  by  using  the  low- field  mobility^  (10  cm^/volt- 
sec)  and  setting  it  equal  to  ermax/me  where  me  is  the  free- 
electron  mass  and  e  is  the  electronic  charge.  We  find  that 
a;rmax  =  i  f°r  10  Mm  radiation.  In  electric  fields  approaching 
breakdown  values,  the  average  electron  drift  energy  can  be 
3hown  to  be  comparable  to  the  longitudinal  optical  (LO)  phonon 


energy,  and  since  the  scattering  is  enhanced  near  the  LO 
c? 

resonance,  r  should  be  less  than  rmax .  Thus  Yablonovitch 
was  working  at  a  frequency  where  cur  was  nrobably  still  a 
small  quantity  compared  to  unity  and  where  simple  impact 
ionization  theories  would  predict  essentially  dc  behavior. 

We  have  extended  Yablonovitch*  s  work  to  1 . 06  pm  by  using 
a  Q-switched  YAf.tNd  laser  to  induce  damage.  The  intent  of 
this  work  was,  first,  to  determine  if  the  relative  intrinsic 
damage  fields  of  the  alkali  halides  are  different  at  1.06  p m 
than  they  are  at  10.6  pm  and  at  dc  and,  second,  to  demonstrate 
the  feasibility  of  conducting  near-IR  and  visible  bulk  optical 
damage  studies  at  low  Dowers  without  catastrophic  self- 
focusing.  It  will  be  shown  below  that  the  relative  breakdown 
strengths  in  fact  change  very  little  at  1.06  pm  and  that  a 
number  of  additional  time-related  observations  support  the 
conclusion  that  intrinsic  breakdown  is  occurring  by  a  time- 
dependent  electron  avalanche. 

Self-focusing  is  ootentially  more  of  a  problem  Jt  1 . 06  p  m 
than  it  is  at  10.6  pm.  As  discussed  in  Chapt,  2,  the  amount 
of  beam  distortion  from  the  self- focusing  nonlinearity 
depends  on  P/Pcr.  Since  it  is  desirable  to  have  as  little 
beam  distortion  as  possible,  P/Pcr  should  be  kept  small.  As 
the  wavelength  increases,  however,  Pcr  decreases  as  A2,  so 
that  to  maintain  the  same  value  of  P/P  ,  P  must  also  scale 

C  i 

2 

as  A  .  Because  the  breakdown  intensity  is  approximately 
independent  of  A,  the  requirement  that  P/Pcr  be  held  constant 
is  equivalent  to  requiring  that  the  area  of  the  focal  spot-- 
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2  2 

-(X/d)  f  where  d  is  the  incident  beam  diameter  and  f  is 

o 

the  lens'  focal  length--also  scale  as  A'.  But  the  diameter  of 
a  TEMqq  mode  depends  on  both  wavelength  and  cavity  design, 
being  on  the  order  of  1  cm  for  the  10.6  pm  laser  of  Ref.  14 
and  0.1  cm  for  our  1.06  pm  system.  In  order  to  maintain  a 
constant  value  of  P/Pcr  as  X  is  decreased,  it  is  necessary  to 
focus  more  tightly  with  an  external  lens  or  to  design  the 
cavity  optics  to  compensate  for  the  dependence  of  d  on  wave¬ 
length. 

In  the  present  work  it  was  not  possible  to  focus  suf¬ 
ficiently  to  maintain  the  same  value  of  P/Pcr  that  Yablonovitch 
had  been  able  to  achieve  at  10.6  pm.  Pcr  changes  by  a  factor 
of  100  as  A  is  decreased  from  10.6  pm  to  1.06  pm,  and  the 
minimum  value  of  f  that  could  be  used  was  limited  by  aber¬ 
rations  to  about  \  inch  as  compared  to  f  =  1  inch  for  the  lens 
that  Yablonovitch  had  used.  As  a  result,  self-focusing  was 
potentially  more  of  a  problem  in  our  measurements,  and  it  was 
necessary  to  conduct  simple  tests  to  confirm  the  absence  of 
significant  self- focusing  effects.  The  designs  of  these  tests 
were  discussed  in  Chapt.  2,  and  their  results  are  summarized 
in  the  next  section. 

B,  The  Experiment 

The  laser  sources  and  experimental  arrangement  for  the 
present  work  are  summarized  in  Section  B  of  Chapt.  3.  In  this 
section  we  will  summarize  the  actual  experiment,  including 
both  the  self-focusing  tests  and  the  measured  breakdown  fields 
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of  nine  alkali  halides  at  1.06  gm  and  one  at  0,69  gm. 

As  discussed  in  detail  in  Chapt.  2,  self- focusing;  effects 
can  lead  to  significant  beam  distortion.  Because  accurate 
measurements  of  breakdown  fields  are  desired,  it  is  necessary 
to  take  great  care  to  avoid  such  distortion  by  restricting 
the  laser  input  power.  Theoretical  self- focusing  parameters 
were  defined  and  derived  in  Chapt.  2  where  quantitative 
corrections  from  the  index  nonlinearity  at  powers  below  PQ 
were  discussed.  Table  III  summarizes  the  numerical  results. 

The  input  power  is  the  experimental  power  at  the  peak  of  the 
laser  pulse  and  is  more  than  one  order  of  magnitude  below  Pc. 
From  a  purely  theoretical  viewpoint,  therefore,  catastrophic 
self- focusing  is  impossible,  and  it  can  be  shown  using  Eq. 
(2-is)  that  beam  distortion  from  the  index  nonlinearity  intro¬ 
duces  at  most  a  few  percent  correction  in  the  measured  electric 
field  strengths.  If  catastrophic  self-focusing  does  occur, 
then  the  breakdown  damage  data  is  a  measure  of  the  critical 
powers  rather  than  the  intrinsic  breakdown  field.  The  measured 
threshold  intensity  will  then  scale  with  the  square  of  the 
calculated  focal  diameter  if  the  process  is  steady-state  and 
will  depend  on  the  pulse  width  if  the  process  is  transient. 

(The  diameter  dependence  in  the  steady-state  results  from  the 
existence  of  a  constant  critical  power  Pc  which  does  not  vary 
with  beam  diameter.) 

Two  experiments  were  conducted  to  test  our  belief  that 
self- focusing  was  absent.  In  the  first  the  relative  field 
strength  threshold  for  damage  in  NaCl  was  measured  with  three 
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TABLE  III 


Calculated  Steady-State,  Self-Focusing  Parameters  and 
Experimental  Values  of  Pulse- Widths  and  Peak  Power 


Wavelength 

t 

P 

v  ,  ,22 

n^X  10 

P 

cr 

r 

c 

P.  . 

input 

(microns) 

(10  ^sec) 

(10  Jsec) 

(mks) 

(lO^watts) 

( 1  Oswalts) 

(lO^watts! 

10.6 

5.5 

200 

48, 000 

175, 000 

120 

NaCl 

1.06 

2.7 

4.7 

2.3 

480 

1,750 

30 

0.69 

2.0 

14 

2  04 

746 

26 

10.  6 

11.2 

200 

13,200 

50, 000 

20 

Rbl 

1.06 

5.4 

4.7 

8.  1 

132 

500 

8.  1 

For  Pjn  <  Pc  catastrophic  self-focusing  will  not  occur. 

The  10.6  prr>  data  is  taken  from  reference  14. 

See  Chapt.  2  for  definitions  of  n„,  and  P  .  To  convert  n„  to  esu  units, 

R  2  c  r  2 

multiply  by  0.91  X  10^. 
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different  focusing  lenses,  corrected  for  spherical  aberrations, 
and  having  focal  lengths  of  1.4,  2.5  and  3.0  cm.  The  experi¬ 
ment  was  conducted  at  1,06  If  steady-state  self- focusing 

were  present,  the  observed  damage  threshold  would  have  scaled 
with  the  inverse  of  the  focal  length.  It  did  not,  and,  in 
fact,  to  within  5  percent  the  field  strength  was  independent 
of  focal  length.  This  effectively  eliminated  the  possibility 
of  steady- state  self-focusing.  Since  tp/tR  >  1  from  Table  III, 
self- focusing  should  not  be  transient.  Eq.  (2-17)  however, 
predicts  the  results  observed  when  transient  self- focusing  is 
present.  For  this  reason  a  measurement  was  made  of  the  damage 
threshold  as  a  function  of  pulse  duration  with  the  beam 
diameter  held  essentially  constant. 

By  changing  the  numbing  level  for  the  YAG  laser,  we  were 
able  to  extend  the  pulse  width  by  a  factor  of  2.3  to  10.8  nsec. 
In  addition,  the  breakdown  strength  at  0.69  pm  was  measured 
with  ruby  laser  pulses  of  14-nanosecor.d  duration  and  a  focused 
diameter  25  percent  smaller  than  that  obtained  with  the  YAG 
laser.  The  same  1.4  cm  focal  length  lens  was  used  in  all  three 
measurements,  and  to  compute  the  ruby  value,  we  assumed  the 
same  transverse  intensity  variation  as  that  present  at  1.06  pm. 
To  within  15  percent  no  change  was  noted  in  the  threshold 
field  despite  the  oulse-width  dependence  in  Eq.  (2-17).  The 
agreement  for  the  ruby  pulses  was  especially  reassuring,  be¬ 
cause  the  critical  power  varies  with  wavelength  squared.  If 
transient  self-focusing  were  present,  we  would  have  seen  a 
change  by  a  factor  of  18  in  the  measured  intensity--an  effect 
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which  would  have  been  quite  dramatic.  A  factor  of  nine 
comes  from  the  pulse-width  dependence  of  the  transient  critical 
power  and  a  factor  of  2  from  the  wavelength  dependence. 

Perhaps  the  best  experimental  check  for  self-focusing:  is 
the  actual  measurement  of  breakdown  strengths.  Self-focusing 
theory  appears  to  be  totally  unable  to  account  for  the  experi¬ 
mental  results  given  below  in  which  both  relative  and  absolute 
values  of  breakdown  strengths  show  striking  similarities  to 
10,6  pm  values.  We  thus  conclude  that  prior  to  the  onset  of 
material  damage,  self-focusing  has  been  effectively  eliminated 
as  a  competing  nonlinearity. 

The  possibility  may  exist  that  self- focusing  occurs  after 
a  sufficient  number  of  electrons  have  been  generated  to  cause 
intense  local  heating  of  the  sample.  We  note,  however,  that 
in  our  measurements  any  late  developing  nonlinearity  is  unim¬ 
portant. 

To  measure  the  breakdown  strengths  of  the  alkali  halides, 
the  laser  beam  was  focused  approximately  2  mm  into  each  sample 
and  the  number  of  laser  pulses  necessary  to  produce  internal 
damage  at  various  power  levels  was  recorded.  In  every  case 
where  damage  occurred,  a  white  spark  was  produced,  and  the 
damage  was  later  carefully  inspected  with  a  microscope.  Be¬ 
cause  of  the  small  volume  damaged  by  our  highly  focused  1.06gm 
pulses  (less  than  2  x  10"  cm;) ,  a  large  number  of  data 
points  (40  to  100)  could  be  taken  with  each  sample. 

Defining  threshold  as  that  value  of  incident  power 
necessary  to  produce  intrinsic  damage  in  a  single  ahot  tor 


50  percent  of  the  Dositions  probed,  we  calculated  the  rms , 
on-axls  electric  field  at  the  measured  threshold  in  NaCl. 
Corrections  were  made  for  reflections  from  various  surfaces 
and  the  changes  in  the  beam  diameter  due  to  the  effect  of  the 
index  nonlinearity.  This  was  the  basic  calibration,  and  all 
other  values  of  threshold  were  measured  relative  to  ENaC1. 

In  order  to  avoid  errors  from  daily  power  fluctuations  and 
possible  alignment  changes,  a  single  sample  of  NaCl  was 
tested  with  each  alkali  halide.  It  was  readily  determined 
that  a  slight  misalignment  of  the  focusing  lens  had  no 
measurable  effect  on  the  relative  breakdown  strengths. 

Visual  inspection  and  the  breakdown  statistics  suggested 
that  soatial  inhomogeneities  from  inclusions  were  not  af¬ 
fecting  the  results  except  in  the  single  case  of  RbCl  (see 
Sect.  D  of  Chant.  3).  Damage  which  we  regarded  as  intrinsic 
consisted  at  each  damage  position  of  a  single  pointed  region 
which  began  at  the  geometrical  focus  and  extended  a  very 
short  distance  back  toward  the  laser,  increasing  in  cross- 
section  to  give  a  tear-drop  appearance.  A  typical  example  is 
indicated  in  Fig.  13  in  Chapt.  3.  In  RbCl,  on  the  other  hand, 
regions  with  low  breakdown  thresholds  consisted  typically  of 
one  or  more  spherical  voioo  randomly  distributed  about  the 
focus  (Fig.  14  in  ChaDt.  3).  A  number  of  points,  however, 
did  aopear  visually  to  have  intrinsic  damage  and  were  con¬ 
sistently  more  difficult  to  breakdown.  These  data  points  were 
used  for  the  RbCl  results. 

Finally,  a  fast  photodiode  detector  system  with  a  C.5  nsec 
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rlsetime  monitored  the  transmitted  light  as  ahown  in  Fig.  i? 
and  was  used  to  confirm  threshold  levels  in  NaCl  and  KC1  as 
well  as  to  establish  the  approximate  time  structure  and 
stability  of  the  laser  output. 

Values  for  the  breakdown  field  obtained  at  1.06  |/ni  are 
summarized  in  Table  TV  and  in  Fie;.  18  alone:  with  both  the 
10.6  pm  data  collected  by  Yablonovitch1  ^  and  a*.  j  >pted  dc 
results.^  These  results  are  normalized  to  the  respective 
values  of  field  necessary  to  damage  NaCl  listed  in  Table  V. 
This  allows  the  striking  similarity  in  trends  of  breakdown 
field  to  be  easily  observed  and  the  possible  systematic  de¬ 
viations  at  1.06  jim  to  be  recognized.  The  quoted  errors  at 
10.6  jim  are  ±10  percent,  and  our  random  experimental  errors 
in  relative  fields  are  estimated  to  be  no  more  than  ±  10  per¬ 
cent  with  possible  errors  due  to  microscopic  strains  adding 
another  ±  5  percent.  Two  different  samples  of  both  NaCl  and 
KBr  from  two  different  manu.  _*•:  curers  a;avo  nearly  Identical 
results. 

The  breakdown  strength  of  NaCl  was  also  measured  with  a 
ruby  laser  to  confirm  the  absence  of  self- focusing  as  noted 
above.  Table  V  records  the  average  of  about  50  damage 
measurements.  Although  the  laser  was  normally  operating  in  a 
single  longitudinal  mode  as  indicated  by  Fabry-Perot  and 
photodiode  studies,  each  laser  shot  during  the  measurement  was 
monitored  with  a  fast  photod'iode  and  recorded. 
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TABLEIV 

RELATIVE  BREAKDOWN  FIELDS  -  NORMALIZED  TO 
ENfln«2  X  106  V/cm 


Nal 

NaBr 

NaCI 

NaF 

DC 

0.460 

0.553 

1 

1.60 

10. 6/im 

0.405 

0.476 

1 

1.06/um 

(0. 29)* 

0.67 

1 

1.64 

K 1 

KBr 

KCI 

KF 

DC 

0.380 

0.460 

0.667 

1.27 

10. 6^m 

0.369 

0.482 

0.713 

1.23 

l.Ofyzm 

0.27 

0.38 

0.57 

Rbl 

RbBr 

RbCI 

DC 

0.327 

0.387 

0.553 

10. 6/im 

0.323 

0.400 

0.472 

1.06/jm 

0.40 

0.55 

0.67 

Crystal  was  extremely  hydroscopic  and  no  final  check  was  made  with 
the  microscope  to  determine  if  inclusions  were  responsible  for  the 
damage  observed. 
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TABLE  V 


ABSOLUTE  BREAKDOWN  STRENGTH 
OF  NaCl 


E  .  (dc) 
paak 


X  10b  V /  cm 


Erms  (10'^m) 


(  1.  95  ±  0.  20)  X  10  V/cm 


E  (1.06  pm) 
rms 


(2.  3  ±  0.46)  X  I0b  V/cm 


E  (0.69  pm) 
rms 


(2.  2  ±  0.44)  X  10b  V/cm 


These  values  are  taken  from  Refs.  54  and  14 
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C.  Discussion  of  Results 

The  experiments  reporter  here  were  performed  under  care¬ 
fully  controlled  conditions  using  stable,  well-characterized 
lasers  and  optical  systems  for  which  aberrations  were  unim¬ 
portant.  Because  we  were  able  to  probe  each  sample  in  many 
different  positions,  random  fluctuations  in  breakdown  strength 
were  averaged  out.  It  was  possible  to  distinguish  between 
inclusion  and  intrinsic  damage  by  inspection  of  the  residual 
damage  sites  and  to  correct  for  the  effects  of  inclusions  in 
the  one  material  for  which  they  were  important.  In  addition, 
experimental  tests  showed  that  catastrophic  self-focusing  was 
absent  and,  consistent  with  theory,  that  the  index  nonlinearity 
did  not  affect  the  results  to  within  experimental  error.  It 
is  therefore  concluded  that  the  results  of  the  1 . 06  and  0.69^  m 
study  as  summarized  in  Fig.  18  and  Table  V  represent  accurate 
measurements  of  intrinsic  bulk  damage. 

Because  the  techniques  of  this  study  are  virtually  iden¬ 
tical  to  those  of  reference  14,  direct  comparison  can  be  made 
to  breakdown  strengths  at  10.6  #im.  It  has  already  been  ob¬ 
served  that  the  damage  thresholds  for  the  alkali-halides  at 
1,06  fim  follow  a  trend  nearly  identical  to  that  observed  with 
the  CO2  laser  and,  in  fact,  to  the  dc  measurements  of 
reference  54.  It  thus  appears  that  the  intrinsic  process  of 
laser-induced  damage  for  the  alkali-halides  has  the  same 
fundamental  character  as  both  ac  damage  in  the  infrared  and 
dc  avalanche  breakdown.  Moreover,  the  consistency  of  the 
optical  breakdown  strengths  at  0,69  |iKi  suggests  that  this 
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same  process  may  dominate  up  to  frequencies  approaching 
4  x  1014  herz. 

Data  from  Fie:.  18  and  Table  have  established  the 
relationship  (Elto6^rms  about  1,5  x  E^  for  ten  different 
compounds.  The  precise  value  of  the  factor  1.5  is  not  im¬ 
portant  since  dc  measurements  are  known  to  be  somewhat  sen¬ 
sitive  to  experimental  techniques. ^  It  is  important,  on 
the  other  hand,  that  consistent  measuring  techniques  have 
measured  roughly  the  same  factor  of  1.5  for  all  ten  alkali 
halides. 

Additional  support  for  an  avalanche  mechanism  comes  from 
three  experimental  observations  concerning  the  time-structure 
of  the  laser  probe  pulses.  The  first  is  that  increasing  tne 
pulse  width  of  the  YAG  laser  outout  by  a  factor  of  2.3  re¬ 
sulted  in  a  14  percent  average  drop  in  threshold  intensity 
for  NaCl.  Averages  were  taken  at  about  20  shots  at  each 
pulse  width.  This  change,  though  small,  is  probably  real, 
because  the  test  was  made  on  a  single  sample  of  high-quality 
NaCl  and  thereby  avoided  a  major  source  of  experimental  un¬ 
certainties  arising  from  material  variations.  The  second 
observation,  noted  at  both  1.06  and  0.69  pm,  is  that  high- 
frequency  time-structure  on  the  pulse  has  little  measurable 
effect  on  the  breakdown  strength.  And  finally,  after  adjusting 
the  power  level  so  that  damage  occurred  regularly  near  the 
top  of  the  laser  pulses,  the  probe  intensities  were  increased 
by  a  factor  of  about  three  by  changing  the  beam  attenuation. 
When  this  was  done,  the  intensity  at  which  the  transmitted 
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light  droDoed  (see  Fie.  l6)  was  higher  by  25  oercent  or  more 
than  it  had  been  with  the  lower  intensity  pulses.  This  was 
interpreted  to  mean  that  increasing  the  effective  risetime  of 
the  optical  field  raises  the  measured  breakdown  strength.  To 
understand  both  this  set  of  observations  and  the  results  from 
Table  v,  some  discussion  of  existing  electron  avalanche 
theories  is  given.  This  discussion  will  be  extended  in 
Chapt.  0. 

An  electron  avalanche  in  solids  is  a  raoid  multiplication 
of  conduction-band  electrons  in  which  an  initially  low  density 
NQ  of  free  carriers  interacts  with  an  intense  electric  field 
in  the  presence  of  phonons.  The  number  of  electrons  increases 
with  time  as 


N(t) 


N  exp 


'  t 

r  a(E)  dt 
0 


(4-1) 


The  gain  coefficient  a(E)  is  a  strongly-varying  function  whose 

value  can  be  inferred  from  dc  measurements  of  breakdown 

strength  as  a  function  of  sample  thickness  for  extremely  thin 

specimens.  For  Eq .  (4-1)  to  be  valid  the  rate  at  which 

electrons  are  lost  by  trapping  and  diffusion  out  of  the  focal 

volume  must  be  small  compared  to  the  rate  at  which  they  are 

generated.  For  Q-switched  laser  pulses,  the  electron  losses 

1  8 

are,  in  fact,  negligible. 

Two  important  conclusions  develop  from  such  an  analysis. 
The  first  is  that  the  entire  process  of  avalanche  and  damage 
involves  energy  exchange  between  the  field  and  the  electrons 
which,  as  has  already  been  discussed,  is  approximately 
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descrlbed  by  the  well-known  formula  for  ac  conductivity  given 
by  Eq.  (1-1).  Eq.  (1-1)  shows  that  the  energy  Input  to  the 
electrons  scales  with  frequency  and  field  as  E  /(I  +w'rc), 
and  because  the  details  of  energy  input  determine  the  electron 
distribution  function  and  hence  N(t),  the  threshold  for  damage 
will  scale  in  the  same  manner.  This  justifies  the  use  of 
root-mean-square  fields  in  Table  V,  It  also  indicates  that 
the  ac  breakdown  strength  will  increase  for  frequencies  near 
l/r .  An  estimate  of  r  for  NaCl*^  indicates  that  frequency 
dispersion  should  begin  to  occur  somewhere  near  that  of  the 
ruby  laser. 

The  second  relevant  conclusion  from  an  analysis  of  ava¬ 
lanche  breakdown  is  that  if  insufficient  time  exists  for  the 
electron  density  to  reach  values  necessary  to  damage,  then 
damage  will  not  occur  even  though  electron  recombination 
losses  have  been  exceeded  and  an  electron  avalanche  is  under¬ 
way.  Damage  will  only  be  produced  when  the  field  is  sub¬ 
sequently  raised  above  the  steady-state  breakdown  field  and 
the  gain  coefficient  a  is  correspondingly  increased.  Such  an 
effect  has  been  observed  in  dc  experiments  by  varying  the 
thickness  of  thin  (  <  50 n  m)  samples.^0,  In  such  experiments 

£  Q 

transit  effects  limit  the  build-up  time.7 

The  field  dependence  of  the  gain  coefficient,  or  ioni¬ 
zation  rate,  has  been  measured  experimentally  in  NaCl  by  using 
subnanosecond  laser  pulses.  This  measurement  will  be  pre¬ 
sented  in  Chapt.  5*  For  present  purposes  we  merely  note  that 
our  observations  of  a  pulse-width  dependence  to  breakdown,  the 
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insensitivity  of  threshold  to  fast  time  structure,  and  the 

Increase  in  breakdown  strength  for  rapidly  rising  pulses  all 

indicate  that  the  intrinsic  damage  process  is  time-dependent 

as  expected  from  an  electron  avalanche  and  that  the  process 

has  an  effective  time  constant  on  the  order  of  1  ns.  Such 

temporal  considerations  may  explain  part  of  the  difference  in 

damage  field  between  dc  and  laser  measurements  as  summarized 
62 

in  Table  V.  A  more  likely  explanation  for  most  of  this 
difference  involves  space  charges  and  will  be  discussed  in 
Chapt.  5. 

The  essential  details  of  the  breakdown  studies  reported 
here  are  thus  fully  explained  by  the  process  of  avalanche 
breakdown. 

In  summary,  careful  measurements  of  laser-induced  bulk 
damage  have  been  made  in  nine  alkali-halides  without  the 
confusing  effects  of  self- focusing.  Comparison  of  the  results 
to  studies  at  dc  and  at  10,6  pm  indicate  that  the  process  op 
ac  avalanche  breakdown,  similar  in  fundamental  character  to 
dc  avalanche  breakdown  is  responsible  for  the  damage  observed. 
Analysis  of  time-related  observations  confirms  this  conclusion. 


CHAPTER  5 


THE  EFFECTS  OF  LASER  FREQUENCY,  PULSE  DURATION 
AND  CRYSTAL  DISORDER  ON  INTRINSIC  OPTICAL  DAMAGE  FIELDS 

A,  Introduction 

In  the  last  chapter  the  intrinsic  breakdown  mechanism  was  identified 

in  the  alkali  halides  for  I.  06  m  radiation.  This  process,  an  electron 

avalanche,  has  been  studied  at  dc  for  over  forty  years.  Despite  the 

extensive  investigation  avalanche  breakdown  has  attracted,  very  little 

progress  has  been  made  in  understanding  the  microscopic  details  of  its 

development.  Not  only  is  the  proper  modeling  of  tnis  highly  complex 

process  beyond  the  present  level  of  solid  state  theory,  ^  but  experiments 

on  dc  breakdown  are  difficult  to  conduct  and  often  yield  very  limited 
58 

information.  No  estimates  can  be  obtained,  for  example,  of  effective 

electron  collision  rates  from  dc  experiments,  and  the  time  dependence  of 

1 8 

the  avalanche  can  only  be  studied  indirectly.  Laser  damage  techniques, 
on  the  other  hand,  are  comparatively  simple,  and  the  great  versitility  of 
the  laser  can  be  exploited  to  probe  aspects  of  an  avalanche  that  are  in¬ 
accessible  to  dc  investigators. 

In  this  chapter  we  extend  the  damage  measurements  in  the  alkali 
halides  to  higher  frequencies  and  shorter  pulse  widths.  It  will  be  shown, 
in  particular,  that  frequency  dispersion  begins  to  develop  at  0.  69  jyim  and 
that  the  damage  field  in  NaCl  increases  as  expected  with  decreasing  pulse 
width.  These  two  observations  provide  an  experimental  estimate  of  the 
high-field  electron-phonon  collision  frequency  and  an  experimental 
measure  of  the  avalanche  ionization  rate.  An  apparent  anomaly 
in  the  frequency  dispersion  suggests  the  possibility 


Preceding  page  blank 
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that  deep  lying  exciton  levels  may  affect  the  breakdown 
strength  of  NaF.  A  comparison  of  the  optical  frequency  ionl- 
zation  rate  to  estimates  based  on  dc  data1  show  at  least 
qualitative  agreement.  The  results  of  a  short  experiment  will 
also  be  described  in  which  the  effects  of  lattice  disorder  on 
the  intrinsic  breakdown  strengths  of  materials  were  studied. 

B.  Avalanche  Breakdown  Induced  by  Ruby  Laser  Light 

We  have  completed  measurements  of  intrinsic  bulk  break¬ 
down  in  nine  single  crystal  alkali  halides  using  a  TEMqq, 
single-longitudinal  mode  ruby  laser.  Self-focusing  was  absent 
in  these  studies,  and  damage  from  inclusion  absorption  was 
distinguished  from  Intrinsic  damage.  It  was  found  that  at 
O.69  pm  the  relative  breakdown  strength*  of  the  alkali  halides 
have  begun  to  differ  from  values  obtained  at  1 . 06  and  10.6  pm 
and  at  dc.  The  onset  of  this  frequency  dispersion  in  the 
avalanche  breakdown  orocess  enables  one  to  estimate  the  high- 
field  electron-phonon  collision  frequency. 

The  laser  system  and  techniques  for  avoiding  and  con¬ 
firming  the  absence  of  self-focusing  have  been  described  in 
detail  in  Chapts.  2  and  3.  Damage  from  inclusions  was  dis¬ 
tinguished  from  intrinsic  damage  by  examining  both  the 
morphology  of  the  damage  sites  and  the  temporal  shape  of  light 
pulses  transmitted  through  the  sample.  (See  Sect.  D  of  Chapt. 
3.)  The  latter  technique  employed  the  fact  that  a  damaging 
light  pulse  is  attenuated  in  a  manner  which  is  characteristic 
of  the  cause  of  damage.  Only  data  obtained  from  intrinsic 
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damage  events  were  considered  in  the  present  work. 

To  measure  the  breakdown  strengths  of  the  alkali  halides, 
we  focused  the  laser  beam  inside  the  sample  about  1.5  mm  from 
the  front  surface.  The  focal  spot  diameter  was  about  15  pm. 
When  the  incident  power  was  sufficiently  high,  a  bright  white 
spark  was  produced  and  a  small  volume  of  the  material 
(  <  10"®  cm^)  was  melted.  Because  a  well-defined  threshold 
could  not  be  found  (see  Chapt.  6)  we  again  defined  the  damage 
field  as  that  value  of  rms  electric  field  inside  the  material 
which  was  necessary  in  order  to  Induce  damage  on  a  simple 
oulse  with  a  probability  of  0.5.  The  incident  power  necessary 
to  reach  this  field  in  NaP  under  our  conditions  of  focusing 
was  146  KW.  For  the  other  crystals,  the  input  power  was 
between  20  KW  (for  NaCl)  and  9.2  KW  (for  Rbl). 

A  large  numbei  of  data  points  were  taken  for  each  sample 
(40  to  100),  but  normally  less  than  half  of  these  damage 
sites  could  be  unambiguously  identified  as  resulting  from 
intrinsic  processes.  The  remainder  were  apparently  caused 
by  the  presence  of  absorbing  inclusions. 

,  As  in  Chapt.  4,  the  rms  on-axis  damage  field  was  deter- 

[  ( 

[  !  mined  for  NaCl  using  the  measured  intensity  distribution  at 

i 

the  focus  of  the  14  mm  focal-length  lens.  This  measurement 
was  the  basic  calibration,  and  in  order  to  avoid  errors  from 
possible  alignment  chane-es,  damage  fields  for  the  other 

i  1 

I  materials  were  measured  relative  to  E^aC^.  The  same  sample 

t  * 

f  of  NaCl  was  tested  with  each  alkali  halide. 

[  | 

|  Figure  19  and  Table  V  (Chapt.  4)  summarize  the  results 

i 
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of  this  study,  the  1.0  6  p  m  data  from  Chapt.  4,  and  accepted 
dc  results.  ^  In  Fig.  19  the  damage  fields  of  the  various 
alkali  halides  are  normalized  to  that  of  NaCl  at  the  appro¬ 
priate  frequency.  This  makes  the  variation  in  trends  between 
the  1.06  urn  and  the  0.69  pm  measurements  more  easily  seen.  The 
relative  breakdown  strengths  at  1.06  pm  are  virtually  identical 
to  those  measured  at  dc.  Although  the  corresponding  data  at 
10.6  are  not  displayed,  they  too  follow  the  same  trend. 

This  is  not  the  case  at  0.69  pm,  however,  even  when  the  10  to 
15  percent  measurement  errors  are  considered. 

Table  V  lists  the  rms  damage  fields  of  NaCl  for  these 
experiments.  The  agreement  found  in  the  four  experiments  is 
heartening  because  there  are  oarticular  difficulties  in 
determining  the  absolute  damage  fields  in  dc  measurements.^®'  ^ 
Root-mean-square  values  of  the  electric  field  strength  are 
given  because  for  laser  pulses  of  this  intensity,  the  build¬ 
up  time  of  an  electron  avalanche  to  damaging  levels  is  on  the 
a  4 

order  of  10 ;  to  10  cycles  of  the  optical  field.  (See  next 
Section.)  Heating  of  the  electron  population  is  thus 
effectively  averaged  over  many  cycles. 

Because  the  measured  breakdown  strength  of  NaCl  was  the 
same  within  experimental  error  at  1.06  pm  and  0.69  pm,  the 
ordinate  in  Pig.  19  could  be  changed  to  read  in  absolute  field 
units  with  virtually  no  change  in  the  relative  positions  of 
the  1,06  pm  and  0,69  pm  data  points.  Absolute  field  units 
were  not  used,  because  the  resulting  error  bars  on  the  altered 
plot  would  be  much  larger.  They  would  include  net  only  the 
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experimental  uncertainty  ir»  field  value*?  relative  to  N«tC  1 ,  as 
displayed  in  Fhr.  19.  but  they  would  also  include  the  coro-- 
rara  lively  lar^-e  experimental  uncertainty  ir.  the  absolute 
damane  fie'd  for  Nad  hr  suwnari  ted  in  "able  V,  It  is  the 
c  ha  nee  in  trends  of  damage  f  \  *  3  d  ?  which  is  most  T.eeninw  ful , 
anc  we  ha.ve-  chosen  to  display  the  data  and  experimental  ur- 
c e r  fa'nfv  in  a  rearm*' r  w.ni  co  will  hart  under root  e  xn.is  i/hnnce. 

A  cc.v  tete  theoretical  destr  r.'tior.  of  avalanche  break¬ 
down  has  not  been  dev*  looe-i ,  T^e  onset  of  frequency  disper¬ 
sion  can  be  u  u*a  1 1  ta  ti  ve  1  y  understood,  however,  by  us  in*  the 
results  jf  models  baaed  on.  relaxation  t  •  me  aprroxi  ^ati  ans 
which  ■ve  oa\  o  brie:!;/  conn  i  d^rc-j  ir.  Ch-vpls.  \  and  u  and  which, 
we  will  consider  l  r  rr.ve  aet-'.il  in  Ch  "r>t.  &.  These  models 
predict  that  the  Harare  field  will  seal  e  with  frequency  as 
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where  o  =  is  t;»e  laser  ridlar.  frequency  and  r  is  an 

effective,  hi  rh-  field  ,  electron-- Phonon  collision  time.  Since 
«gcf.  me  ter  {ml  h~..n  a  di  ffer  valv*  of  ••  ,  the  relative  break¬ 
down  fields  for  the  alkali  hal i de?  should  beedn  to  change  s  t 

hi  rK  f  ^ «. i  j  *  i.  ;•  1 1—  ,  O' .  '  ft  iij  h  c'lr^r  i  ri  ^ 1  i*y**  19  Can  be 

used  to  1  rfer  aonre >: i  ra  f e  rela  ti "e  vs i u e s  of  r  for  a  variety 
o  f  t ho ore 1 1  c-a  ]  iron  e Is  . 
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tud trial  or';:  cal  photon  erercies.’'4 
tr.e  results  o*  there  calculations 
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the  change  in  relative  breakdown  strengths  observed  at  ruby 
frequencies  for  most  of  the  alkali  halides.  The  collision 
time  for  NaCI  under  this  assumption  is  about  2  x  lO”^  sec. 

For  NaF  these  calculations  predict  that  the  relative  break¬ 
down  field  (ENaF^NaCl  19)  will  decrease  at  ruby 

frequencies,  contrary  to  the  change  which  is  experimentally 
observed.  This  discrepancy  may  be  the  result  of  the  inade¬ 
quacy  of  the  perturbation  calculations  or  it  may  indicate  that 
the  frequency  dependence  of  the  electron  avalanche  is  not 
determined  by  the  electron-phonon  collision  frequency  alone. 

Seitz^  hag  suggested  that  the  presence  of  deep-lying 
exciton  bands  may  influence  the  dielectric  strength  of  alkali 
halide  crystals.  If  this  is  the  case,  then  as  the  field 
frequency  cj  is  increased,  direct  excitation  out  of  these  bands 
becomes  possible  and  the  damage  field  will  decrease.  NaF, 
which  has  the  deepest  lying  bands  (1.5  -  2,0  eV)  of  the 
materials  studied, ^  will  experience  this  effect  at  a  higher 
frequency  than  the  other  aj-kali  halides.  Such  considerations 
of  the  relative  importance  of  the  exciton  bands  may  explain 
the  observed  large  increase  in  the  relative  NaF  damage  field 
at  0.69  M  m. 

Multiphoton  absorption  directly  across  the  bandgap  can¬ 
not  explain  the  changes  in  relative  breakdown  strength  which 
have  been  observed.  In  addition,  theoretical  calculations  of 
the  fields  at  0.69  necessary  to  induce  damage  from  multi¬ 
photon  ionization  or  from  its  low-frequency  limit,  tunnel 
ionization,  give  damage  fields  which  are  about  an  order  of 
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magnitude  larger  than  those  measured.  ’  ~ 

In  conclusion,  we  have  measured  the  intrinsic  optical 
breakdown  fields  of  nine  alkali  halides  using  a  ruby  laser. 
Although  the  absolute  and  relative  thresholds  fields  are  com¬ 
parable  to  the  thresholds  observed  at  1.06  pm  and  at  dc, 
differences  are  observed  which  suggest  that  at  v  =  4.3  x  lO1^ 
sec”1,  avalanche  breakdown  is  no  longer  identical  to  dc  di¬ 
electric  breakdown.  Current  theories  of  avalanche  breakdown  do 
not  appear  to  explain  the  details  of  this  observed  difference. 

C.  The  Pulse-Width  Dependence  of  Optical  Avalanche  Breakdown 

Measurements  are  reported  here  of  optical  damage  induced 
by  subnanosecond  laser  pulses.  These  measurements  were  per¬ 
formed  by  focusing  mode-locked  YAG«Nd  laser  pulses  having 
durations  of  1 5  and  300  picoseconds  inside  single  crystal 
NaCl.  Because  the  experimental  procedures  used  in  the  present 
work  were  identical  to  those  used  in  the  studies  of  Chapt.  4 
with  a  Q-switched  YAGiNd  laser,  the  subnanosecond  measurements 
can  be  directly  compared  to  the  results  of  those  studies.  It 
was  found  that  the  intrinsic  breakdown  field  increased  by 

7 

almost  an  order  of  magnitude  to  over  10  volts/cm  as  the  laser 
pulsewidth  was  decreased  from  10  ns  to  15  ps.^  The  de¬ 
pendence  of  the  damage  field  on  laser  nulse  duration  is  used 

to  calculate  a  field-dependent  ionization  rate  which  is  com- 

1  8 

oared  to  the  predictions  of  Yablonovitch  a.'d  Bloembergen 
who  estimated  the  ionization  rate  from  Published  measurements 
of  the  dc  dielectric  strength  of  NaCl  for  thin  samples  with 
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varying  thickness.  Qualitative  agreement  is  found. 

The  laser  used  for  the  present  work  (9ee  Fig.  12  in 
Chaot.  3)  was  a  passively  mode-locked  YAGiNd  laser  operating: 
in  a  TEMqo  mode  at  1.06  n m.  Without  intercavity  etalons,  this 
oscillator  produced  bandwidth-limited  light  pulses  of  15 
picosecond  duration.  When  the  cavity  output  mirror  was  re¬ 
placed  with  a  sapphire  etalon,  the  oulsewidth  was  lengthened  to 
about  300  picoseconds.  Two-photon-fluorescence  measurements 
failed  to  detect  substructure  with  pulses  of  either  duration. 

A  laser-triggered  spark  gap  was  used  to  select  a  single  light 
pulse  which,  after  attenuation,  was  focused  through  a  14  mm 
focal  length  lens  about  2  mm  into  the  sample.  Care  was  taken 
to  insure  that  spherical  aberrations  from  both  the  lens  and 
the  Plane  entrance  surface  of  the  sample  being  tested  were 
unimportant.  An  energy  monitor  recorded  the  energy  in  each 
laser  Dulse. 

Although  the  intrinsic  damage  process  in  transparent 

12 

materials  is  an  inherently  statistical  process  (see  Chapt.  6), 
it  is  virtually  threshold-like  in  NaCl.  Consistent  with  Chapt. 
b •  a  damage  field  can  be  defined  as  that  value  of  root-mean- 
square  electric  field  inside  the  sample  which  produced  damage 
on  a  single  shot  with  a  probability  of  0.5.  Damage  was  iden¬ 
tified  by  the  occurrence  of  a  faint  spark  and  was  accompanied 
by  a  small  melted  region  (<  2  x  10“9  cm^)  inside  the  crystal. 

At  least  20  data  points  were  taken  for  each  pulse  duration  at 
the  0.5  Probability  point. 

Beam  distortion  from  self- focusing  was  avoided  by 
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confining  the  laser  inDut  cowers  to  well  below  the  calculated 
critical  powers  for  eatastronhic  self-focusing.  (See  Table 
VI.)  To  verify  the  absence  of  self- focusing,  two  different 
lenses  (focal  lengths  of  14  mm  and  2  5  mm)  were  used  to  focus 
the  laser  radiation.  As  expected  from  diffraction  effects 
alone,  the  input  damage  powers  scaled  with  the  square  of  the 
focal  lengths.  If  eatastronhic  self- focusing  had  been  present 
with  the  subnanosecond  pulses,  the  input  damage  cower  would 
have  been  independent  of  focal  length.  The  absence  of  in¬ 
clusion  damage  was  confirmed  in  the  manner  discussed  in  Sect.  D 
of  ChaDt.  3. 

Table  VI  summarizes  the  results  of  the  present  measure¬ 
ments  and  those  of  ChaDt.  4.  An  increase  in  breakdown 
strength  was  observed  as  the  duration  of  the  laser  pulse  was 
decreased.  As  the  pulse  duration  was  changed  from  10.3  ns 
to  15  PS,  there  was  a  total  change  by  a  factor  of  5.8  in 
damage  field  strengths  or  a  factor  of  33  in  damage  intensity. 
The  exoerimental  points  are  plotted  in  Fig.  20  along  with  the 
semi-empirical  predicted  curves  from  Ref.  18. 

The  existence  of  a  pulse-width  dependence  to  intrinsic 

damage  is  qualitatively  explained  by  classical  theories  of 

67—68 

electron  avalanche  ionization.  *  "  These  theories,  which 

are  summarized  in  Chant.  8,  predict  that  the  density  of  con¬ 
duction-band  electrons,  N(t),  increases  with  time  as 
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Experimental  Breakdown  Fields  and 
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’ig.  20.  THE  RATE  OF  IONIZATION  IN  NaCl  AS  A  FUNCTION  OF 
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1  8 

Eq,  (5-2)  is  valid  when,  as  is  the  case  with  laser  breakdown, 
electron  diffusion  and  trapping  can  be  ignored.  N0  is  the  low 
density  of  conduction  electrons  before  the  application  of  the 
electric  field  and  a  is  the  ionization  rate  which  increases 
monotoni cally  with  increasing  electric  field.  Breakdown 
occurs  when  the  density  of  electrons  becomes  high  enough  to 
cause  a  material  irreversibility  such  as  a  phase  change.  (See 
Chapt.  1.)  As  the  time  available  for  the  avalanche  to  develop 
to  damaging  proportions  decreases,  the  rate  of  ionization  and 
hence  the  electric  field  must  be  increased  in  order  to  produce 
damage. 

It  is  desireable  to  compare  the  laser  data  to  dc  results. 
Such  a  comparison  cannot  be  made  directly,  because  impulse 
dc  measurements  with  subnanosecond  impulse  durations  have  not 
been  made.  As  Yablonovitch  and  Bloembergen1  have  suggested, 
however,  dc  measurement!*  on  samoles  with  varying  thickness 
provide  an  indirect  comparison  because  the  maximum  duration 
of  the  dc  avalanche  is  limited  to  the  electron  drift  time  from 

cQ 

cathode  to  anode,  Ey  considering  limits  on  the  electron 
drift  velocity,  Yablonovitch  and  Bloembergen  have  calculated 
a(E)  for  dc  fields  using  Eq,  (5-2)  and  previously  published 
measurements  of  dc  damage  fields  in  thin  samples  of  NaCl.^0, 
Breakdown  was  assumed  to  occur  when  in  Eq .  (5-2)  reached  a 

Q 

value  of  10  .  The  ionization  rates  for  the  laser  data  can  be 
found  by  replacing  the  integral  in  Eq.  (5-2)  by  a(Erms)tp 
where  Ermg  is  the  root-mean-square  field  on  axis  at  the  peek 
of  the  laser  pulse  and  tp  is  the  laser  pulsewidth.  Then 
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a (firms )  is  given  by 

a(Erms)  =  T~  £nMc  *  18/tn  '  (5-3) 

P  p 

This  relation  has  been  U3ed  to  convert  the  quantity  a(E)  used 
along  the  vertical  axis  in  the  figure  of  Ref.  18  to  our  figure 
which  used  tD“*.  We  have  shifted  the  curves  along  the  hori¬ 
zontal  axis  to  obtain  agreement  with  the  experimental  values 
for  the  breakdown  field  firmg  for  the  long  pulses. 

It  should  be  noted  that  the  damage  fields  for  the  dc 
studies  of  Refs.  60-6l  are  approximately  a  factor  of  2.3  lower 
than  the  corresponding  values  for  laser-induced  damage.  There 
are  two  reasons  to  suspect  that  this  factor  is  a  systematic 
error  in  determining  absolute  field  strength  at  dc  rather  than 
an  indication  of  a  fundamental  difference  between  dc  and  laser- 
induced  breakdown. First,  this  factor  of  2.3  is  nearly  the 
same  for  all  nine  alkali  halides  studied  at  10.6  ^m1*4  and 
1.06  Jim,  and  second,  dc  field  values  are  average  values  given 
by  the  voltage  difference  between  anode  and  cathode  divided 
by  their  separation  without  regard  for  field  inhomogeneities 
from  effects  such  as  space  charges  which  can  be  important  in 
dc  experiments.  Seme  evidence  in  fact  exists  in  the  litera¬ 
ture  that  sDace  charge  effects  are  influencing  dc  breakdown. 

For  example,  dc  impulse  experiments  have  been  performed^0  in 
which  the  duration  of  the  applied  field  was  reduced  to  about 
10  ns,  a  time  interval  too  short  to  allow  the  development  of 
ionic  space  charges.^1  It  was  found  that  the  damage  field  in 
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NaCl  increased  from  about  1  MV/cm,  its  steady-state  value,  to 
.lust  over  2  MV/cm,  a  value  which  agrees  closely  with  the  laser 
measurements.  The  time-related  considerations  discussed  in 
this  section  do  not  explain  this  change.  In  other  experiments 
the  start  of  an  avalanche  has  been  observed  to  vary  in  a 
statistical  manner  with  average  lags  of  1  |is  or  more.''’2'  ^ 

This  effect  has  been  interpreted  as  resulting  from  space  charge 
development.^  A  complete  understanding  of  space  charge 
effects,  electrode  effects,  and  other  factors  which  can  affect 
dc  measurements  of  breakdown  has  not  emerged,  and  thus  we 
should  regard  absolute  measurements  of  dc  breakdown  with  some 
caution,  preferring  instead  to  compare  trends  in  absolute 
field  strengths  as  different  materials  are  investigated  or  as 
parameters  are  varied. 

In  Fig.  20  the  four  laser  measurements  are  plotted  with 
the  computed  curve  from  Ref.  18.  The  two  branches  to  the 
computed  curve  correspond  to  two  limits  on  the  high-field 
electron  drift  velocity.  Within  experimental  error,  the  laser 
data  overlap  the  upper  curve  of  Ref.  l8  which  was  derived  on 
the  assumption  that  the  mobility  in  the  hot  electron  gas  is 
independent  of  Erms.  Quantitative  agreement  should  not  be 
emphasized,  however,  because  the  present  analysis  is  based  on 
at  least  two  important  assumptions  which  may  not  be  valid 
over  the  range  of  damage  fields  considered.  The  first  assump¬ 
tion  is  that  factors  in  the  dc  experiments  such  as  space 
charges  and  electrode  effects  do  not  change  as  the  sample 
thickness  is  reduced  to  approximately  a  micron.  And  the 
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second  assumption  is  that  the  same  intrinsic  mechanism 

dominates  over  the  range  of  laser  pulsewldths  in  Table  VI, 

Another  intrinsic  mechanism--multlphoton  ionisation  --may 

cause  damage  at  lower  fields  than  required  for  avalanche 

breakdown  when  the  laser  pulsewidth  is  extremely  short. 

Estimates  for  1.06  p m  radiation  in  NaCl  indicate  that  when 

the  laser  pulsewidth  is  less  than  about  a  picosecond,  multi- 

1 8 

photon  ionization  is  responsible  for  intrinsic  damage. 

Since  the  shortest  pulsewidth  considered  in  the  present  work 
is  15  ps ,  the  neglect  of  multiphoton  ionization  appears  to  be 
justified.  If  the  estimates  of  Yablonovitch  and  Bloembergen 
are  inaccurate,  however,  and  damage  from  multiphoton  ionization 
is  occurring,  the  ionization  rate  determined  from  the  1 5  ps 
pulse  is  an  upper  bound  for  the  actual  value  of  a  at  Erms  = 

12.4  MV/ cm. 

In  summary,  intrinsic  laser-induced  damage  has  been  shown 
to  be  a  time-dependent  process.  As  the  laser  pulsewidth  was 
decreased  to  15  ps,  the  damage  field  in  NaCl  increased  to 
over  10?  v/cm.  From  the  pulsewidth  dependence  of  the  optical 
damage  field,  a  field-dependent  ionization  rate  was  determined 
and  found  to  agree  at  least  qualitatively  with  experiments 
using  dc  fields.  The  agreement  underscores  the  basic  simi¬ 
larity  between  intrinsic  laser-induced  damage  at  1.06  pm  and 
dc  dielectric  breakdown. 

D.  Effects  of  Disorder  on  the  Intrinsic  Damage  Field 


Measurements  are  reported  here  of  optical  bulk  damage  in 
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three  disordered  system3--polycrystal  KC1,  a  single-crystal 
KBr-KCl  (6??S- 33d£)  alloy,  and  fused  quartz.  In  each  case  the 
damage  field  for  the  disordered  system  is  compared  to  the 
optical  strength  of  the  corresponding  crystal.  These 
measurements  were  made  in  order  to  determine  if  the  optical 
breakdown  field  increases  with  severe  lattice  disorder  as  had 

n  c  i  h 

been  observed  in  dc  breakdown  experiments  '  and  as  pre¬ 
dicted  by  simple  theories  of  avalanche  breakdown. 

The  laser  system  and  the  experimental  techniques  used 
here  were  identical  to  those  of  Chapt,  4  except  that  in¬ 
spection  of  the  transmitted  laser  light,  rather  than  micro¬ 
scopic  insoection  of  the  damaged  crystal,  was  used  to  dis¬ 
tinguish  inclusion  damage. 

It  was  found  that  the  damage  field  of  the  large-grain 
(20  pm)  polycrystal  was  the  same  as  that  measured  in  the 
single  crystal  and  that  the  alloy  damage  field  was  about  20 
percent  larger  than  the  damage  field  measured  in  the  majority 
constituent,  KBr .  In  the  quartz  system,  on  the  other  hand, 
the  disordered  (amorphous)  phase  was  noticeably  stronger  than 
the  crystal,  the  ratio  of  damage  intensities  being  5  ±  1. 

This  ratio  is  identical  to  the  corresponding  ratio  of  surface 
damage  fields  measured  by  Bass  and  Barrett, 

It  is  to  be  expected  that  the  large-grain  polycrystal 
should  have  the  same  damage  field  as  the  single  crystal.  The 
average  grain  diameter  (20  pm)  and  the  laser  focal  diameter 
are  comparable  so  that  in  the  high  intensity  region  near  the 
beam  axis  where  breakdown  is  observed  to  initiate  (see  Fig.  13 
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in  Chapt.  3),  the  sample  looks  like  a  single  crystal. 

By  a  simplified  argument  we  can  predict  the  approximate 
crystallite  size  necessary  to  affect  the  breakdown  strength. 
Classical  theories  of  avalanche  breakdown  (see  Chapt.  8) 
predict  that  the  dynamics  of  electrons  with  energies  greater 
than  the  longitudinal  optical  (LO)  energy  determine  the 
characteristics  of  the  avalanche.  The  LO  energy  in  the  alkali 
halides  corresponds  to  electron  momenta  of  about  0.1  times 
the  reciprocal  lattice  vector  G  so  that  the  important  electrons 
have  k  >0.1  G.  Phonons  with  values  of  q  >  0.1  G  will  inter- 

a/ 

act  most  strongly  with  these  electrons.  Because  such 
lattice  vibrations  have  wavelengths  equal  to  10  lattice  con¬ 
stants  or  less,  we  expect  that  unless  crystal  disorder  appears 
on  the  scale  of  about  10  lattice  constants  (  -  lOoif)  or  less 
the  damage  field  should  be  unaffected  by  disorder. 

Amorphous  systems  are,  of  course,  disordered  on  such  a 
scale.  Our  observation  that  the  glass  is  more  resistent  to 
damage  is,  therefore,  consistent  with  the  argument  .just  pre¬ 
sented.  The  behavior  of  the  quartz  system  can  be  explained 
in  somewhat  more  quantitative  terms  by  an  equivalent  argument. 
Referring  to  Eq.  (l-l)  on  p.  1-3*  it  is  seen  that  the  rate  of 
energy  Input  into  the  electron  population  decreases  with  de¬ 
creasing  electron  mobility  ji.  It  is  known  that  the  mobility 
of  disordered  systems  is  less  than  the  crystal  mobility'  so 
that  it  should  be  more  difficult  tn  heat  the  electron  dis¬ 
tribution  in  the  glass.  As  we  have  observed,  then,  the 
damage  field  for  glass  should  be  higher  than  the  crystal 
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damage  field. 

The  breakdown  field  of  the  KBr-KCl  alloy,  measured 
relative  to  the  breakdown  of  NaCl,  was  0.45  ±  0.06.  This 
value  comnares  to  0.38  for  pure  KBr ,  0.57  for  pure  KC1,  and 
0.44  for  the  average  breakdown  field  Eave.  weighted  according 
to  the  alloy  composition!  Eavff  =  0.6?  +  0.33  %ci* 

may  be  fortuitous  that  the  measured  breakdown  strength  of  the 
alloy  and  Eavg.  agree.  On  the  other  hand,  the  average  lattice 
constant  of  this  alloy  determined  by  X-ray  diffraction  scale 
with  comoosition,  and  other  material  parameter  such  as 
bandgap,  dielectric  constant,  and  reststrahlen  frequency 
probably  scale  in  the  same  linear  manner.  Since  the  break¬ 
down  field  depends  on  these  various  material  parameters , ^ 
it  is  reasonable  that  the  breakdown  field  should  also  scale 
with  composition — provided  that  the  crystal  disorder  resulting 
from  the  mixed  composition  does  not  affect  the  breakdown 
strength.  I^ore  alloy  systems  must  be  studied  before  firm 
conclusions  can  be  drawn  concerning  the  effects  of  alloying 
on  the  breakdown  strength  of  materials. 

It  thus  appears  that  only  extreme  lattice  disorder  such 
as  present  in  amorphous  systems  has  a  measurable  effect  on 
intrinsic  damage  fields.  Future  work  is  needed,  however,  to 
ascertain  any  general  correlations  between  lattice  disorder 
and  breakdown  fields.  Such  work  is  important  not  only  to  an 
understanding  of  avalanche  breakdown  but  it  is  also  important 
from  a  practical  viewpoint  to  the  design  of  more  damage- 
resistent  optics. 
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new  experimental  data  is  compatible  with  the  results  of 
Ref.  12.  Some  new  arguments  will  be  presented  in  i>ect.  C  to 
support  the  original  conclusion  that  the  surface  damage 
statistics  are  intrinsic  in  origin,  and  possible  sources  of 
statistics  within  the  avalanche  process  will  be  discussed. 


B •  Experimental  Measurements  of  Laser  Damage  Statistics 

1  2 

Bass  and  Barrett1  observed  that  a  precisely  defined 
threshold  for  laser  induced  surface  damage  does  not  exist  in 
the  ten  different  solids  which  they  studied.  Instead,  there 
is  a  range  of  power  levels  within  which  damage  can  develop  on 
each  shot  with  some  finite  probability  p1  such  that  0  <  p1  <  1. 
The  damage  probability  p^  at  some  power  level  was  defined  as 
the  ratio  cf  the  total  number  of  damage  sites  to  the  total 
number  of  laser  shots.  (See  Eq.  (3-1).)  By  varying  the  power 
level  of  the  highly-stabilized  laser  source,  they  were  able 
to  measure  p^  as  a  function  of  the  laser  power.  Fig.  21,  which 
summarizes  their  experimental  results,  shows  that  the  data 
appears  to  satisfy  the  relationship 


Pj  «  exp(-  K/E)  . 


(6-1) 


An  analysis  of  the  statistical  distribution  f^  (see  Eq .  (3-2)) 
confirmed  that  the  statistics  were  not  resulting  from  either 
material  homogeneities  fixed  in  position  or  from  laser 
fluctuations. 

Additional  evidence  for  the  probabilistic  nature  of 
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laser  damage  was  obtained  in  experiments  in  which  an  imago 
converter  streak  camera  was  used  to  measure  the  distribution 
of  breakdown  starting  times  for  surface  damage  to  two  different 
materials.  (The  time  at  which  a  visible  spark  first  appeared 
was  assumed  to  be  the  breakdown  starting  time.)  If  the  break¬ 
down  process  were  completely  described  by  a  well-defined 
threshold,  then  this  distribution  should  be  very  narrow.  The 
data,  however,  showed  broad  distributions,  particularly  when 
the  applied  field  was  such  that  p  <  1.  In  all  cases  the  most 
probable  starting  time  for  breakdown  occurred  before  the  time 
of  maximum  field. 

Damage  statistics  develop  because  there  exists  a  statis¬ 
tical  time  lag  to  the  damage  process  as  confirmed  by  Bass  and 
Barrett’s  streak  camera  measurements.  When  lasers  are  used 
to  induce  damage,  the  electric  field  is  applied  for  only  a 
short  time  interval  so  that  late-starting  avalanches  may  not 
be  able  to  develop  to  damaging  proportions  before  the  laser 
pulse  has  passed.  Avalanche  breakdown  induced  by  a  continuous 
laser  source,  on  the  other  hand,  should  have  a  well-defined 
threshold  for  damage  instead  of  a  statistical  threshold  as 
indicated  by  Eq.  (6-1), 

During  the  Q-switched  damage  experiments  on  the  alkali 
halides  (see  Chapts.  4  and  5).  evidence  was  recorded  which 
supports  the  statistical  viewpoint.  By  monitoring  the  trans¬ 
mitted  laser  light  with  a  fast  photodiode-oscilloscope  combin¬ 
ation  (risetime  =  0.5  ns),  we  observed  that  the  laser  light 
is  attenuated  when  damage  develops.  The  first  instant  of 
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measurable  attenuation  can  occur  before,  at,  or  after  the 
peak  of  the  laser  pulse,  so  that  no  well-defined  relationship 
exists  between  the  laser  intensity  and  the  first  instant  of 
attenuation.  Fig.  22  shows  examples  of  such  observations 
made  with  a  ruby  laser  beam  focused  to  produce  damage  inside 
NaCl.  The  laser  pulses  are  fully  time-resolved  as  verified 
by  Fabry-Perot  studies.  In  Fig.  23  another  ruby  laser  pulse 
was  focused  into  NaCl  but  did  not  cause  damage  whereas  a 
second  pulse,  apparently  identical  to  the  first,  did  induce 
damage  when  focused  into  the  same  volume  of  the  crystal. 

No thine  was  moved  between  the  two  laser  shots,  and  the  laser 
was  firing  automatically  at  a  repetition  rate  of  about 
1  pulse/5  sec,  A  similar  observation  was  recorded  at  1.06/im 
in  Fig.  17  of  Chapt.  U  where  the  automatic  firing  rate  was 
jus t  over  1  pps. 

The  relationship  between  laser  light  attenuation  and 

the  size  of  the  electron  avalanche  is  difficult  to  establish. 

A  reasonable  estimate  indicates  that  the  transmitted  light  is 

unaffected  by  the  avalanche  until  the  density  of  conduction- 

1 8  -3 

band  electrons  reaches  a  level  of  about  10A  cm  •  It  is  not 
strictly  correct,  therefore,  to  associate  the  first  instant 
of  attenuation  with  the  breakdown  or  avalanche  starting  time. 
But  as  noted  below,  there  should  be  essentially  no  random 
character  to  the  time  delay  between  the  avalanche  starting 
time  and  the  first  instant  of  attenuation.^  We  will  assume 
as  a  simplification  that  the  two  instants  of  time  are  identical 
within  experimental  error.  It  should  be  noted,  incidentally, 
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THE  OCCURRENCE  OF  INTERNAL  DAMAGE  IN  NaCI 
ATEMqq  mode  ruby  laser  with  total  pulse  energy 
of  0.3  mJ  was  focused  inside  the  inclusion  free 
sample  with  a  14mm  focal  length  lens 

The  laser  intensity  transmitted  tin  the  sample  is  shown  in  these  photographs. 

Fig.  22.  THE  OCCURRENCE  OF  INTERNAL  DAMAGE  IN  NaCI 
DUE  TO  RUBY  HAS  UR  IRRADIATION 
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10  nsec 


Ruby  laser  pulses  in  NaCI  with  no  self-focusing. 
The  first  pulse  caused  no  damage.  Five  seconds 
later  the  second  pulse  caused  damage  before 
peak  field  was  reached. 

>.  uout  5  seconds  passed  between  the  apparently  identic  al  laser  pulses.  Because 
nothing  was  moved  during  this  time  interval,  the  same  volume  of  material  was 
irradiated  by  both  pulses. 


Fig.  23. 


INTRINSIC  BULK  DAMAGE  OCCURRING  IN 
NaCI  ON  THE  SECOND  OF  TWO  RUBY 
LASER  PULSES 


Reproduced  Irom 
best  available  copy. 


-124- 


that  the  time  interval  over  which  the  transmission  drops  from 
effectively  unity  to  nearly  zero  may  have  no  simple  relation¬ 
ship  to  the  ionization  rate  a(E)  defined  in  Chapt.  5«  Because 
the  density  of  electrons  is  quite  high  by  the  time  measurable 
attenuation  occurs,  electron-electron  collisions  may  alter  the 
avalanche  development  and  reflection  from  the  plasma  may  be¬ 
come  important,  particularly  at  longer  laser  wavelengths. 

Encouraged  by  our  experimental  observations,  we  decided 
to  conduct  a  careful  measurement  of  the  statistical  time  lag 
in  a  number  of  materials.  Pigs.  24-26  summarize  meiasurements 
made  in  fused  quartz  at  1.06  pm,  in  NaP  at  0.69  pm,  and  in 
sapphire  at  1.06  pm.  As  usual,  self- focusing  was  found  to  be 
absent  under  the  conditions  of  our  measurements.  In  each 
case  30  to  50  damage  events  were  recorded,  and  wa3  defined 
by  Eq.  (3-1)'  The  position  of  the  peak  of  the  laser  pulse 
could  be  identified  to  within  intervals  of  0.5  ns  for  the 
1.06  pm  pulses  and  1.0  ns  for  the  0.69  p m  pulses.  Fig.  15  of 
Chapt.  3  and  the  techniques  of  Chapt.  3  (Sect.  D)  were  used  to 
verify  that  the  statistics  were  not  occurring  because  of  the 
presence  of  point-to-point  material  inhomogeneities  such  as 
inclusions. 

The  distributions  of  breakdown  starting  times  for  surface 
and  bulk  damage  to  fused  quartz  are  shown  in  Fig.  24  to  be 
virtually  identical  for  the  same  value  of  p^.  These  distri¬ 
butions  have  the  same  characteristics  as  those  reported  by 
Bass  and  Barrett  for  entrance  surface  damage.  The  breakdown 
can  begin  at  any  time  over  a  relatively  long  interval  which 
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24  MEASURED  DISTRIBUTIONS  OF  BULK  BREAKDOWN  STARTING  TIMES  FOR 
FUSED  QUARTZ  USING  YAG;Nd  LASER  IRRADIATION 


t(nsec) 

Fig.  25.  MEASURED  DISTRIBUTION  OF  BULK  AND  SURFACE 

BREAKDOWN  STARTING  TIMES  FOR  NaF  USING  RUBY 
LASER  IRRADIATION 
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Q-Switched  YAG:Nd  Laser 
Pulse  Waveform 


Sapphire 


Fig.  26.  MEASURED  DISTRIBUTION  OF  BULK  BREAKDOWN 
STARTING  TIMES  FOR  SAPPHIRE  USING  YAG:Nd 
LASER  IRRADIATION 
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Includes  times  after  the  peak  intensity  In  contrast  to  the 
very  sharply  defined  breakdown  starting  time  expected  for  a 
threshold-like  process.  The  moat  probable  time  for  damage  is 
in  all  cases  before  or  at  the  instant  of  maximum  intensity  or 
optical  field,  and  as  the  applied  field  is  reduced  (Pj  is 
lowered),  the  time  of  maximum  damage  probability  shifts  to 
the  peak  of  the  pulse.  (Compare  Figs,  24e  and  d.)  These 
qualitative  properties  were  shown  in  Re*1.  12  to  be  explained 
by  the  probabilistic  point  of  view. 

Quantitative  comparison  can  be  made  between  the  present 
results  and  the  earlier  studies  by  using  ar  expression  for 
the  distribution  of  breakdown  starting  timss  g(t)  which  was 
derived  by  Bass  and  Barrett.  It  was  shown  in  Ref.  12  that 
g(t)  is  given  by  the  compound  probability  that  breakdown  occurs 
at  a  particular  instant  given  that  it  has  not  occurred  before 
that  time  or 

At 

g(t)  =  h(t)  exp(-  j  h(t')  dt ’ )  (6-2' 

-CO  ' 

where  h(t)  is  the  probability  per  unit  time  that  the  applied 
laser  field  U(t)  causes  breakdown.  The  probability  p1  in 
Eq.  (3-1)  is  related  to  g(t)  by 

<20 

Pi  =  J  g(t)  dt  (6-3) 

-00 

and  the  histograms  in  Figs,  24-26  are  experimental  measures 
of  g(t).  If  the  damage  rate  h(t)  can  be  expressed  quanti¬ 
tatively,  then  g(t)  can  be  computed  and  compared  to  our 
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expert  mental  data. 

Eq.  (6-1)  suggests  that  the  damage  rate  should  be  of 
the  form12 

h(t)  =  Aq  exp(-K/E(t)) 

(6-4) 

Aq  5. 8  a  normalization  constant,  and  K  is  the  experlmen tally 

determined  quantity  in  Eq.  (6-l).  E(t)  is  determined  from 

the  peak  laser  field  E  and  the  electric  field  waveform. 

o 

If,  as  Bass  and  Barrett  had  assumed,  a  damaging  avalanche 
could  form  for  all  values  of  E(t),  then  Eq.  (6-4)  would  be 
an  appropriate  choice  for  h(t).  It  was  demonstrated  in 
Sect.  C  of  Chapt.  5,  however,  that  the  avalanche  is  time 
dependent,  when  E(t)  is  less  than  some  minimum  field  Em^n, 
the  time  required  for  the  avalanche  to  reach  damaging  pro¬ 
portions  will  exceed  the  pulse  duration  and  the  damage  rate 

must  vanish.  In  a  oulsed  field  which  reaches  E_, _  at  t  . 

mm  mm 

and  falls  below  Em|n  at  t  ,  h(t)  is  therefore  effectively 
zero  when  t  <  tm^n  and  t  >  tmax.  Assuming  that  the  avalanche 
formation  time  is  negligibly  small  (<  0,1  ns)  when  E(t)  >  Em^n, 
then  Eq.  (6-4)  should  be  replaced  by 

0  ;  t  <  tmin 

Aq  exp(-K/E(t))  ,  S-nin  *  1  *  ‘max  (6"5) 

0  »  *  >  tmax 

An  estimate  for  Em^n  can  be  made  from  the  measured  form  of 
g(t)  by  setting  equal  to  the  lowest  electric  field  value 

at  which  breakdown  has  been  observed  to  occur  in  a  particular 
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material.  The  times  t  and  t  „„  can  then  be  determined 

mm  max 

from  the  pulse  waveform. 

In  order  to  complete  the  numerical  comparison,  values 
are  needed  for  the  normalization  constant  Aq  and  the  slope  X. 
We  have  found  K  experimentally  in  fused  quartz  by  measuring 
p^  at  three  different  field  values.  The  value  so  obtained — 
71  x  108  v/cm — compares  favorably  to  the  earlier  measurement 
of  6l  x  108  v/cm  which  Bass  and  Barrett  obtained  from  their 
surface  studies  in  fused  quartz. 

The  computation  of  Aq  is  simplified  by  introducing  the 
function  v(t),  defined  as  uhe  .fractional  number  of  tests 
where  no  damage  has  occurred  by  the  time  t.  When  the  laser 
pulse  is  passed,  v(t)  must  equal  1  -  Pj »  the  probability  that 
no  damage  was  produced.  In  Ref.  12,  v(t)  was  shown  to  be 


i* 

v(t)  =  exp(-  J  h(t')  dt ' ) 


(6-6) 


Bass7"  has  calculated  AQ  from  Eqs.  (6-5)  and  (6-6)  and 
the  final  condition  on  v(t)  and  has  found  that 


1 1 1CLJV 

A0  =  -  In  d-pj)/(  Jt  exP(-K/E(t))  dt 


(6-7) 


Aq  can  thus  be  evaluated  from  experimental  parameters.  If 
the  damage  rate  assumed  in  Eq.  (6-5)  is  consistant  with  the 
experimental  data,  then  AQ  should  be  a  constant  for  a  parti¬ 
cular  material. 

Eqs.  (6-3),  (6-5),  and  (6-7)  can  now  be  used  to  calculate 
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tri*  predicted  distribut ion  of  breakdown  starting  times  for 
bulit  dv'tcufi  In  fused  quarts  which  car.  be  compared  to  the 
experimental  results  of  Fie.  ?*.  !*sing  r.u^eri cal  techniques, 
rns.93  obtained  the  curves  in  Fie.  27.  The  dar.A  used  to  com¬ 
pute  these  curve?  are  summarized  in  Table  VII. 

'ft*.  2"  shows  good  qualitative  agreement  between  the 
model  baa  2d  an  Baer  sr.i  Barrett's  ex  peri mental  respite  as 
sururi  i  jo  by  Eq  •  f  t  - 1  ;  and  our  own  experimental  data  it.  fused 
q.artz.  Th“re  is,  however,  some  deviation  between  the  two 
for  times  after  tri.-*  ic.-.k  of  the  las'5:  pulse,  that  ir,  for 
times  when  E(t)  decreases  to  £m^n.  The  form  used  for  h(t) 
when  t_.  <  t  <  t_  assumed  that  the  avalanche  formation 

time  was  negligibly  small  whenever  the  field  was  greater 
than  Ec.  As  E(t)  decreases,  however,  the  formation  time  must 
increase  at  an  increasing  rate.  This  means  that  at  such 
times  h(t)  will  be  decreased  over  the  value  in  Eq.  (6-5)  at 
a  rate  which  increases  as  t— x,  and  the  computed  values 
of  g(t)  will  be  reduced  accordingly.  The  effect  of  allowing 
for  finite  formation  time,  then,  would  make  the  computed 
curves  in  Fig.  27  decrease  more  rapidly  for  times  past  the 
peak  of  the  pulse  and  thereby  improve  the  qualitative  agree¬ 
ment  in  that  figure. 

The  agreement  between  model  and  experiment  which  has 
already  been  acnieved  is  strengthened  by  the  results  in 
Table  VII.  There  we  find  that  both  cases  were  computed  with 
nearly  the  same  value  of  AQ.  The  +  20  percent  experimental 
error  in  K  and  EQ  can  easily  account  for  the  10  percent 
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The  parameters  relevant  to  these  figures  are  listed  in  Table  VII. 

Fig.  27.  MEASURED  AND  CALCULATED  DISTRIBUTIONS  OF  BREAK¬ 
DOWN  STARTING  TIMES  FOR  BULK  DAMAGE  IN  FUSED 
QUARTZ  USING  A  YAG :  Nd  LASER 
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TABLE  VII 


Parameters  Used  to  Calculate  g(t)  in  Fig.  27 
Materials  Studied:  Fused  Quartz 


Experimentally 


Pi 

(a) 

E 

o 

K 

A 

o 

t 

max 

Jg(t)dt 

lmin 

Fig.  27a 

0.  380 

7.  3  X  108  V/m 

71  X  )08  V/m 

2599.8 

0.373 

Fig.  27b 

0.032 

5.  7  X  108  V/m 

71  X  108  V/m 

2776.4 

0.031 

(a)  E  is  the  peak  value  of  the  applied  optical  field  within  the  medium. 
Tne  electric  field  waveform  is  the  square  root  of  the  waveform  in 
Fig.  24a. 
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difference  between  the  A0’s.  As  a  further  check  on  the  con¬ 
sistency  of  the  model,  we  also  computed  the  time  integral 
of  g(t)  since,  by  definition,  it  should  equal  the  experi¬ 
mentally  determined  value  of  p1 .  in  Table  VII  we  find  that 
this  condition  is  also  satisfied.  These  computations  do  not 
involve  any  fitted  parameters  and  as  such  are  a  useful  test 
of  the  consistency  of  the  model  with  the  experimental  data. 

Finally,  it  should  be  noted  that  the  breakdown  statistics 
should  depend  on  the  size  of  the  focal  volume.  This  con¬ 
clusion,  discussed  in  Ref.  13,  follows  from  a  general  model 
which  assumes  that  the  statistics  arise  because  of  fluctuations 
in  the  formation  of  the  first  few  hot  electrons  in  the  small 
focal  volume.  As  larger  volumes  are  probed,  the  number  of 
potential  starting  electrons  increases  and  breakdown  statis¬ 
tics  should  become  less  important.  Comparable  statistics 
can  only  be  obtained,  therefore,  from  experiments  in  which 
approximately  the  same  volumes  of  material  are  irradiated. 

In  conclusion,  the  present  work  confirms  the  existence 
of  a  statistical  nature  to  the  laser  induced  intrinsic  damage 
process  and  supports  the  notion  that  the  damage  mechanism 
both  on  the  surface  and  in  the  bulk  is  an  electron  avalanche 
with  statistical  starting  properties. 

C,  Sources  of  Damage  Statistics 

It  was  shown  in  the  last  section  that  the  optical  break¬ 
down  process  has  a  statistical  character  which  leads  to  an 
experimentally  measurable  probability  of  the  form  of  Eq.  (6-1). 
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Basa  and  Barrett  explained  this  behavior  by  using  a  simpli¬ 
fied  model  for  avalanche  breakdown.  In  their  model  it  was 
assumed  that  statistics  for  the  entire  process  were  determined 
by  the  heating  of  the  first  electron  in  the  avalanche.  This 
assumption  was  based  on  the  expectation  that  the  secondary 
electrons  produced  by  an  ionizing  collision  have  energies  much 
greater  than  the  thermal  energy  kT  and  are  therefore  more 
easily  accelerated  to  ionizing  energies  than  is  the  starting 
electron.  Using  a  classical  collision  model,  similar  in 

Q  A 

conceptual  form  to  a  dc  model  by  Shockley,  they  were  able 
to  derive  Eq.  (6-1)  for  a  laser  pulse  of  constant  intensity. 

The  model  of  Bass  and  Barrett  is  greatly  oversimplified 
in  computational  details  because  of  the  use  of  a  classical 
relaxation  time  approximation  with  an  energy-independent 
collision  rate.  In  addition,  the  most  pronounced  damage 
statistics  in  their  study  are  observed  with  the  pathological 
crystal  SrTiO^  which  is  known  to  possess  large  surface  strains 
and  enormous  local  field  effects.®1  Yet  the  idea  that  the 
heating  of  starting  electrons  determines  the  observed  damage 
statistics  i 8  an  intriguing  one,  and  their  model  is  consis¬ 
tent  with  experimental  results.  We  will  reconsider  the  matter 
of  starting  electrons  after  looking  for  other  sources  of 
damage  statistics  in  addition  to  the  dynamics  of  the  electron 
avalanche. 

It  was  observed  by  tests  such  as  those  summarized  in 
Sect.  D  of  Chapt,  3  that  the  damage  statistics  do  not  result 
from  inclusions  and  fixed  material  inhomogeneities.  A 
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measurement  has  been  made  of  the  density  of  inclusions  on 
the  surface  of  fused  quart2  to  show  that  this  conclusion  is 
reasonable.  At  our  request  Dr.  0.  J.  Guentert  of  Raytheon 
Corp.  conducted  an  electron  microprobe  survey  of  the  surface 
of  a  piece  of  conventionally  polished  fused  quartz  used  in 
the  actual  damage  measurements  of  Ref.  12.  He  found  that 
there  were  no  more  than  about  100  impurities  per  cm.  The 
microprobe  could  sense  impurities  with  diameters  >  1  gm  and 
to  a  deDth  of  5  pm.  In  the  surface  damage  measurements  of 
the  last  section  and  those  of  Bass  and  Barrett,  the  high 
intensity  (l/e)  region  of  the  focused  laser  beam  was  about 
10  tim  in  diameter,  and  damage  was  observed  to  commence  within 
a  few  gm  of  the  surface.  On  the  average,  then,  the  pro¬ 
bability  of  intercepting  a  lar?e  inclusion  was  about  10“^  per 
shot.  If  we  assume  that  the  density  of  potentially  damaging 
inclusions  with  diameters  d  such  that  0.1  gm  <  d  <  1  gm  is 
about  100  times  as  great  as  Guentert  had  measured  for  large 
inclusions,  then  the  probability  increases  to  10  ,  a  value 

which  is  still  much  less  than  unity.  We  thus  see  that  the 
experimental  confirmation  of  the  absence  of  inclusion  damage 
on  the  surface  is  consistant  with  the  measured  density  of 
inclusions.  A  similar  result  holds  for  bulk  damage.  It 
should  be  noted,  incidentally,  that  the  use  of  small  beam 
diameters  resulting  from  tightly  focusing  external  optics  is 
necessary  in  order  to  avoid  inclusion  damage.  Many  studies 
of  optical  damage  have  been  conducted  with  ‘laser  beam  diameters 
<  100  ^5*  ^9  ancjf  therefore,  probably  reflect 
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characteristics  of  inclusion  da.., age  rather  than  characteristics 
of  Intrinsic  processes. 

It  has  been  suggested  that  dust  particles  settling  on 

op 

the  surface  could  lead  to  surface  damage  statistics.  While 

dust  may  affect  the  results  of  surface  studies,  it  is  difficult 

to  see  how  statistics  from  the  settling  of  dust  can  have  the 

form  of  Eq.  (6-l).  Also,  statistics  have  now  been  observed 

in  the  bulk  where  surface  dust  does  not  present  a  problem. 

Another  possible  source  of  statistics  is  from  thermally- 

induced  migration  of  absorbing  impurities.®-*  This  orocess  can 

lead  to  statistical  distributions  at  fixed  laser  power  of  the 

form  of  Fig.  15  in  Chapt.  3«  Let  us  attempt  to  model  such  a 

process.  It  might  be  expected  that  some  activation  energy 

< „  exists  for  the  migration.  The  driving  force  is  a  heat 
g 

pulse  injected  by  a  small  residual  absorption  of  the  laser 

light,  If  large  volumes  are  considered,  the  temoerature  rise 

5T  -  I  tp,  such  as  we  found  for  inclusions  in  Chapt,  1,  and 

if  small  volumes  are  considered,  ST  -  the  laser  intensity  I 

2 

or,  equivalently,  ST  ~  E  .  The  probability  that  such  a 
migration  will  occur  by  laser  irradiation  is  proportional  to 
a  Boltzman-like  factor  exp  (-<g/k  T) .  At  first  glance  this 
factor  appears  to  be  of  the  form  of  Eq.  (6-1)  but  with  closer 
study  we  realize  that  the  present  exponential  function  depends 
on  E2  and  not  E  as  observed  in  Fig.  21  and  as  summarized  by 
Eq.  (6-1).  This  is  a  general  resulti  Physical  observables 
from  a  thermal  process  should  be  functions  of  E  and  not  E. 

The  damage  Drobability  p^  is  a  function  of  E,  however,  so  that 
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we  can  oonclude  that  no  simple  thermal  process  is  causing 
the  damage  statistics. 

There  are  other  possible  sources  of  damage  statistics. 

In  dc  experiments  it  is  very  likely  that  damage  statistics 
result  from  space  charge  development  which  intensifies  the 
local  field.  Rut  space  charges  should  not  develop  in 
optical  fields  or,  for  that  matter,  in  quasi-dc  fields  of 
nanosecond  duration  so  that  this  model  cannot  explain  our 
results.  A  cumulative  effect  from  the  trapping  of  electrons 
in  low  lying  states  may  in  principle  influence  the  experi¬ 
mental  results, 55  put  it  was  found  that  the  damage  statistics 
were  not  a  function  of  repetition  rate  as  that  rate  was 
varied  from  about  1  pps  to  h-  pps.  Also,  the  life-time  of  such 
trapped  states  is  less  than  a  second. 55 

Statistics  in  laser-induced  gas  breakdown  can  occur 

because  there  may  be  no  free  electrons  in  the  irradiated 
QU 

volume.  Let  us  consider  this  mechanism  as  a  source  for 
damage  statistics  in  solids.  In  an  insulator  the  density  of 
thermally  excited  conduction  electrons  in  zero  field  is  of 
the  order  of  1012/cm5  or  more.  The  high-intensity  region  of 
the  focused  laser  beam  where  bulk  damage  commences  (see  Fig.  13 
in  Chapt,  3)  is  about  10”^  cm^  in  volume  so  that  about  10^ 
quasi-free  electrons  should  exist  in  that  volume.  This 
number  is  adequate  to  insure  that  at  least  one  ionized  electron 
is  present  in  the  focal  volume.  Multiphoton  ionization  of 
impurities  should,  in  fact,  increase  this  number.  It  thus 
appears  that  an  initially  ionized  (conduction-band)  electron 
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wlll  always  be  present  In  the  focal  volume. 

Bass  and  Barrett  have  recently  made  an  observation  which 
gives  further  support  to  the  conclusion  that  statistics  are 
present  in  the  avalanche  itself.1^  When  LHO3  and  SrTU^ 
were  damaged  with  ruby  laser  light,  no  intrinsic  statistics 
could  be  seen.  This  was  not  the  case,  however,  with  three 
other  materials  studied  at  0.69  gm.  The  interpretation  given 
for  the  abrupt  change  in  damage  characteristics  was  that 
another  intrinsic  mechanism,  multiphoton  absorption,  was 
responsible  for  damage  in  LiI03  and  SrTi03  by  ruby  laser  light. 
Damage  from  multiphoton  absorption  is  not  expected  to  have  a 
measurable  statistical  character  because  multiphoton  ioni¬ 
zation  does  not  require  a  sequence  of  improbable  events.  This 
is  not  the  case  with  avalanche  ionization,  however,  because 
the  development  of  the  avalanche  requires  some  electrons  to 
be  heated  to  energies  well  in  excess  of  the  average  electron 
energy.  (See  Chapt.  8.)  In  support  of  their  interpretation, 
Bass  and  Barrett  noted  that  two-photon  absorption  is  possible 
in  LiNbO^  and  SrTiC^  at  0.69  /im  but  it  is  not  possible  at 
1.06  ix m.  The  other  materials  studied  at  both  wavelengths  have 
band-gap  energies  in  excess  of  3.6  ev  so  that  two-photon 
absorption  is  not  possible  in  them  at  either  1.06 /xm  or 
0.69  MW*  If  this  interpretation  is  correct,  these  measure¬ 
ments  represent  the  first  experimental  confirmation  of  the 
competition  between  multiphoton  absorption  and  avalanche 
breakdown  in  solids,  (See  the  discussion  at  the  end  of 
Chapt.  1.) 
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The  source  of  the  damage  statistics  thus  appears  to  be 
in  the  electron  avalanche  itself-- through  the  details  of  the 
initial  electron  heating,  because  of  the  chance  termination 
of  incipient  avalanches,  or  because  of  fluctuations  in  the 
ionization  rate  a.  Yablonovitch55  has  shown  that  the  latter 
process  cannot  lead  to  measurable  statistics*  The  second 
possibili\  is  virtually  identical  to  the  first  because  it  is 
unlikely  that  an  avalanche  which  has  produced  more  than  a  few 
generations  can  terminate  with  any  reasonable  probability. 

We  are  thus  left  with  the  basic  hypothesis  of  Bass  and 
Barrett,  that  measurable  damage  statistics  can  develop  from 
the  first  generations  in  the  avalanche. 

Bass  and  Barrett  have  suggested  a  description  of  those 
first  generations.  A  relatively  simple  improvement  of  their 
model  would  be  to  use  an  energy-dependent  collision  rate  such 
as  discussed  in  Chapt.  8,  Alternately,  it  may  be  possible 
to  show  that  the  statistics  of  the  first  events  do  not  result 
from  the  heating  of  the  first  electron  to  the  ionization 
energy  I  but,  rather,  they  result  from  the  heating  of  the 
electron  to  energies  just  greater  than  the  longitudinal  optica] 
phonon  (L0)  resonances  at  about  0,01  I.  Thornber  and 
Feynman^  have  shown  that  electron  energy  losses  to  the 
lattice  are  extremely  high  near  the  LO  energies.  These 
resonances,  therefore,  represent  barriers  to  electron  heating 
and  may  effectively  produce  low-energy  electron  traps.  Be¬ 
cause  there  are  probably  only  of  the  order  of  10^  quasl-free 
electrons  in  the  focal  volume  at  the  start  of  the  laser 
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pulse,  there  may  be  no  electrons  initially  present  with 
energies  past  the  LO  resonances.  Once  an  electron  moves  past 
this  barrier  because  of  a  sequence  of  favorable  but  improbable 
collisions  with  the  lattice,  the  heating  may  be  rapid.  Pro¬ 
ducts  of  the  ionization  will  probably  have  energies  which  are 
much  greater  than  the  LO  energy  so  that  the  avalanche,  once 
begun,  will  not  terminate. 

We  thus  see  that  there  is  ample  evidence  that  optical 
avalanche  breakdown  has  a  measurable  statistical  character. 

From  a  practical  viewpoint,  however,  it  is  not  clear  that 
damage  statistics  such  as  we  have  measured  are  important  to 
laser  system  design.  It  is  the  author's  experience  that  when 
inclusion  damage  develops,  variations  in  optical  strength  are 
normally  much  more  pronounced  than  the  statistics  we  have 
observed  from  intrinsic  damage.  In  fact,  the  observation  of 
a  reasonably  well-defined  damage  threshold  is  itself  an  indi¬ 
cation  that  intrinsic  damage  is  occurring.  (See  Sect.  D  of 
Chapt.  3*)  Intrinsic  damage  statistics  have  been  observed 
because  inclusions  have  been  largely  avoided  in  our  measure¬ 
ments  and  in  those  of  Ref,  12,  stable  single-mode  laser  sources 
have  been  used,  great  care  has  been  taken  in  recording  damage 
levels  and  only  small  focal  volumes  have  been  irradiated. 

From  a  fundamental  viewpoint,  on  the  other  hand,  the  study 
of  avalanche  statistic3  may  give  us  valuable  insight  into  the 
details  of  the  avalanche  formation,  insight  which  could  not 
be  gained  by  observation  of  other  damage  parameters. 


CHAPTER  7 


OPTICAL  SURFACE  DAMAGE 

A.  Introduction 

12  13 

In  the  experiments  of  Bass  and  Barrett,  *  the  laser  beams 

2 

were  focused  to  about  25  pm  in  diameter  {l/e  in  intensity)  in  order  to 

induce  damage  on  the  surface.  For  such  focal  diameters  the  laser 

beam  was  effectively  collimated  over  a  distance  of  about  50  pm  near 

86 

the  focus,  and  yet  damage  was  consistently  observed  to  develop 

within  the  firs'.  0.  25  pm  of  material.  Because  surface  contaminants 

were  apparently  not  responsible  for  damage  in  this  study  (see  Chapt.  6), 

this  surprising  result  indicated  that  the  surfaces  of  conventionally 

polished  solids  have  a  lower  intrinsic  threshold  than  the  bulk.  Giuliano 

87 

observed  in  a  later  work  with  a  weakly  focused  ruby  laser  beam  that 

the  surface  damage  field  for  sapphire  was  increased  when  the  surface 

was  ion-beam  polished,  a  procedure  that  removed  most  but  not  all  of 

the  polishing  scratches  on  his  samples. 

A  quantitative  relationship  between  surface  damage  fields  and 

those  of  the  bu'k  was  not  established  experimentally  in  these  or  in  any 

other  published  study,  because,  to  the  author's  knowledge,  until 

14 

the  work  of  Yablonovitch  and  the  present  work,  no  accurate  measure¬ 
ments  of  intrinsic  optical  bulk  damage  had  every  been  made.  Using 
the  bulk  damage  techniques  of  C.hapts.  4-6  and  the  surface  focusing 
scheme  of  Ref.  12,  we  have  completed  a  study  of  surface  damage 
at  1.  06  pm  in  which  surface  damage  fields  were  measured  for  a 
variety  of  surface  preparations  and  compared  to  bulk  damage  fields. 
A  simple  model  proposed  by  Bloembergen^  is  found  to  be  fully 

Preceding  page  blank 


-154- 


caoable  of  explaining  the  experimental  results. 

B .  The  Experiment 

The  laser  source  and  experimental  procedures  for  studying 
intrinsic  bulk  damage  were  discussed  in  detail  in  Chapts.  3 
and  4.  A  simple  scheme  for  accurately  focusing  on  the  surfaces 
of  solids  was  developed  in  Ref.  12  and  applied  in  the  present 
work.  In  this  scheme  the  focusing  lens  for  the  laser  beam  is 
used  as  a  microscope  objective.  Laser  light  enters  the  micro¬ 
scope  between  the  eyepiece  and  objective  and  is  reflected  off 
a  beam-soittter  towards  the  objective  lens.  By  inducing  a 
spark  with  laser  light  on  the  surface  of  a  piece  of  steel  shim 
stock,  it  is  possible  to  accurately  place  the  shim  stock  at 
the  focus  of  the  objective  lens.  When  the  separation  of  the 
two  lenses  in  the  eyepiece  is  adjusted  to  bring  the  surface  of 
the  shim  stock  into  a  sharp  visible  focus,  a  calibration  is 
established  such  that  the  focused  object  plane  for  the  observer 
who  uses  the  total  eyepiece-objective  system  is  the  same  as 
the  focal  plane  for  the  laser  light  which  enters  only  the 
objective  lens.  In  this  manner  the  separation  of  the  objective 
lens  and  the  sample  was  adjusted  by  simple  visual  inspection 
to  accurately  focus  the  laser  light  on  the  sample  surface. 

For  each  sample  studied,  the  surface  damage  field  Es  was 
measured,  and  by  bringing  the  sample  about  1  mm  closer  to  the 
objective  lens,  the  bulk  damage  field  Eg  was  determined  for 
that  same  sample.  All  field  measurements  were  made  relative 
to  the  bulk  damage  field  of  NaCl  which  was  known  from  the 
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experiment  of  Chapt.  4  to  an  accuracy  of  ?0<.  ’‘el.vtive 
measurements  of  damage  fields  wero  made  +  ■.;  c»n  accuracy  of 
better  than  1(X.  The  damare  field  itself  was  defined  in  the 
usual  manner  as  that  value  of  rms  electric  field  which  pro¬ 
duced  damage  half  of  tho  time  on  R  single  shot  ( l^e . ,  t  5  11  c,. ;-} , 

Table  VIU  summaries  the  measured  ratios  Kg/%  for  three 
materials  whose  surfaces  have  hern  pr e t« red  tv  various  rmns. 
These  data  show  that  the  clean  but  conventionally  polished 
surface  (Fig.  28a)  of  a  transparent  medium  Is  gcnernlly  more 
easily  damaged  than  the  hulk.  Or.  the  other  hand,  when  cave  fu 
taken  to  achieve  imperfection.- free  surface  finishes  such  us 
shown  in  Fig,  28b,  the  bulk  and  si:rfcc*»  damage  fields  are 
equal. 

The  momhology  of  tie  reel  dual  entrance  surface  damage 

gives  additional  support  for  the  preceding  conclusions.  When 

a  spark  was  observed  during  laser  irradiation,  very  small  pit 
17 

damage  was  found  on  a  convent i onally  finished  surface.  Im¬ 
perfection-free  surfaces  showed  mure  extensive  damage  in¬ 
cluding  a  region  of  clacks  snrroundin  the  melted  irradiated 
volume  (Fig.  29).  Such  damage  extends  many  microns  into  the 
material  and  is  very  similar  to  a  cross-section  of  bulk  damage 
as  recorded  in  Fig.  1.3  of  Chapt.  3. 

Since  the  focused  laser  light  was  effectively  collimated 
over  several  microns,  near  the  surface,  high  intensity  light 
extended  into  the  bulk,  and  it  was  Dossible  for  the  light  to 
cause  internal  or  bulk  damage.  A  conventionally  finished 
surface,  however,  damages  more  easily  than  the  bulk.  Because 


Very  faint  pit  ~  20  pm  in  dia.  and  ~  0.25  pm  deep 

Extensive  cracking  out  to  -  150  pm  with  central  hole  -  20  pm  dia.  and  ~  3-10  pm  deep 

TABLE  VIII 

Comparison  of  Bulk  and  Surface  Damage  Fields  for  Different  Samples  and  Surface  Finishes 
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Fig.  29.  RESIDUAL  SURFACE  DAMAGE  ON  ION-BEAM 


POLISHED  FUSED  QUARTZ 


Surface  which  were  effectively  imperfection-free  produced  the  form  of 
damage  morphology  pictured  here.  Cracks  surround  a  small  molten  region 

37  .  .  7 

about  10  pm  in  diameter.  Bass  has  observed  a  similar  morphology  when 
damage  occurs  from  self-focusing  inside  the  bulk  and  breaks  out  through 
the  surface. 
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the  onset  of  damage  is  accompanied  by  a  strong  attenuation  of 

the  transmitted  beam,  surface  damage  on  conventionally  finished 

samples  is  localized  on  the  surface.  If,  on  the  other  hand,  a 

more  highly  finished  surface  does  not  damage  any  more  easily 

than  the  volume  of  the  medium,  it  is  possible  for  intense  light 

to  reach  the  bulk.  Damage  can  then  commence  anywhere  within 

the  effective  focal  volume.  Occasionally  damage  will  begin 

within  0.25  fi  m  of  the  surface  and  produce  a  faint  pit,  but 

more  often  it  will  begin  somewhere  in  the  bulk  several  microns 

from  the  surface.  When  bulk  damage  occurs,  the  molten  region 

grows  back  towards  the  laser  and  breaks  out  of  the  entrance 

surface  producing  the  morphology  recorded  in  Fig.  29.  These 

morphological  observations,  which  are  valid  only  when  inclusion 

induced  damage  is  absent,  show  that  if  very  small  pit  damage 

is  normally  produced,  then  that  surface  is  more  easily  damaged 

than  the  bulk.  The  observations  that  damage  on  imperfection- 

free  surfaces  appears  to  be  bulk  damage  which  erupted  from 

within  the  medium  might  explain  Giuliano's  report  of  more 

extensive  surface  damage  on  ion-beam  polished  sapphire  than 

87 

on  a  conventionally  polished  sample. 

The  morphology  of  damage  on  the  bowl  feed  ( " superoolished" ) 
surface  of  BSC-2  glass  was  generally  of  the  small  pit  type 
though  more  extensive  damage  was  detected  occasionally.  The 
measured  ratio  of  bulk  to  surface  damage  fields,  1.0  ±  0.1, 
implied  a  perfect  surface,  but  based  on  the  damage  morphology 
we  decided  that  the  surface  was,  in  fact,  not  as  uniformly 
imperfection-free  as  oossible.  Not  long  after  this  conclusion 
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was  reached,  we  received  the  electron  micrograph  shown  in 
Fic.  28c  which  confirmed  our  interpretation  of  the  surface 
quality.  Faint  residual  scratches  can  still  be  seen  in  that 
figure. 

Bloembergen  has  completed  a  simple  but  intriguing  model* ^ 
which  predicts  that  the  apparent  optical  strength  of  surfaces 
will  normally  be  less  than  the  optical  strength  of  the  bulk. 

He  recognized  that  the  dielectric  discontinuity  occurring  at 
structural  defects  can  enhance  the  electric  field  over  its 
average  value.  The  following  values  are  predicted  for  the 
ratios  Ep/Ec;t 

Spherical  void  3</(2f  +  l) 

Cylindrical  groove  2*/(c  +  1) 

"Vee"  groove  e 

where  t  =  optical  frequency  dielectric  constant.  These  ex¬ 
pressions  were  used  to  obtain  the  predicted  ratios  in  Table 
VIII. 

It  is  seen  that  except  for  sapphire  the  measured  ratios 
for  conventionally  oolished  samples  lie  within  the  range  of 
values  which  Bloembergen  predicts  for  spherical  voids  and 
cylindrical  grooves.  Actual  defects  on  conventionally 
polished  surfaces  probably  have  more  complicated  shapes  than 
Bloembergen  has  considered.  But  a  high  density  of  subsurface 
voids  and  surface  pits  does  result  from  standard  abrasive 
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polishing8^  as  seen  in  FI*.  28a,  so  that  the  reasonable 
agreement  between  the  first  two  predicted  ratios  and  experi¬ 
ment  is  impressive.  The  sapphire  sample  had  visible  scratches 
that,  as  Bloembergen  predicted,  result  in  a  large  damage 
ratio.  Finally,  the  bowl-feed  and  ion-beam  polished  samples 
appear  to  have  only  gentle  variations  in  surface  topography 
(Fig.  28b)  which  should  lead  to  little  if  any  measurable 
field  enhancement  according  to  Eloembergen' s  model.  Again 
there  is  experimental  agreement.  It  should  be  noted,  inci- 
dently,  that  good  quality  conventional  polishes  (Fig.  28a) 
appear  spatially  uniform  to  a  20  ptm-diametor  laser  beam.  In¬ 
trinsic  damege  on  such  a  surface  will  produce  damage  statis¬ 
tics  of  the  form  of  Fig.  1 5  in  Chapt.  3. 

Another  ooint  of  interest  in  Table  VIII  is  the  fact  that 
the  first  and  third  fused  quartz  samples  were  purchased  from 
the  same  manufacturer  at  the  same  time.  These  were  described 
as  Hl/20  wave  interferometer  flats, "  Ion  beam  polishing  of 
the  fourth  samole  removed  material  to  a  depth  greater  than 
3  times  the  smallest  grit  size  used  in  the  mechanical  polishing 
process,  and  a  large  increase  in  the  surface  damage  field 
resulted.  The  ion  beam  polisher  employed  Ar+  ions  at  10  kV 
and  incident  at  60°  to  the  surface  normal  to  remove  1.25  gm 
of  material. 

Table  IX  gives  the  measured  bulk  damage  fields  for  the 
samdes  in  Table  VIII.  These  values  were  measured  relative 
to  the  damage  field  for  NaCl,  and  most  of  the  assumed  errors 
in  Table  IX  are  due  to  the  experimental  uncertainty  in  that 
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TABLE  IX 

Measured  Bulk  Damage  Fields 

at  1.06  pm 

Material 

Eg(MV/cm)(a) 

Fused 

Quartz 

Conventional  #1 

Conventional  #2 

Bowl  feed  finish 

Ion  beam  polished 

5.2  ±  1.3 

Sapphire 

6. 3  ±  1.6 

BSC-2  Gl&ss 

4.7  ±  1.2 

(a)  RMS  field  values  are  listed 

Measurements  were  made  relative  to  Eg  in  NaCl 


latter  value. 

In  conclusion,  we  see  that  the  incident  laser  fields  for 
entrance  surface  and  bulk  damage  to  transparent  solids  are  the 
same  fo-  imperfection-free  surfaces.  When  structural  imper¬ 
fections  are  present,  they  enhance  the  optical  field  locally 
so  that  the  surface  damages  more  easily  than  the  bulk.  Im¬ 
proved  surface  finish  quality  is  recognized  as  essential  to 
minimizing  the  problem  of  scattering  in  optical  systems.  The 
results  of  this  work  show  that  it  is  also  critical  in  the 
problem  of  laser  induced  surface  damage. 


CHAPTER  8 


THEORIES  OF  DC  AND  OPTICAL  ELECTRIC  BREAKDOWN 

A.  Introduction 

Electric  breakdown  is  one  of  the  most  difficult  problems  of  modern 
solid  state  theory.  As  is  probably  the  case  with  most  hot  electron  phen¬ 
omena  in  solids,  the  process  is  understood  only  in  qualitative  terms. 

The  basis  of  the  difficulty  in  modeling  breakdown  is  the  fact  that  the 
relevant  interactions  cannot  be  treated  separately.  It  is  not  strictly 
correct,  for  example,  to  speak  of  electron  excitations  as  being  distinct 
from  those  of  the  lattice.  The  electron-phonon  interaction  in  NaCl,  for 
example,  is  so  strong  that  the  effective  collision  time  appears  to  be 
~2  X  10  ^  sec.  (Sect.  B  of  Chapt.  6)  which  implies  an  electron  energy 
uncertainty  of  the  order  of  the  ionization  energy.  In  addition,  distortion 
from  electric  fields  having  rms  values  ~10b  volts /cm  cannot  properly 
be  treated  as  a  perturbation  on  the  energy  spectrum  of  the  electron- 
lattice  system. 

A  number  of  theories  have  been  formulated  which  largely  ignore 
these  difficulties.  Many  such  theories  use  a  relaxation  time  approximation 
with  the  electron-phonon  collision  time  calculated  from  perturbation 
theory.  Despite  their  lack  of  rigor,  such  models  appear  to  explain  many 
features  of  the  process. 

In  this  chapter  we  will  briefly  review  some  of  the  published  models 

of  dc  and  optical  electric  breakdown  and  compare  their  predictions  with 

6  7  68 

the  experimental  results  of  Chapts.  4  and  5,  A  recent  work  by  Holway 
will  be  considered  in  which  the  high-field  electron  distribution  function  is 


Preceding  page  blank 


166 


calculated  from  a  classical  Fokker-Planck  equation.  Although 
Holway’s  model  greatly  simplifies  the  optical  breakdown 
process,  it  demonstrates  that  electric  breakdown  may  not  be 
an  average  electron  phenomenon.  The  model,  in  addition,  can 
provide  a  useful  framework  for  evaluating  many  untested 
theoretical  ideas  concerning  electric  breakdown. 

B.  Historical  Review  of  Breakdown  Theories 

The  literature  on  breakdown  distinguishes  between  two 
general  types  of  dielectric  failure  called  intrinsic  breakdown 

£  O 

and  avalanche  breakdown.-  In  intrinsic  breakdown  theories, 
the  quasi-free  electron  density  is  seen  to  jump  discontinuously 
to  damaging  Droportions  from  a  low  steady-state  level  as  the 
external  field  is  raised  by  a  small  amount.  Calculated 
thresholds  are  independent  of  sample  thickness  (dc)  or  duration 
of  the  apolied  field.  Avalanche  theories,  on  the  other  hand, 
consider  a  time-dependent  electron  multiplication  whose 
threshold  is  sensitive  to  the  time  available  for  the  avalanche. 
As  the  sample  thickness  is  decreased  or  as  the  time  duration 
of  the  applied  field  is  reduced,  the  electric  field  necessary 
to  induce  damage  will  increase. 

Experimental  evidence  suggests  that  the  idealized  intrin¬ 
sic  mechanism  does  not  exist.  Materials  such  as  the  alkali 
halides  and  mica,  which  were  originally  believed  to  fail  by 
intrinsic  breakdown,  have  been  shown  to  have  a  time-dependent 
damage  field  at  dc  and,  as  we  saw  in  Chapt.  5,  at  optical 
frequencies.  Yet  the  conceptual  distinction  between  an 
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intrinsic  instability  and  a  time-dependent  avalanche  persists 
in  the  recent  theoretical  literature,  and  in  this  section  we 
will  group  breakdown  theories  according  to  this  distinction. 

Intrinsic  breakdown  theories  consider  a  balance  between 
electron  energy  gain  A  from  the  field  and  electron  energy  loss 
B  to  the  lattice.  Equilibrium  exists  when  the  relationship 


A  E  ,  g)  =  B(g! 


(8-1) 


can  be  satisfied  in  a  physically  acceptable  manner.  0  E  is 
the  external  electric  field  and  /3  symbolizes  various  material 
parameters  that  determine  A  and  B.  When  the  field  is  raised 
to  a  critical  value  Ec,  Eq.  (8-1)  can  no  longer  be  satisfied 
in  a  physically  acceptable  manner  and  breakdown  is  assumed  to 
develop.  Ec  is  identified  as  the  breakdown  field.  Various 
theories  of  intrinsic  breakdown  differ  in  two  basic  aspects-- 
the  modeling  of  the  manner  in  which  the  electron  population  is 
coupled  to  the  lattice  ( i .e. ,  in  the  calculation  of  B)  and  the 
explicit  assumption  of  what  constitutes  a  physically  acceptable 
condition.  We  will  consider  important  examples  of  such 
theories . 

One  of  the  most  frequently  referenced  theories  of  dc 
breakdown  is  the  original  model  of  FrBhlich.^  In  this  work 
Frbhlich  considered  ionic  insulators  and  assumed  that  the 
electron-lattice  coupling  arises  entirely  from  the  interaction 
between  electrons  and  longitudinal-optical  (LO)  phonons.  A 
perturbation  calculation  of  the  electron-phonon  interaction 
leads  to  an  effective  collision  frequency  =  l/r  which 
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decreases  with  increasing:  electron  energy  r  for  t  greater 
than  the  LO  phonon  energy  1vjLq •  In  particular,  when 
«  > 


k<G  (8-2) 

e'3/2  k>G 

where  k  is  the  free  electron  wave  vector  and  G  is  the  reciprocal 
lattice  wavevector.  For  NaCl,  n  G/2m  =  0. 5  I  where  m  is  the 
free  electron  mass  and  I  is  the  ionization  energy.  The  term 
B  (/9)  in  Ea .  (8-1),  which  describes  the  energy  transfer  to  the 
lattice,  is  proportional  to  v  times  a  factor  depending 
logarithmically  on  the  electron  energy. ^  That  is. 


A(E,/3)  has  the  form  of  Eq.  (1-1)  with  n  -  er/m.  At  dc  ( (o  =  0) 
the  energy  gain  for  each  electron  is 


(8-4) 


FrOhlich  defined  the  breakdown  field  Ec  as  that  value  of 
electric  field  for  which  A  and  B  intersect  at  the  ionization 


energy  I.  When  E  >  Ec,  the  point  of  intersection  is  at 
<  I,  and  when  E  <  Ec,  the  point  of  intersection  is  at 
<2  >  3»  This  behavior  is  shown  schematically  in  Fig.  30a. 
At  E  =  Ec,  electrons  with  energies  greater  than  I  gain  more 


energy  on  the  average  from  the  field  than  they  lose  to  the 
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lattice,  and  electrons  with  <  <  I  lose  more  energy  on  the 
average  to  the  lattice  than  they  gain  from  the  field. 

The  basis  for  defining  Ec  in  this  manner,  is  the  in¬ 
tuitive  notion  that  a  steady-state  can  only  exist  when 
collision  ionization  is  balanced  by  the  inverse  process  of 
recombination.  If  E  <  Ec,  then  a  high  energy  electron  (  *  =  I) 
created  by  recombination  will  lose  its  energy  to  the  lattice, 
but  if  E  >  Ec,  the  high  energy  electron  created  by  recombina¬ 
tion  will  continue  to  gain  energy  on  the  average  from  the 
field  until  it  produces  another  ionization.  The  precise 
manner  in  which  low  energy  electrons  reach  t  =  T  and  cause 
ionization  is  not  soecified. 

Frfthlich  has  thus  considered  a  possible  electron  insta¬ 
bility,  but  he  did  not  relate  this  instability  to  the  actual 

dynamics  of  electron  multiplication  and  eventual  material 

<4 

damage.  An  earlier  theory  by  von  Hipoel,  which  was  developed 

88 

further  by  Callen,  is  based  on  the  far  more  stringent  re¬ 
quirement  that  E  must  be  sufficiently  large  to  move  the  point 
of  intersection  of  the  curves  for  A  and  B  to  r  =  0,01  I. 

The  Frbhlich  model  predicts  breakdown  fields  about  an  order  of 
magnitude  larger  than  experimentally  observed. 

Zverev  et  al. have  considered  the  Frbhlich  model  in  the 
high  frequency  limit  by  setting  gj  »v  in  Eq.  (1-1 ).  This 
changes  A  (E,/3)  to 


zrz 

e  E 


mu; 


v  (e) 

b 


(8-5) 
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B(/3)  is  again  given  by  Eq.  (8-3).  In  the  energy  range  of 
interest,  the  collision  frequency  vt(<)  is  given  by  the  second 
(high  momentum)  form  of  Eq.  (8-2). 

Fig.  30b  shows  the  functional  form  of  A  and  B.  It  is 
seen  that,  in  this  high  frequency  limit,  when  E  =  Ec,  electrons 
with  <  >1  lose  energy  to  the  lattice  while  those  with  *  <  I 
gain  from  the  field  on  the  average--the  opposite  of  the  dc 
benavior.  The  intuitive  basis  for  defining  E,,.  therefore,  is 
not  valid  in  the  limit  Zverev  et  al.  have  considered.  In 
addition,  their  calculated  breakdown  strength  for  sapnhire 
is  about  an  order  of  magnitude  larger  than  we  have  measured 
at  1.06  gm.  (nee  Table  IX  in  Chapt.  7.) 

A  second  group  of  intrinsic  breakdown  theories  use  a 
Boltzmann  equation  to  calculate  an  electron  distribution 
function  f ( O .  We  will  consider  these  theories  only  briefly. 

As  before,  the  perturbation  results  of  Eqs.  (8-2)  and  (8-3) 
are  used  to  model  the  energy  transfers  among  the  field,  the 
electrons,  and  the  lattice.  A  general  solution  for  f(0  is 
exoanded  in  terms  of  spherical  harmonics  with  only  the  first 
two  terms  retained.  Breakdown  is  assumed  to  occur  when  either 
an  effective  electron  temperature  increases  without  limit  or 
when  the  distribution  function  can  no  longer  be  normalized. 

The  minimum  electric  field  which  produces  these  conditions  is 
defined  to  be  the  damage  field  Ec. 

A  number  of  authors,  including  Frfthlich,®^"*^  have  used 
this  formalism  to  describe  dc  breakdown.  The  simplification 
is  normally  made  that  f(<)  is  essentially  a  Boltzmann 
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distribution  even  though  this  assumption  is  valid  only  at  low 
fields  or  under  conditions  where  the  electron  density 
n  >  101-®  cnT®.^  A  high  frequency  (optical)  extension  of  this 
formalism  was  carried  out  by  Wasserman  who  assumed  that 
o»r  >>1  in  Rq .  (1-1).^  He  found  that  an  effective  electron 
temperature  could  not  be  found  in  NaCl  when  E  >  2  x  lO'7  volts/cm. 
This  field,  assumed  to  be  the  damage  field,  is  the  same 
magnitude  as  the  critical  fields  obtained  in  the  dc  theories, 
but  it  is  a  factor  of  10  larger  than  the  optical  breakdown 
field  we  measured  in  Chapts.  4  and  5* 

Avalanche  breakdown  theories  consider  damage  to  develop 
when  a  time-dependent  electron  multiplication  has  increased 
the  conduction  band  densitv  to  damaging  proportions.  As  was 
the  case  with  intrinsic  breakdown  theories,  the  modeline 
eenerally  takes  one  of  two  forms,  either  calculations  of 
average  electron  dynamics  or  calculations  of  the  high-field 
electron  distribution  function, 

ft  0 

A  phenomenological  dc  breakdown  model  due  to  Shockley 
disolavs  some  simple  features  of  an  electron  avalanche. 

According  to  Shockley,  the  ionization  rate  at  dc  (see  Rq .  (5-2) 
is  proportional  to  the  probability  P(tj)  that  a  low  energy 
electron  avoids  ohonon  collisions  for  a  sufficiently  long 
period  of  time  tj  to  be  accelerated  by  the  dc  field  to  the 
ionization  energy  1. 

We  can  find  P  by  the  following  argument.*’®  The  probability 
that  an  electron  will  avoid  collisions  for  the  period  (t  +  dt) 
is  eiven  by  the  compound  Probability  that  it  avoids  collision 


for  the  period  t  and  the  probability  that  it  does  not  collide 
with  a  nhnnon  during  the  interval  dt.  That  is 


P(t  +  dt)  =  P(t)(  1  -  v  (D  dt) 

1/ 


(8-6) 


The  solution  to  Eq.  (8-6)  is 


P(tj)  =  exp 


v  (s)  dt 

h 


(8-7) 


Time  and  energy  are  related  for  a  free  electron  in  a  dc  field  by 


_  /  dv  \  i.  ,1/ 

dt  ‘  mv(d t)  =  Umc) 


.1/2  eE 
m 


(8-8) 


which  can  be  inserted  into  Eq .  (8-7)  to  give 

P(t1>  «  exp  (-K/E)  *  i  ^  £  a  .  (8-9) 

N  is  the  density  of  conduction  band  electrons.  The  constant 
K,  which  has  the  dimension  of  a  field,  can  be  calculated  from 
Eqs,  (8-7)  and  (8-8).  This  simple  (.»**„.?  '.once  has  the  same 
functional  form  as  the  Townsend  c^e  fC  *  c  •  *»r  t  in  eas  breakdown.^2 
Baraff1^  calculated  a(E)  from  a  nonlinear  Boltzmann 
equation  by  assuming  that  v  ( t )  is  a  constant  independent  of 
energy.  He  found  that  Eq.  (8-9)  approximately  describes  the 
functional  dependence  of  a(E)  at  low  field3  but  that  K  is 
actually  higher  than  Shockley  had  calculates  The  physical 
reason  for  this  descrepancy  is  that  Shockley  implicitly  assumed 
that  each  electron  ccllidinc  at  intermediate  energies  returns 
to  zero  energy  before  attemptine  another  flieht  to  hieh 
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energy.*^  This  assumption  leads  to  an  underestimate  of  the 
number  of  electrons  which  reach  ionizing  energies.  In  high 
fields  the  Ionization  coefficient  was  found  by  Baraff  to  oe 
approximately  proportional  to  exp  (-k’/E^).  Baraff s  cal¬ 
culations  were  Intended  to  describe  avalanche  multiplication 
in  semiconductors  such  as  Ge  although  his  choice  of  an  energy- 
independent  v,  is  probably  more  suggestive  of  breakdown  in 
gases. 

1  2 

Pass  and  Barrett  have  shown  that  Shockley's  model  can 
be  applied  to  breakdown  from  optical  electric  fields  by  con¬ 
sidering  a  sequence  of  " luckv"  collisions  that  reverse  the 
electron  momentum  at  about  the  instant  when  the  field  changes 
direction.  In  this  manner  a  “lucky"  electron  will  continually 
remain  in  phase  with  the  field  and  be  accelerated  to  ionizing 
energy.  A  probability  was  calculated  which  has  the  same 
functional  dependence  on  field  as  does  Eq.  (8-9).  Although 
this  simple  model  was  advanced  to  explain  laser  breakdown 
statistics  (see  Chapt.  6),  it  leads  to  a  prediction  of  the 
ionization  rate  a  (E)  in  the  same  manner  that  Shockley's 
calculation  of  P(tj)  leads  to  a(E).  This  model  suffers  from 
the  same  simplifications  as  Shockley ' s--the  assumption  of  a 
constant  vL  and  an  underestimate  of  the  number  of  electrons 
which  reach  ionizing  energy. 

Seitz  has  written  one  of  the  classic  papers  on  dc  electron 
avalanche  breakdown. His  description  of  the  avalanche 
contains  the  famous  "40  generations"  model  and  is  generally 
based  on  physical  arguments  rather  than  on  formal  mathematical 
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modeline.  This  work  contains  many  important  ideas  that  had 
not  previously  been  considered  in  breakdown  theory. 

Seitz  recoenized  that  the  LO  phonon  interaction  (or 
Frfthlich  interaction)  may  not  realistically  model  the  electron- 
lattice  coupling.  In  all  solids  the  electron  can  be  coupled 
to  the  lattice  by  a  deformation  potential  interaction.  This 
interaction  increases  with  increasing  energy  and  may  dominate 
the  LO  Dhonon  interaction  in  ionic  solids  when  *  >  = 

0.1  1.  A  perturbation  calculation  that  considers  both  types 
of  coupling  can  be  aoplied  to  calculate  an  effective  w{  ( < )  as 
shown  schematically  in  Fig.  31. 

The  "40  generations"  model  is  an  estimate  of  the  number 
of  orogeny  which  a  single  electron  must  produce  in  order  to 
create  a  damaging  avalanche.  With  simple  arguments,  Seitz 
demonstrated  that  if  a  single  electron,  starting  its  flight 
at  the  cathode,  produces  40  secondary  electrons  before  the 
avalanche  is  terminated  at.  the  anode,  then  enough  energy  will 
be  deposited  from  joule  heating  into  the  medium  near  the  anode 
to  destroy  the  material.  This  description  leads  naturally  to 
the  prediction,  later  verified, that  the  breakdown  field 
will  depend  on  sample  thickness.  If  the  specimen  thickness  is 
reduced,  the  ionization  rate  oer  unit  distance  must  increase 
in  order  to  compensate  for  the  reduced  spatial  extent  (or  time 
duration)  available  for  40  secondaries  to  be  produced.  In 
Chaot.  6  we  in  effect  considered  the  basic  concept  of  the 
"40  generations"  model  by  defining  breakdown  as  the  growth 
of  the  multi oil  cat ion  factor  Mc  to  a  critical  size.  (See 
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Fig.  31.  SCHEMATIC  OF  THE  ELECTRON-PHONON  COLLISION 

FREQUENCY  AS  A  FUNCTION  Or  ENERGY  FOR  BOTH 
POLAR  AND  NONPOLAR  INTERACTIONS  (Seitz  1949) 
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Eq.  (5-2).)  Our  choice  of  Mc  =  107  corresponds  to  23  genera- 

Qii 

tions. 

The  damage  field  for  NaCl  predicted  by  the  "40  generations** 
model  is  about  an  order  of  magnitude  larger  than  experimentally 
observed.  Seitz  explained  this  descrepanc.y  by  showing  that 
avalanche  breakdown  may  not  be  an  average  electron  phenomenon. 
Fluctuations  in  electron  energies  can  substantially  lower  the 
breakdown  threshold.  This  same  result  was  demonstrated  several 
years  later  in  the  calculations  of  Baraff^”  and  Holway. 

Seitz  considered  the  effects  of  exciton  states  on  break¬ 
down,  a  concept  introduced  in  Sect.  B  of  Chapt.  5.  He  also 
attempted  to  explain  the  origin  of  noisy  dc  pre-breakdown 
currents  that  had  been  observed  experimentally,  and  he  con¬ 
sidered  the  Dossibility  that  Bragg  reflection  at  the  Brillouin 
zone  boundary  (Unklapu  processes)  might  affect  the  breakdown 
strength  of  some  materials. 

67 

A  recent  model  by  Holway,  developed  by  analogy  to  gas 
breakdown,  orovides  a  relatively  simple  mathematical  formalism 
for  evaluating  many  theoretical  ideas  suggested  by  Seitz's 
work.  This  model,  which  will  be  considered  again  in  the  next 
section,  can  be  aoplied  at  dc  and  at  optical  frequencies. 

Holway  calculated  a  high-field  electron  distribution  function 
by  using  a  classical  Fokker-Planck  equation  and  a  relaxation¬ 
time  model  for  the  electron-phonon  interaction.  The  calculation 
used  Eq.  (8-5)  for  A(<)  and  Eq .  (8-3)  for  B  (O,  and  only  the 
high  energy  form  for  the  collision  frequency  in  Eq .  (8-2). 

Two  important  results  were  obtained.  First,  the  growth  of 
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the  electron  density  was  found  to  proceed  in  the  manner  in¬ 
dicated  by  Fq .  (5-2)  so  that  an  ionization  rate  could  be 
calculated.  Second,  the  term  in  the  Fokker-Planck  equation 
describing  the  average  electron  enerey  was  found  to  be  unim¬ 
portant  to  the  avalanche  dynamics.  As  Seitz  had  suggested^ 
and  as  Baraff  had  found  with  a  more  comdicated  Boltzmann 
equation  calculation,  the  avalanche  appears  to  develop  be¬ 
cause  of  electron  diffusion  to  higher  energies  and  not  be¬ 
cause  the  average  electron  has  been  heated  sufficiently  to 
produce  ionization. 

Finally,  we  mention  two  more  calculations  related  to 

CO 

breakdown.  In  the  first,  Thornber  and  Feynman  treated  the 
FrBhlich  interaction  quantum  mechanically  to  a  high  degree  of 
precision  over  the  energy  extending  to  about  twice  the  LO 
phonon  energy.  Their  results,  which  are  quite  involved 
mathematically,  are  not  directly  applicable  to  avalanche 
breakdown.  As  Seitz  has  ooserved,  the  LO  interaction  described 
approximately  by  the  FrBhlich  Hamiltonian,  may  not  completely 
describe  the  electron-ohonon  coupling  over  the  energy  range 
of  Interest.  We  also  realize  that  avalanche  breakdown  is  a 
statistical  orocess,  and  it  aooears  to  be  necessary  to  cal¬ 
culate  an  electron  distribution  function  before  any  realistic 
predictions  of  breakdown  strength  can  be  made.  The  details  of 
the  electron-ohonon  interaction  near  the  LC  resonance  may  not 
even  be  important  to  an  understanding  of  the  avalanche  develop¬ 
ment.  Since  the  cross-section  for  impact  ionization  will  not 
be  large  until  <  =1.2  I,  the  Droducts  of  the  ionization 
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event  may  have  energies  much  greater  thanti^Q . 

The  last  theory  of  breakdown  we  will  consider  is  Zener 

tunneling  which  at  high  frequencies  becomes  multiphoton  ab- 
1  8 

sorption.  Even  though  this  Drocess  is  not  an  electron 

avalanche,  it  can  cause  breakdown,  and  we,  therefore,  consider 

20 

it  briefly  here.  Zener  originally  introduced  the  concent  of 
field  emission  breakdown,  and  Franz^  nut  the  idea  on  a  more 
quantitative  basis  for  dc  fields.  An  aoproximation  to  Franz' 

t-Q 

results,  given  by  O'Dwyer,  predicts  a  breakdown  field  in 
volts/cm  of 


E 
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(8-10) 


where  I  is  the  band-gap  (ionization)  energy  in  eV  and  tp  is 
the  duration  of  the  applied  field  in  microseconds.  Eq .  (8-10) 
is  Probably  a  useful  approximation  to  optical  tunneling  field 
at  1  gm  in  NaCl,  For  tQ  =  10  ns,  Eq.  (8-10)  predicts 
Ec  =  2  x  107  volts/cm,  and  for  tc  =  15  ps,  Ec  =  5  x  107  volt3/cm. 
These  values  compare  to  0.2  x  10?  volts/cm  and  1.2  x  10? 
volts/cm  measured  at  1 . 06  pm  in  Chapt.  <.  This  process  thus 
appears  to  have  a  weaker  pulse-width  dependence  than  optical 
avalanche  breakdown  and  should  therefore  be  observed  when 
extremely  short  laser  pulses  (t  -  10“12  sec.)  are  used  to 
induce  damage. 


C.  Implications  of  Experimental  Results  to  Theory 


In  Chapts.  4  and 


C 

J 


we  observed  that  the  avalanche  break- 
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down  process  has  a  time  and  frequency  dependence  and  that  it 
apnears  to  be  affected  by  severe  crystalline  disorder.  We 
were  also  able  to  determine  accurate  values  of  breakdown  fields 
which  were  not  confused  by  the  many  experimental  uncertainties 
of  dc  measurements. 

The  observation  of  a  time  dependence  to  damage  suggests 
that  the  breakdown  models  which  we  called  "intrinsic  breakdown 
theories”  in  the  last  section  are  not  applicable  to  the 
materials  we  have  studied.  It  appears,  in  fact,  that  the 
distinction  between  intrinsic  breakdown  and  avalanche  break¬ 
down  is  entirely  artificial.  Breakdown  in  materials  must  have 
a  time  dependent  aspect. 

Our  data  on  the  time  dependence  to  the  damage  field  can 
be  compared  to  the  simple  prediction  o^  Rq .  (8-9)  usine  the 
definition  of  a  eiven  in  Eq,  (5-3).  This  comparison  is  made 
in  Fig.  32  where  it  is  seen  that  the  functional  form  of 
Shockley’s  model  is  not  an  unreasonable  approximation  to  the 
experimental  data. 

A  weak  frequency  dependence  to  the  damage  field  was 
observed  in  the  alkali  halides  in  Chapt.  5.  In  Sect.  F  of 
Chapt.  2,  a  comparison  between  the  sapphire  damage  field  at 

pi 

0.69  gm  measured  indirectly  by  Cuiliano  et  al.  1  and  our  own 
measurement  at  1.06  gm  (Table  X  in  Chrpt.  7)  indicates  that 
the  sapphire  damage  field  is  still  effectively  in  its  dc 
limit  at  ruby  frequency.  These  surprising  results  suggest 
that  the  effective  electron-ohonon  collision  times  in  these 
materials  are  of  order  10"10  sec,  a  value  so  short  that  the 
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Fig.  32.  COMPARISON  OF  SUBNANOSECOND  LASER 
DATA  TO  SHOCKLEY  MODEL 
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electron  energy  uncertainty  is  comparable  to  the  band-gap 
energy.  If  the  effective  collision  frequency  v  is  determined 
by  the  electron  dynamics  at  energies  much  larger  than^^Q, 
as  Holway  has  found,  then  the  Frbhlich  interaction  is  not 
capable  of  explaining  this  result.  We  are  left  with  the 
suggestion  of  Seitz  that  the  non-polar  or  deformation  potential 
interaction  may  be  important  in  describing  avalanche  breakdown. 
In  addition,  the  details  of  the  frequency  dispersion  summarized 
in  Fig.  19  of  C hapt.  5  give  evidence  that  deep-lying  exclton 
levels  may  play  a  role  in  the  avalanche  at  least  for  wide 
band-gap  Insulators  such  as  NaF. 

We  thus  realize  that  the  high  frequency  limit  of  the 
Fr&hlich  interaction  which  was  used  by  Zverev  et  al. . 

ffrj  6*7  6  A 

Wasserman,  •*  and  Holway  '  '  does  not  describe  the  energy 
gain  in  intense  optical  fields  at  least  when  the  frequencies 
are  less  than  about  4  x  10**4,  herz  in  the  materials  we  have 
studied.  It  should  be  noted,  incidently,  that  the  laser  fre¬ 
quency  dependence  displayed  by  Eq.  (1-1)  is  not  the  true 
prediction  of  rnlaxation-time  models  of  avalanche  breakdown. 
That  prediction  can  only  be  obtained  by  averaging  Eq .  (1-1) 
over  the  appropriate  electron  distribution  function. 

The  dependence  of  damage  field  on  severe  crystalline 
disorder  was  considered  in  Chapt.  5,  and  the  observations  of 
breakdown  statistics  were  summarized  in  Chapt.  6,  Here  we 
merely  note  that  these  results  further  support  the  conclusion 
that  the  intrinsic  breakdown  mechanism  involves  ionization  from 
electron  heating  and  not  from  tunnel  ionization  or  multiphoton 
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absorption. 

Virtually  all  breakdown  theories  predict  damage  fields 
about  an  order  of  magnitude  larger  than  we  have  observed.  An 
important  reason  for  this  discrepancy  appears  to  be  that 
most  models  have  treated  avalanche  breakdown  as  an  average 
electron  phenomenon.  It  can  be  seen  very  simply,  however, 
that  the  average  electron  energy  is  probably  not  much  greater 
than  about when  damage  develops.  Consider  NaCl,  for 
example.  The  low-field  mobility  of  about  10  cm  /volt-sec  can 
be  treated  as  an  upper  bound  for  the  high-field  mobility.  If 
we  u3e  the  free  electron  mass,  a  field  of  2  x  10^  volts/cm 
produces  an  average  electron  energy  of  less  than  about  0.2  eV 
as  compared  to  =  kT  =  0.03  eV  and  a  band-gap  energy  of 

about  8  ev.  Because  the  mobility  should  decrease  in  high 
fields,  ;  the  actual  average  electron  energy  may  be  much  less 
than  0.2  eV.  As  Baraff  and  Holway  have  shown,  the  diffusion 
of  electrons  towards  the  high  energy  tail  of  the  energy  dis¬ 
tribution  thus  appears  to  be  more  important  than  the  average 
electron  energy  in  determining  the  avalanche  characteristics. 
Including  the  diffusion  term  in  the  calculation  lowers  the 

68 

predicted  damage  field  by  about  a  factor  of  three  or  more. 

Another  source  of  the  discrepancy  may  be  in  the  choice 
of  modeling  for  v(<).  Based  on  Seitz's  discussion  as  sum¬ 
marized  in  Fig.  31,  it  appears  that  a  constant  vL  may  be  a  more 
realistic  approximation  to  reality  than  the  perturbation 
modeling  of  the  Frbhlich  interaction.  It  may  also  be  that 
the  classical  relaxation-time  approximation  normally  used  to 
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predict  the  electron  dynamics  may  not  be  capable  of  describing 
avalanche  breakdown. 

A  comDlete  modeling  of  avalanche  breakdown  alone  the  lines 
of  the  Thornber-Peynman  calculation  is  probably  many  years 
distant.  Nonetheless  we  believe  that  it  is  possible  to  under¬ 
stand  many  qualitative  features  of  the  avalanche  by  using 
classical  models  such  as  that  developed  by  Holway.  In  closing 
this  section  we  will  consider  possible  extensions  of  Holway's 
calculations  which  may  be  useful  in  answering  a  number  of 
theoretical  questions  raised  in  the  preceding  discussion. 

The  Fokker-Planck  equation  approach  appears  to  be  equiva¬ 
lent  to  using  a  first-order  nonlinear  Boltzmann  equation.9^ 

The  advantages  it  provides  over  the  Boltzmann  approach  are 
that  the  Fokker-Planck  equation  is  simpler  from  a  computational 
viewpoint  as  well  as  easier  to  interpret  physically. 

By  using  Holway’s  formalism,  a  number  of  problems  can  be 
treated  in  a  direct  manner.  It  is  possible,  for  example,  to 
determine  the  importance  of  the  details  of  the  electron- 
phonon  interaction  to  the  development  of  the  avalanche.  Seitz's 
suggestion  that  the  non-polar  interaction  may  be  important  can 
be  tested  within  this  formalism  by  choosing  different 
functional  forms  for  v(0,  including  one  that  takes  account  of 
both  the  polar  and  the  non-polar  interactions. 

A  second  problem  to  be  resolvedi  Is  the  relaxation  time 
approximation  capable  of  explaining  our  experimental  results 
in  the  alkali  halides?  In  particular,  can  Holway's  formalism 
with  a  reasonable  choice  for  ^(O  predict  the  relative  damage 
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fields  we  have  measured  in  Chapts.  4  and  5*  and  can  it  predict 
a(E)  in  NaCl?  It  may  be  possible  to  fit  i'(<)  to  the  experi¬ 
mental  results. 

Seitz's  suggestion  that  deep-lying  exciton  states  may  be 
important  to  breakdown  can  be  tested  by  this  model,  and  it 
may' be  possible  to  develop  generalized  predictions  for  a(E) 
based  on  material  parameters.  The  role  of  impurity  states 
can  also  be  investigated . 

In  conclusion,  we  see  that  although  the  purist  may  find 
avalanche  theory  to  be  in  a  dismal  state,  relatively  simple 
phenomenological  models  may  be  caoable  of  explaining  many 
features  of  avalanche  breakdown.  In  this  area  of  work, 
however,  experiment  must,  lead  the  theory,  and  much  more 
experimental  data  is  needed  to  Provide  both  a  basis  and  a 
test  for  theoretical  models. 


CHAPTER  9 

SUGGESTIONS  FOR  FUTURE  WORK 


The  study  of  laser  induced  damage  that  has  been  presented  in  the 
preceding  chapters  has  considered  damage  from  absorbing  inclusions, 
beam  distortion  from  self-focusing,  and  electron  avalanche  breakdown. 

Many  of  the  basic  problems  concerning  laser  damage  have  been  resolved 
by  this  work  and  by  other  recent  studies,  but  a  great  many  problems 
remain.  In  this  concluding  section  we  will  briefly  consider  some  areas 
for  future  experimental  work.  The  following  list  is  not  intended  to  be 
exhaustive. 

1 .  High-Frequency  Avalanche  Breakdown  and  Multiphoton 

Absorption:  The  first  evidence  of  frequency  dispersion  in  the  damage 
fields  of  the  alkali  halides  has  been  obtained.  Because  this 
dispersion  does  not  follow  simple  predictions,  measurements  of 
high  frequency  avalanche  breakdown  should  be  extended  to  wave¬ 
lengths  shorter  than  0.  69  pm.  It  would  be  of  great  interest  to  track 
the  frequency  dependence  of  intrinsic  breakdown  over  a  sufficient 
range  of  frequencies  to  show  the  eventual  development  of  multiphoton 

absorption  as  the  intrinsic  damage  mechanism.  Except  possibly  for 

1  3 

the  recent  work  of  Bass  and  Barrett,  no  unambiguous  measure¬ 
ments  have  yet  been  made  of  multiphoton  absorption  and  its  frequency 
and  pulse-width  dependences.  As  the  laser  frequency  is  increased, 
complications  from  self-focusing  may  get  more  serious  and  tighter 
external  focusing  may  be  required.  The  effective  limit  to  external 
focusing  will  be  dtermined  by  spherical  aberrations.  Laser 
damage  studies  at  frequencies  near  the  two-photon  absorption 
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resonance,  therefore,  may  he  more  profitably  conducted 
on  surfaces  rather  than  in  the  bulk,  orovided  that 
essentially  imoerfect ion- free  surfaces  can  be  produced. 
Initial  work  in  this  area  might  be  conducted  in  the  alkali 

i 

halides  with  a  nitrogen- laser-Pumoed  dye  l*aser.  It  would 

♦' 

be  of  interest  to  compare  the  frequency  dependence  of 
intrinsic  damage  in  crystalline  quartz  to  that  of  fused 
quartz. 

Subnanosecond  Damage  Measurements t  In  Chapt.  6  it  was 

.  » 

» 

shown  that  an  approximate  ionization  rate  can  be  found 
from  the  dependence  of  the  intrinsic  damage  field  on  laser 
oulse  width.  These  techniques  can  be  applied  to  other 
materials  iri  order  to  determine  if  the  ionization  rate  has 
a  universal  form  in  insulators.  The  frequency  dependence 
of  the  ionization  rate  can  also  be  determined  from  sub¬ 
nanosecond  studies.  It  should  be  noted  that  dc  measure¬ 
ments  of  dielectric  breakdown  cannot  directly  duplicate 
these  studies  because  high  voltage  quasi-dc  impulses  of 
picosecond  duration  have  not  been  achieved.  In  addition, 
dc  measurements  are  difficult  to  conduct  and  may  be 
plagued  by  a  number  of  poorly  understood  experimental 
complications--surface  and  electrode  effects,  space 
charges,  internal  strains,  etc. 

Effect  of  Disorder  on  Avalanche  Breakdown!  Measurements 
were  Dresented  in  Chaot.  6  which  suggest  that  severe 
microscopic  disorder  can  increase  the  breakdown  strength 
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of  solids.  This  work  can  be  extended  to  other  trans¬ 
parent  materials  which  can  exist  in  both  crystalline  and 
amorphous  forms  and  to  alloys  such  as  YAG  (yttrium 
aluminum  garnet) -YGG  (yttrium  gallium  garnet)  which  can 
be  oreoared  in  arbitrary  proportions. 

Exciton  Fluorescence!  Evidence  has  been  obtained  which 

shows  that  excitons  can  be  excited  in  insulators  by  high- 

9? 

oower  laser  beams  confined  to  subd amarine;  intensities. 
Prof.  H.  Ehrenreich  has  suggested  to  the  author  that  a 
study  of  exciton  fluorescence  during  damage  may  help 
clarify  the  role  of  deep  lying  exciton  states  in  the 
avalanche  Drocess.  Such  work  is  only  useful  in  large 
band-gap  materials  such  as  NaP,  because  the  smaller  band- 
gap  alkali  halides  have  exciton  states  which  are  lifetime- 
broadened  sufficiently  in  damaging  electric  fields  to 
overlap  the  conduction  band.  Vacuum  ultraviolet  tech¬ 
niques  would  therefore  be  required  for  such  studies. 

Temperature  Dependence  of  Avalanche  Hreakdownt  A  great 
deal  of  work  has  been  published  on  the  temperature  de¬ 
pendence  of  dc  avalanche  breakdown.  Much  of  this  work, 
however,  is  confused  and  inconsistent.  Lasers  provide  a 
useful  means  for  measuring  the  temperature  dependence  of 
intrinsic  damage  thresholds.  A  careful  study  of  the 
temperature  dependence  of  intrinsic  breakdown  has 
important  implications  to  avalanche  theory. ^ 


Effects  of  Impurities  on  Intrinsic  Damage  Fieldst  In 
ChaDt,  1  we  discussed  the  role  of  imourity  absorption  and 
found  that  unless  the  imourity  electrons  can  be  excited 
into  the  conduction  band  (or  imourity  holes  into  the 
valence  band),  impurities  will  have  little  effect  on 
damage  fields.  Wide  band-cap  insulators  are  thus  not 
particularly  attractive  for  impurity  studies.  Semicon¬ 
ductors  such  as  OaAs  or  Ge,  however,  can  be  used  with  a 
CO2  laser  to  study  impurity  effects  in  optical  damage. 

One  can  imagine  strongly  doping  the  semiconductor,  and 
measuring  the  intrinsic  damage  threshold  as  a  function  of 
imourity  concentration.  If  shallow  thermalized  states  are 
introduced,  the  gain  coefficient  in  Eq .  (5-2)  can  in 
principle  be  varied  and  directly  measured  as  function  of 
field . 

Damage  to  Thin  Film  Orticsi  In  the  present  work  there  has 
been  no  discussion  of  thin  film  optical  damage.  Very 
little  fundamental  information  about  optical  damage  can 
be  obtained  from  thin  film  studies,  but  thin  film  coatings 
are  nonetneless  quite  imoortant  to  the  practical  design 
of  laser  systems.  Because  we  now  understand  the  role  of 
contaminants  and  structural  imperfections  in  the  damage 
Drocess  and  because  we  now  have  careful  measurements  of 
intrinsic  bulk  damage  fields  to  nerve  as  standards,  the 
comD'licated  problem  of  thin  film  damage  can  be  approached 
with  some  confidence.  Thin  film  preparation  techniques 
are  needed  which  minimize  structural  defects  and 


-191- 


contaminatlon,  and  diagnostic  tools  must  be  developed  to 

monitor  the  quality  of  optical  coatings.  As  Dr.  J.  M.  Khan 
fi  1 

has  observed.  J  thin  film  damage  is  fundamentally  a 
materials  problem  which  probably  requires  sophistocated 
vacuum  and  surface  equipment. 

8,  Gas  and  Liquid  Breakdown!  Laser  induced  gas  breakdown 
can  be  rrofitablv  studied  with  the  techniques  of  the 
present  work.  Because  the  pressure  of  eases  is  an  experi¬ 
mental  variable,  it  is  relatively  easy  to  study  the  com¬ 
petition  between  multiphoton  ionization  and  avalanche 

breakdown.  Work  along  these  lines  has  been  reported  in 

o 

the  literature  with  fixed  laser  pulse  durations,-7  but 
the  pulse-width  dependence  of  intrinsic  processes  has  not 
yet  been  carefully  determined  and  interpreted  in  terms  of 

Qfl 

dc  breakdown.  Gas  breakdown  studies  are  particularly 
attractive  because  the  breakdown  process  is  much  better 
understood  in  gases  at  dc  and  at  optical  frequencies  than 
it  is  in  solids.  Liquids  also  offer  advantages  with 
respect  to  laser  damage  studies.  Atomic  impurities  with 
short  excited  state  relaxation  times  can  be  added  to  a 
liquid  so  that  the  effects  of  absorption  on  laser  damage 
thresholds  can  be  determined. 

Optical  damage  is  thus  seen  to  be  an  area  where  much 
fundamental  and  practical  work  remains.  It  is  hoped  that  the 
present  study  is  not  viewed  as  a  closing  of  another  door  of 
experimental  physics  but  rather  that  it  is  viewed  as  the 
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first  series  of  stens  in  solving  a  large  and  Comdex  set 
of  oroblems. 


APPENDIX 


A  practical  limit  to  the  decree  of  focusing  that  can  be 
achieved  with  any  real  lens  is  set  by  spherical  aberrations. 
This  distortion  comes  about  because  the  position  of  the  focus 
for  light  rays  at  varying  distances  above  the  axis  of  the  lens 
deDends  on  the  value  of  that  distance.  A  perfectly  collimated 
light  beam,  then,  which  would  in  principle  be  focused  to  a 
ooint  at  the  geometrical  focus,  produces  instead  a  blurred 
image  with  a  finite  radius.  Spherical  aberrations  will  also 
be  introduced  when  a  light  beam  is  focused  through  any  plane 
surface. 

In  this  appendix  we  will  show  that  the  spherical  aber¬ 
rations  introduced  by  the  focusing  lens  and  plane  front  sur¬ 
face  of  the  samples  under  study  in  the  bulk  damage  studies 
were  negligible  comoared  to  diffraction  effects.  Aberrations 
of  the  lens  will  be  characterized  by  using  manufacturer’s 
SDecif ications  on  the  lens  along  with  well-known  results 
available  in  most  standard  optics  textbooks.  A  formula 
describing  the  aberrations  introduced  by  focusing  through  a 
single  plane  surface  will  be  derived  from  Snell’s  law. 

When  a  focusing  lens  is  aligned  symmetrically  with  respect 
to  a  diffraction-limited  light  beam,  the  intensity  distribu¬ 
tion  at  the  focus  is  the  far-field  diffraction  pattern  of  the 
beam  at  the  entrance  to  the  lens  convoluted  with  a  function 
describing  the  lens'  spherical  aberrations.^  For  any  par¬ 
ticular  lens,  the  effective  focal  diameter  from  spherical 
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aberrations  is  approximately  proportional  to  R ^ ,  where  R  is 
the  incident  beam  radius,  and  the  diffraction-limited  focal 
diameter  is  -(A/2R)f  where  f  is  the  focal  length  of  the  lens. 
The  relative  size  of  the  spherical  aberrations  for  a  particular 
lens  design  as  compared  to  diffraction  is  thus  driven  by 

Spherical  Aberrations  d4a  it  1 1 

Diffraction  Diameter  '  ' 

and  it  is  seen  that  spherical  aberrations  become  rapidly  more 
important  as  the  incident  beam  diameter  is  increased.  An 
optimum  value  of  R  exists  which  produces  a  minimum  effective 
focal  diameter.  For  any  larger  value  of  R,  the  spot  size 
increases  and  spherical  aberrations  limit  the  focal  diameter. 
For  smaller  values  the  spot  size  also  increases,  but  now 
diffraction  limits  the  focal  diameter. 

The  focusing  lens  used  in  the  damage  studies  was  designed 
for  minimum  spherical  aberrations  by  properly  selecting  the 
radii  of  curvature  for  the  two  surfaces.  According  to  the 
manufacturer's  specifications,^00  a  0.633  /nn  diffraction- 
limited  laser  beam  incident  on  this  lens  will  produce  the 
smallest  possible  focal  diameter  when  the  incident  beam  has  a 
diameter  of  3  mm.  It  can  be  assumed,  therefore,  that  spherical 
aberrations  and  diffraction  effects  are  equal  when  R  =  1 . 5  mm 
and  A  =  0,633  fun, 

In  our  experiments  at  1.0 6  (.a n,  we  used  a  beam  with  a 
diameter  of  just  under  1  mm  so  that  R  0.  5  mm  and  A  =  1.06  n m. 
The  relative  size  of  the  spherical  aberrations  as  compared  to 
diffraction  can  then  he  calculated  from  Eq .  (A-l)  for  our 
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conditions  of  measurement.  It  is  found  to  be  about  0.01. 
Because  the  true  beam  shape  is  a  convolution  of  the  aberration 
function  and  the  diffraction  pattern,  the  size  of  the 
spherical  aberrations  and  the  diffraction-limited  diameter  do 
not  add  to  produce  the  true  beam  size  at  the  focus.  Rather, 
they  combine  approximately  as  the  square  root  of  the  sum  of 
their  respective  squares.  Spherical  aberrations,  therefore, 
introduce  a  necliaible  error  of  about  0.01??  into  the  measured 
breakdown  intensities. 

Spherical  aberrations  introduced  by  a  plane  surface  can 
be  important  for  tierhtly  focused  laser  beams.  Fig.  33  illus¬ 
trates  the  source  of  aberrations.  An  incoming  light  ray, 
which  intercepts  the  entrance  plane  at  a  height  h  and  forms 
an  angle  9  with  respect  to  the  horizontal,  focuses  at  a 
distance  d2  behind  the  entrance  plane.  It  will  be  shown  that 
d 2  depends  on  h  so  that  different  rays  will  focus  at  different 
points  in  the  medium,  and  the  focal  area  will  consequently 
be  blurred. 

From  Snell’s  law  we  can  relate  $'  to  $.  In  particular. 


sin  9 


n' 

n 


sin  A' 


N  sin  c 


(A-2) 


where 

sin  G  -  •  -  h-  -  ,  sin  9-  ,  - [2 - 

.:.Z  2  ~2 - 2 

vh  +dl  'h  +d2 

Sq.  (A-2)  can  be  used  to  find  d2  in  terms  of  d^.  With  the 
substitution  h/dj_  =  tan  9,  we  find  that 


-197- 


d2(N  =  Nd, 


tan  0  +  1 


(A-3) 


The  longitudinal  spherical  aberration  (LSA)  is,  by 
definition,  d2(h)  -  d2(0).  According  to  Eq,  (A-3), 


LSA  --  N 


t'i|  (^x) tan2  6  +  1 


The  quantity  of  interest  is  the  transverse  spherical  aber¬ 
ration  (TSA)  which  is  equal  to  the  radius  of  the  blurred 
focal  spot.  Its  value  is  LSA  x  tan  Q'  or,  since 


tan 


sin  & 


sin  6  1 


cos  0  ( 1  - s in  3  ' ) 


_  si;'  & _ 

i sinks') ^  ■  2~ 


the  radius  of  the  blurred  focal  spot  is 

/n2_m  2  1/2  | 

TSA  -  Ndj  tan  6  r  1  -11- 


X2  .  2  1/2 
(N  -sin  0 >  ' 


(A-4) 


Under  the  conditions  of  focusing  used  in  the  bulk  damage 
studies,  tan  6  -  sin  d  -  0.  Eq .  (A-4)  then  becomes 


TSA  »  d 


(A-5) 


If  the  plane  entrance  face  is  in  the  far  field  of  the 
focused  laser  beam,  then  <9~-A/2aQ  where  aQ  is  the  minimum 
beam  radius  inside  the  material.  We  have  focused  approximately 
2  mm  inside  the  samples,  so  that  d^  -  1.3  mm.  From  Eq.  (A-5) 


4 
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the  calculated  value  of  TSA  is  less  than  0.1  ,u m  as  compared  to 
a  diffraction-limited  beam  waist  aQ  of  about  10  /tm  for  th'i 
YAGiNd  Q-switched  laser  experiments.  We  thus  see  that,  as' 
was  the  case  with  the  focusine  lens,  the  distortion  intro¬ 
duced  by  the  plane  entrance  face  is  negligible  for  our 
conditions  of  measurement 
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